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Share  And  Modtran  Merged 

SHARC 

Strategic  High-altitude  Atmospheric  Radiation  Code 

SWIR 

Short  Wave  Infrared 

SPIRE 

Spectral  Infrared  Rocket  Experiment 

SSGM 

Strategic  Scene  Generation  Model 

SUSIM 

Solar  Ultraviolet  Spectral  Irradiance  Monitor 

SSI 

Spectral  Sciences,  Inc. 

SSBUV 

Shuttle  Solar  Backscatter  Ultraviolet 

SZA 

Solar  Zenith  Angle 

UARS 

Upper  Atmosphere  Research  Satellite 

uss 

User  Supplied  Spectra 

uv 

Ultraviolet 

VNBM 

Very  Narrow  Band  Model 

WBM 

Wide  Band  Model 

1-D 

One  Dimensional 
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1.0  INTRODUCTION 


The  calculation  of  infrared  (IR)  radiance  and  transmittance  spectra  is  of  fundamental  importance 
in  many  areas  of  atmospheric  science.  These  include  modeling  the  atmospheric  energy  budget  for 
global  climate  change  predictions,  analyzing  data  from  remote  sounding  experiments,  and 
understanding  the  impact  of  the  atmosphere  on  target  detection.  The  usefulness  of  the  models  and 
codes  used  in  simulations  lies  in  the  fundamental  understanding  of  the  underlying  physics  and 
chemistry  that  give  rise  to  these  atmospheric  effects.  The  overall  objective  of  this  effort  is  to 
evaluate  our  understanding  of  IR  atmospheric  phenomenology  through  comprehensive  analysis  of 
spectral  and  in-band  radiance  field  data,  and  to  transition  these  results  into  practical  models  for 
inclusion  into  existing  atmospheric  simulation  codes. 

Available  computer  codes  for  atmospheric  IR  radiation  that  have  gained  wide  use  include  several 
developed  and  maintained  by  AFRL/VS.  MODTRAN  [Berk  et  ah,  1998]  is  a  rapid,  local 
thermodynamic  equilibrium  (LTE)  code  for  low-  to  moderate-resolution  spectra.  FASCODE 
[Clough  et  ai,  1988]  is  a  high-resolution,  line-by-line  code  that  can  be  used  in  both  LTE  and 
non-LTE  applications;  however,  the  non-LTE  molecular  state  populations  must  be  externally 
generated.  SHARC,  the  Strategic  High-altitude  Atmospheric  Radiation  Code  [Sundberg  et  al,  1 995] 
calculates  emission  and  transmittance  spectra  in  the  2-40  pm  region  for  arbitrary  line-of-sight  (LOS) 
paths  between  30  and  300  km  with  sufficient  accuracy,  speed,  and  resolution  (1.0  cm  ')  to  be  useful 
for  a  wide  range  of  applications.  It  incorporates  the  IR  bands  of  NO,  C02, 03,  H20,  OH,  CO,  and 
CH4  found  in  the  quiescent  atmosphere,  including  minor  isotopic  bands  of  C02  and  H20.  It  also 
accounts  for  auroral  production  and  excitation  of  C02,  NO,  and  NO+  caused  by  the  flux  of  energetic 
solar  electrons.  SAMM  [Sharma  et  ai,  1 996],  SHARC  and  MODTRAN  Merged,  was  developed  by 
integrating  MODTRAN  and  SHARC  into  a  single,  seamless  code  with  a  unified  and  correlated 
radiative  transport  algorithm  applicable  at  altitudes  from  ground  to  300  km. 

This  Final  Report  summarizes  the  work  performed  by  Spectral  Sciences,  Inc.  (SSI)  under  the 
Infrared  Airglow  and  Contamination  Modeling  and  Interpretation  program.  The  overall  objectives  of 
the  program  are  to  interpret  and  model  recent  infrared  airglow  observations  using  first-principles 
radiation  codes.  The  main  tasks  of  this  effort  include:  (1)  analyzing  data  from  the  CIRRIS-IA 
Shuttle  experiment  through  comparisons  with  calculations  using  the  SHARC,  SAMM  and/or 
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MODTRAN  codes,  (2)  providing  MODTRAN  upgrades  which  address  the  data/model  comparisons, 
and  (3)  investigating  the  accuracy  of  MODTRAN  in  selected  radiation  transfer  applications. 

The  report  is  organized  as  follows:  Section  2  presents  a  summary  of  a  detailed  spectral  and 
radiometric  analysis  of  the  CIRRIS-IA  data  using  SAMM;  section  3  describes  the  addition  of  a 
correlated-k  capability  to  MODTRAN;  a  summary  of  Band  Model  Calculations  of  Atmospheric 
Fluxes  and  Cooling  Rates  is  given  in  Section  4;  the  creation  of  a  MODTRAN  input  file  from 
radiosonde  measurements  is  described  in  Section  5;  Section  6  discusses  improvements  to 
MODTRAN  accuracy  for  remote  sensing;  the  new  solar  irradiance  data  sets  incorporated  in 
MODTRAN  is  presented  in  Section  7;  next,  the  incorporation  of  NOVAM  aerosol  model  into 
MODTRAN  3.7/4.0  is  described  in  Section  8;  finally,  upgrades  to  the  MODTRAN  3.7/4.0  aerosol 
and  cloud  models  are  discussed  in  Section  9. 
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2.0  SPECTRAL  AND  RADIOMETRIC  COMPARISONS  OF  CIRRIS-1 A 

DATA  AND  SAMM 

Atmospheric  infrared  (IR)  radiation  constitutes  the  background  against  which  targets  are  viewed 
by  many  space-based  optical  sensors.  The  operability  of  these  system  elements  is  limited  by  the 
ability  to  discriminate  target  signatures  from  the  atmospheric  optical  background.  The  cost-effective 
acquisition  of  IR  sensors  for  tactical,  theater,  and  strategic  weapon  systems  requires  accurate, 
validated  atmospheric  IR  radiance  models  for  generating  environment  simulations  over  a  wide  range 
of  wavelengths,  lines  of  sight,  and  atmospheric  conditions. 

The  highly  successful  April  1991  CIRRIS-IA  mission  [Ahmadjian  et  al,  1990]  on  the  space 
shuttle  STS-39  and  the  development  of  the  SAMM  atmospheric  optical  background  code  [Sharma  et 
al. ,  1 996]  provide  a  unique  opportunity  for  a  comprehensive  comparison  of  high-quality  spectral  and 
radiometric  data  with  code  calculations.  This  section  describes  a  systematic  comparison  of  CIRRIS 
atmospheric  IR  radiance  measurements  with  SAMM  code  calculations,  and  the  generation  of  a 
database  of  matched  overlay  plots.  The  CIRRIS  spectral  data  were  taken  with  a  Michelson 
interferometer  and  a  Si:As  detector  array  (focal  plane  3),  and  include  limb  views  for  tangent  heights 
ranging  from  near  the  ground  to  over  200  km  during  different  times  of  the  day  and  night  and  over  a 
wide  latitude  range  in  both  hemispheres.  Medium-wavelength-IR  spectra  were  also  taken  in 
near-nadir,  mostly  sunlit  viewing  geometries.  In  addition,  in-band  radiance  profiles  were  measured 
using  filtered  radiometers  (in  focal  planes  1  and  2)  for  eight  wavelength  bandpasses. 

The  CIRRIS-1  A-SAMM  database  contains  approximately  70  1-em  ' -resolution  spectra  and  70 
radiometer  profiles  measured  by  the  CIRRIS  experiment  for  non-auroral  atmospheric  conditions. 
The  SAMM  calculations  used  to  compare  with  these  data  are  based  on  24  model  atmospheres, 
corres|x>nding  to  eight  geographic  locations  along  the  shuttle  orbit  and  three  levels  of  geomagnetic 
activity.  The  calculations  properly  account  for  the  CIRRIS- 1 A  instrument  wavelength  response, 
resolution,  and  field  of  view.  The  CIRRIS- 1 A  data  were  evaluated  for  near-field  contamination  and 
other  experimental  artifacts.  The  agreement  between  the  calculated  and  observed  data  is  found  to  be 
generally  good,  typically  within  a  factor  of  two  for  in-band  radiances.  The  discrepancies  are  due 
mainly  to  the  assumed  profiles  of  temperature  and  certain  minor  species  and  to  the  current  lack  of 
mechanisms  in  SAMM  for  the  quiescent  daytime  chemiluminescence  from  NO  high-v,J  states  and 
from  NO\  which  are  important,  new  emission  process  discovered  in  the  CIRRIS-1  A  data. 
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This  section  discusses  the  overall  methodology  and  key  results  of  the  study.  A  more  detailed 
presentation  of  the  results,  including  basic  technical  documentation  on  the  database  files  and 
software,  and  complete  sets  of  radiometric  and  spectral  comparison  plots  are  contained  in  a  technical 
report  by  [Adler-Golden  et  al.  1994].  The  CIRRI S-IA-SAMM  database  is  available  from  the  Air 
Force  Research  Laboratory  (AFRLA^S).  It  is  planned  that  the  database  will  be  incorporated  into  the 
BMDO  Strategic  Scene  Generation  Model  (SSGM)  as  well  as  the  AFRLA^S  integrated  optical 
background  code  PLEXUS. 

2.1  SAMM  Calculations 

2.1.1  Code  Overview 

SSI  and  AFRL  developed  the  SAMM  code  [Sharma  et  al. ,  1 996]  for  self-consistently  calculating 
atmospheric  radiance  and  transmittance  from  0  to  300  km  for  the  2  to  40  pm  spectral  range.  The 
code  smoothly  integrates  the  high-altitude,  non-LTE  SHARC  [Sundberg  et  al,  1995]  algorithms 
with  the  MODTRAN  [Berk  et  al,  1998]  model  for  lower  altitudes.  SAMM  includes  a 
first-principles  vibrational  non-LTE  calculation  above  30  km,  auroral  and  terminator 
phenomenologies  with  multiple  atmospheric  regions,  aerosol,  cloud,  and  rain  models,  single  and 
multiple  scattering,  and  1  cm  '  spectral  binning  for  an  effective  resolution  of  about  2  cm‘.  The  input 
model  atmosphere  profile  is  selected  by  the  user.  The  spectra  and  radiance  profiles  in  the 
CIRRIS-SAMM  data  base  were  calculated  using  a  pre-released  Version  of  SAMM. 

2.1.2  Inputs 

For  the  CIRRIS-SAMM  database,  SAMM  calculations  were  performed  using  model  atmosphere 
profiles  generated  using  the  December  1993  version  of  the  SHARC/SAMM  Atmosphere  Generator 
(SAG)  [Adler-Golden,  1993].  SAG  incorporates  the  NASA  MSISE-90  empirical  model  [Hedin, 
1991]  for  temperature  and  major  species,  AFGL  [Anderson  et  al ,  1 986]  and  NRL  [Summers,  1 993] 
climatology  databases  for  minor  species,  and  additional  photochemical  and  empirical  models  for 
COj,  HiO,  odd  oxygen,  and  NO.  The  atmosphere  profiles  are  a  function  of  the  local  time,  latitude, 
Julian  day,  solar  flux  indices  FI 0.7  and  F10.7A,  and  the  geomagnetic  activity  index  Ap.  The  solar 
zenith  angle  is  a  function  of  the  local  time  and  latitude  and  is  input  separately  to  SAMM. 
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To  span  the  range  of  conditions  covered  during  the  CIRRIS  flight,  SAG  and  SAMM  were  run  for 
eight  geographic  conditions  and  three  Ap  values,  for  a  total  of  24  model  atmospheres.  The 
geographic  conditions  were  selected  to  evenly  sample  the  latitudes  and  solar  zenith  angles  along  the 
vehicle  orbit  (see  Figure  1).  For  each  atmosphere,  the  line  of  sight  (LOS)  spectral  radiances  from 
200  cm  '  to  5000  cm  '  (2  to  50  pm)  were  calculated  for  38  tangent  views  between  -10  km  and  240 
km  and  two  downward  views  (0°  and  30°to  the  nadir)  using  the  same  observer  location  and  LOS 
azimuth  angle.  Table  1  lists  the  SAMM  input  variables  pertaining  to  the  geographic  conditions  and 
the  LOS  geometries  at  a  tangent  height  of  1 00  km.  These  are  based  on  actual  CIRRIS  viewing 
conditions  during  Blocks  12A  and  16.  Table  2  lists  additional  input  parameters  required  for  SAG 
and  SAMM.  Initial  SAMM  calculations  showed  that  using  the  Mie  phase  function  for  aerosol 
scattering  gives  an  artificial  discontinuity  in  the  radiance  profile  at  30  km,  due  to  the  transition  in 
MODTRAN  between  the  stratospheric  and  thermospheric  aerosol  models.  For  subsequent 
calculations  the  Henyey-Greenstein  phase  function  with  the  default  asymmetry  parameter  value  of 
0.9  was  used,  and  these  results  are  presented  in  the  database. 


O  CIRRIS  100  KM  VIEWS 

m.  SAMM  CALCS . 


SOLAR  ZENITH  ANGLE  (DEG) 

Figure  1.  CIRRIS  Viewing  Geometries  at  100  km  Tangent  Height  During  Non-Auroral  Data 
Blocks 
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Table  1.  Geographic  and  LOS  Geometry  Inputs  for  SAMM  Calculations  (Tangent  Height  = 
100  km) 


Atmosphere 

Description 

Local 

Time 

(hr) 

Tangent  Pt 
SZA  (deg) 

Observer® 

Tangent  PT 

Sun 

GMT 

(hr) 

Lat 

(deg) 

Long 

(deg) 

LOS 

Az 

(deg) 

Lat 

(deg) 

Long 

(deg) 

MidLat  Night 

1.92 

149.9 

-32.3 

30.7 

-66.6 

-26.6 

17.7 

0.66 

LowLat  Night 

3.59 

122.9 

5.0 

55.1 

-94.9 

3.8 

42.5 

0.85 

MidLat  Day 

5.31 

88.7 

40.2 

82.1 

-94.7 

38.0 

66.1 

1.03 

HighLat  Day 

10.37 

54.7 

56.2 

151.6 

-25.9 

66.9 

137.5 

1.23 

MidLat  Day 

14.78 

45.9 

32.4 

40.0 

14.7 

44.5 

44.4 

11.86 

LowLat  Day 

16.88 

73.6 

-5.2 

246.6 

30.7 

5.7 

253.0 

0.11 

MidLat  Night 

19.03 

110.8 

-40.6 

274.0 

28.2 

-29.2 

280.8 

0.29 

HighLat  Night 

21.90 

141.6 

-57.1 

326.9 

-15.0 

-44.9 

322.4 

0.45 

Table  2.  Additional  Inputs  for  SAMM  and  SAG 


Parameter 

Value 

Ap 

18  (Low),  42  (Mean),  80  (High) 

Julian  day 

119 

F10.7 

l58  ^ 

F10.7A 

221 

Aerosol  model 

Maritime 

Scattering  phase  function 

Henyey-Greenstein,  g=0.9 

2.1.3  Outputs 

Each  SAMM  run  for  a  given  atmosphere  generates  a  file  containing  the  spectral  radiance  versus 
wavenumber  for  each  LOS  specified  in  the  input  list.  For  comparison  with  the  CIRRIS  data,  the 
SAMM  output  spectra  were  interpolated  between  tangent  heights,  integrated  over  the  detector 
vertical  footprint,  and  degraded  in  spectral  resolution,  as  described  in  Subsection  2.5.  The 
calculation  of  the  in-band  radiance  profiles  is  also  discussed  below. 
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2.1.4  In-Band  Radiance  Computations 


To  compare  with  the  CIRRIS  radiometer  data,  in-band  radiance  profiles  were  computed  by 
convolving  each  SAMM  spectrum  with  the  effective  radiometer  response  functions,  which  were 
measured  by  Utah  State  University  in  collaboration  with  AFRL  [Bingham,  1991].  The  approximate 
wavelengths  for  25%  of  peak  response  are  listed  in  Table  3. 


Table  3.  CIRRIS-IA  Radiometer  Filter  Bandpasses 


Filter 

Bandpass  (pm)* 

Major  Radiators 

0 

8.5-18 

CO^Cvj),  OjCVj) 

1 

4.1  -4.5 

COjfVj),  CO(Av=l ),  NO"(Av=l ) 

2 

4.9  -  7.0 

NO(Av=l),  H20(v2) 

3 

6.0 -8.8 

H20(V2),03(V,),  CH4(V4) 

FP2 

2.6 -3.3 

OH(Av=l),  COj,  aerosols 

5 

11.1  - 12.8 

CO2,  03(V3),  HNO3 

6 

17.9-23 

HjOfRot) 

7 

8.0-11.8 

03(^3)?  CO2 

®  25%  spectra 

response  limits 

The  detector  response  function  for  each  bandpass  is  the  product  of  a  relative  response  function, 
normalized  to  a  peak  value  of  1,  and  a  scale  factor.  The  detector  responsivity  values  used  in 
processing  the  raw  data  were  determined  from  absolute  calibrations  with  a  blackbody  source, 
assuming  that  the  relative  response  functions  were  those  of  the  filters  alone.  However,  actual 
relative  response  functions  were  subsequently  measured  and  found  to  be  significantly  different  due 
to  the  non-uniform  spectral  response  of  the  detectors.  The  scale  factor,  which  serves  to  compensate 
for  these  differences,  is  the  ratio  of  the  integral  of  the  blackbody-weighted  filter  response  function  to 
the  integral  of  the  equivalently  weighted  total  response  function.  We  performed  approximate 
calculations  of  the  scale  factors,  assuming  a  300  K  blackbody  calibration  source,  and  found  them  to 
agree  to  within  10-20%  with  results  obtained  from  comparisons  of  radiometer  and 
bandpass-integrated  interferometer  data  taken  simultaneously.  Taking  into  account  the  likely 
uncertainties  in  the  interferometer  calibration  and  in  the  calculations,  our  best  estimates  for  the  scale 
factors  are  given  in  Table  4. 


Table  4.  Scale  Factors  for  the  CIRRIS-IA  Radiometer  Response  Functions 


Filter 

Scale  Factor 

Uncertainty 

0 

1.1 

30% 

1 

1.0 

20% 

2 

1.3 

20% 

3 

1.2 

20% 

FP2 

1.2 

30% 

5 

1.1 

20% 

6 

1.5 

25% 

7 

1.0 

20% 

2.1.5  Spectral  Radiance  Computations 

The  comparisons  of  SAMM  spectra  with  CIRRIS  data  for  a  specified  field  of  view  (FOV)  take 
into  account  both  the  tangent  height  and  the  detector  vertical  footprint.  CIRRIS  spectra  for  tangent 
heights  that  lie  between  the  pre-calculated  values  are  logarithmically  interpolated.  The  integration 
over  the  detector  footprint,  which  on  average  is  16  km  for  the  largest  detector  (number  3-2)  in  focal 
plane  (FP)  3,  is  performed  using  a  three-point  average  of  interpolated  spectra.  Tests  of  this  approach 
using  an  example  of  extreme  altitude  variation  (a  2.5  km  scale  height,  which  occurs  for  O3  and  OH 
emissions  in  the  lower  thermosphere  above  90  km)  indicate  a  maximum  error  of  30%,  which  is 
comparable  to  the  tangent  height  measurement  uncertainty  of  about  ±1  km.  Both  the  effect  of  the 
footprint  size  and  the  error  in  the  computation  decrease  very  rapidly  as  the  scale  height  increases, 
and  are  typically  less  than  10%. 

The  CIRRIS  interferometer  spectra  are  1.0-cm  '  resolution  FWHM,  output  at  approximately 
0.5-cm  '  increments,  whereas  the  SAMM  calculated  spectra  are  actually  radiance  sums  for 
non-overlapping  l-cm  '-wide  bins.  Thus,  a  quantitative  comparison  of  the  two  data  sets  is  possible 
only  at  a  reduced  spectral  resolution.  The  CIRRIS  and  SAMM  spectra  were  both  degraded  using  the 
same  resolution  function.  The  resolution  FWHM  was  chosen  based  on  the  desired  final  resolution, 
assuming  that  the  resolution  FWHM  and  the  intrinsic  1-cm'  resolution  of  the  undegraded  spectra 
add  in  quadrature,  as  would  be  the  case  for  Gaussian  resolution  functions.  For  example,  nominal 
3-cm  '  resolution  spectra  were  calculated  using  a  resolution  FWHM  of  V(3--l^)  =  2.83  cm  '. 
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Several  different  resolution  functions  were  tried,  including  square,  triangular,  Gaussian,  and  sech 
(=l/cosh)  functions.  Of  these,  the  last  two  gave  the  smoothest  line  shapes.  The  spectral 
comparisons  plots  in  the  data  base  were  calculated  using  the  sech  function. 

Typically,  a  modest  spectral  degradation  resulted  in  a  very  good  match  between  the  observed  and 
calculated  spectral  shapes.  An  exception  to  this  is  when  the  spectrum  consists  of  a  series  of  evenly 
spaced  lines,  such  as  are  found  in  the  CO  and  NO  bands.  The  depths  of  the  valleys  between  the  lines 
are  sensitive  to  the  valley  widths,  which  in  turn  depend  on  the  particular  bins  into  which  the  lines 
fall.  For  certain  line  spacings,  SAMM's  relatively  coarse  binning  can  result  in  artificial  oscillatory 
structure  in  the  degraded  spectrum,  even  at  a  resolution  sufficient  to  wipe  out  the  genuine  line 
structure.  The  sech  resolution  function,  which  has  less  steep  wings  than  a  Gaussian,  provides  some 
reduction  in  this  artificial  structure  and  was  therefore  used  for  the  final  calculations.  For  the 
calculations  in  the  database,  the  final  FWHM  ranged  from  2.5  to  4  cm''  and  the  degraded  spectra 
were  generated  at  1  -cm  '  intervals. 

2.2  CIRRISDATA 

2.2.1  Radiometer  Profiles 

Files  of  a  representative  subset  of  CIRRIS  radiometer  profiles  were  generated  using  the  OPUS 
software  package  developed  by  Boston  College  [Hegblom,  1991].  The  data  were  selected  on  the 
basis  of  both  quality  and  coverage  of  the  full  range  of  atmospheric  conditions.  Generally,  all  six  FP 
1  detectors  and  two  of  the  three  FP  2  detectors  (2-1  eind  2-3)  were  sampled.  Only  four  of  the  FP  1 
detectors  (1-1, 1-2,  1-3,  and  1-5)  were  sampled  with  filter  6,  since  the  other  two  detectors,  1-4  and 
1-6,  showed  much  higher  background  contamination  levels  with  that  filter. 

Many  of  the  profiles  used  in  the  database  comparisons  were  obtained  from  the  “staircase” 
scanning  modes  which  feature  rapid  vertical  jumps  in  between  constant-altitude  stares.  Therefore, 
frequent  sampling  of  the  radiometers  was  required  for  the  in-between  altitudes.  This  resulted  in  very 
large  data  files,  typically  2  Mb  each.  To  keep  the  total  database  manageable  while  maintaining  the 
overall  integrity  of  the  data,  an  altitude-dependent  algorithm  was  implemented  for  thinning  the 
number  of  data  points  in  each  scan.  The  data  were  placed  in  1  -km  bins  and  every  VN-th  point  was 
retained,  where  N  is  the  number  of  points  in  the  bin.  In  addition,  data  falling  below  a  minimum 
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radiance  threshold  (in  the  10  ’  W/cmVsr  range)  were  eliminated,  since  these  data  are  potentially 
affected  by  near-field  shuttle  contamination  or  off-axis  telescope  leakage. 

2.2.2  Spectra 

Using  OPUS,  files  of  spectral  radiance  versus  wavenumber  were  generated  from  measurements 
taken  with  an  individual  FP  3  detector  during  a  full-length  interferometer  scan  at  a  constant  tangent 
height  or  angle  to  the  nadir.  The  FWHM  resolution  of  the  full-length  scans  is  1.0  cm  '  with 
Kaiser-Bessel  apodi2ation  (a=2).  The  corresponding  spectral  noise  files  were  also  compiled,  and 
were  used  to  flag  spectral  regions  of  poor  signal-to-noise  for  exclusion  from  the  plotting  routine. 
The  noise  level  is  calculated  by  OPUS  from  signal  fluctuations  in  the  uncalibrated  spectrum  at 
wavelengths  where  the  detectivity  is  negligible. 

Representative  CIRRIS-1 A  data  were  selected  primarily  from  the  staircase”  altitude-scanning 
modes  of  data  blocks  PC  1 2  A  and  PC12B,  which  sampled  the  atmospheric  radiance  in  approximately 
12-km  tangent  height  steps  during  mid-latitude  nighttime  and  high-latitude  daytime  conditions. 
Different  staircase  scans  used  different  bandpass  filters  in  order  to  isolate  a  particular  spectral  region 
from  much  stronger  interfering  emissions,  such  as  from  CO^  (v^)  and  O3  (Vj),  that  can  generate 
photon  noise  and  harmonic  distortion  components.  The  approximate  filter  bandpasses  are  listed  in 
Table  5.  For  a  given  staircase  scan,  a  spectrum  was  selected  every  1  or  2  altitude  steps  or, 
equivalently,  about  every  12  to  24  km.  Typically,  the  data  from  the  largest  and  most  sensitive 
detector,  3-2,  were  chosen.  When  detector  3-2  was  saturated  or  when  a  smaller  vertical  footprint  was 
desired,  data  taken  by  one  of  the  smaller,  less-sensitive  detectors  (3-3  to  3-5)  were  used.  The  vertical 
footprint,  which  is  tangent  height-dependent,  is  about  13  to  1 8  km  for  detector  3-2  and  about  2  to  3 
km  for  detectors  3-3,  3-4,  and  3-5.  The  downward-looking  measurements  were  performed  in  a 
horizontal  stare  mode  during  block  PC5 1 ;  a  few  representative  spectra  from  detector  3-2  were  chosen 
to  include  nadir  views  in  the  database. 
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Table  5.  Nominal  CIRRIS-1 A  Interferometer  Bandpasses 


Filter  No. 

Bandpass  (pm) 

Major  Radiators 

0 

All  (no  filter) 

CO,iv,),  03(V3) 

1 

<4.9 

C02(v3),C0(Av=1),03 

2 

<3.7 

OH(Av=l),  COj,  aerosols 

3 

4.8-13 

H20(v3),03(v,),N0(Av=1) 

5 

11  - 13 

O3,  OH(Rot),  CO2 

6 

>  17 

H20(Rot) 

7 

8-13 

O3,  OH(Rot),  CO2 

2.3  RESULTS  AND  DISCUSSION 
2.3.1  Day  and  Night  Comparisons 

The  CIRRIS-1  A-SAMM  database  of  in-band  radiance  profiles  and  spectra  provides  the  first 
comprehensive  and  systematic  comparisons  of  real-world  upper  atmospheric  IR  background  data 
with  code  calculations  for  a  wide  range  of  molecular  bands  and  geographic  and  atmospheric 
conditions.  This  is  significant  because  the  database  provides  the  basis  for  assessing  and  enhancing 
the  predictive  capabilities  of  atmospheric  radiance  and  transmission  codes  required  by  systems 
engineers  to  optimize  system  design  and  operability.  The  major  results  of  the  in-band  radiance 
profiles  and  spectra  comparisons,  particularly  absolute  radiance  comparisons,  are  summarized  below 
according  to  wavelength  region. 

Radiance  profile  and  altitude  dependent  spectral  plots  comparing  CIRRIS  data  with  SAMM  were 
generated  using  the  SAG  model  atmospheres  "mOl "  (1 .92  hrs  local  time,  mean  Ap,  27°  S  latitude) 
and  "mlO"  (10.37  hrs  local  time,  mean  Ap,  67°N  latitude),  which  are  described  in  detail  in  Tables  1 
and  2.  The  major  uncertainties  in  the  radiometer  data  are  20-30%  in  the  absolute  radiances,  reflecting 
the  instrument  calibrations,  and  ~2  km  in  the  tangent  height  determinations.  To  avoid  plotting  very 
noisy  spectral  data,  the  spectral  noise  file  from  OPUS  was  compared  with  the  radiance  file,  and  when 
the  RMS  signal-to-noise  ratio  was  less  than  a  pre-defined  value  (in  this  case,  4),  the  data  were 
replaced  with  a  non-plottable  (negative)  number. 
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Longer  than  17  p,m  (below  600  cm  ')  The  main  source  of  emission  at  wavelengths  longer  than  1 7 
\im  is  pure  rotation  transitions  of  H^O.  Representative  night  and  day  radiance  profiles  are  shown  in 
Figures  2  and  3,  respectively.  The  measured  and  calculated  in-band  radiance  profiles  (filter  6)  agree 
to  within  a  factor  of  two  or  better  up  to  80  km.  Below  about  60  km  the  emission  is  very  optically 
thick  and  thus  fairly  insensitive  to  the  H^O  concentration.  Above  80  km,  a  quantitative  comparison 
is  more  difficult  due  to  the  rapid  falloff  in  atmospheric  radiance  and  the  presence  of  near-field 
contamination  from  H2O  outgassing  from  the  shuttle  [Zhou  et  al,  1993]. 

13-17  pm  (600-800  cm  ')  Radiation  in  the  13-17  pm  spectral  region  dominates  the  1-20  pm  IR 
region  at  most  altitudes  and  originates  from  several  CO^  v^-type  bands.  Over  most  of  the  altitude 
range  from  0  to  240  km,  the  agreement  between  the  radiance  measurements  and  calculations  is  good 
(see  Figures  4-6),  and  the  systematics  of  the  observed  variations  in  the  radiometer  data  (filter  0)  near 
140  km  are  reasonably  well  reproduced  by  the  calculations  [Wise  et  al,  1993].  However,  near  the 
mesopause  (80-100  km)  there  is  typically  more  fine  structure  in  the  radiance  profile  than  is 
calculated.  While  the  structure  is  quite  variable,  often  there  is  a  narrow  layer  of  enhanced  radiance 
near  85  km,  especially  between  midnight  and  dawn  in  the  southern  hemisphere,  as  shown  in  Figures 
4-6.  This  shows  up  in  spectra  from  the  smaller  focal  plane  3  detectors  as  well  as  in  the  radiometer 
profiles,  but  not  in  the  spectra  from  the  large  detector  3-2.  Since  the  CO^  bands  are  at  or  near 
thermal  equilibrium  and  the  mixing  ratio  is  well  known  at  85  km,  vertical  structure  in  the 
temperature  profile  which  is  not  modeled  by  SAMM/SAG  is  the  best  explanation  for  the  observed 
radiance  layering.  This  would  be  consistent  with  recent  LIDAR  data  [Chanin,  1 992]  that  often  show 

a  "double  mesopause"-shaped  temperature  profile  that  is  considerably  different  from  the  MSIS 
model  profiles  used  by  SAG. 

The  spectral  features  of  C02(Vj)  under  daytime  conditions  are  illustrated  in  Figures  7  and  8  for 
tangent  altitudes  of  1 3  km  and  59  km,  respectively,  while  Figures  9  and  1 0  show  C02(V2)  at  night  for 
tangent  altitudes  of  1 3  and  67  km,  respectively.  Again,  the  agreement  between  the  CIRRIS-1 A  data 
and  SAMM  is  quite  remarkable. 
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Figure  2.  Comparisons  of  the  Nighttime  CIRRIS-IA  Radiance  Profile  in  the  18-23  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of  -43° 


Figure  3.  Comparisons  of  the  Daytime  CIRRIS-IA  Radiance  Profile  in  the  18-23  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of  65° 


13 


Figure  4.  Comparisons  of  the  Nighttime  CIRRIS-IA  Radiance  Profile  in  the  8-18  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of -29° 


Figures.  Comparisons  of  the  Terminator  CIRRIS-IA  Radiance  Profile  in  the  8-18  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of  10° 
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Figure  6.  Comparisons  of  the  Daytime  CIRRIS-1 A  Radiance  Profile  in  the  8-18  |im  Bandpass 
with  SAMM  Calculations  at  a  Latitude  of  28° 


400  500  600  700  800  900  1000  1100  1200 

Wavenumber  (cm-1) 

Figure  7.  Daytime  Comparisons  of  CIRRIS-IA  Data  in  the  8-25  ^m  Spectral  Region  with 
SAMM  Calculations  at  a  Tangent  Height  of  13  km  and  a  Latitude  of  63° 
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Figure  8.  Daytime  Comparisons  of  CIRRIS-IA  Data  in  the  8-25  pm  Spectral  Region  with 
SAMM  Calculations  at  a  Tangent  Height  of  59  km  and  a  Latitude  of  65° 


Figure  9.  Nighttime  Comparisons  of  CIRRIS-1 A  Data  in  the  8-25  pm  Spectral  Region  with 
SAMM  Calculations  at  a  Tangent  Height  of  13  km  and  a  Latitude  of -42° 
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Figure  10.  Nighttime  Comparisons  of  CIRRIS-IA  Data  in  the  8-25  pm  Spectral  Region  with 
SAMM  Calculations  at  a  Tangent  Height  of  67  km  and  a  Latitude  of  -43° 

11-13  pm  (800-900  cm'')  In  the  11  -  13  pm  atmospheric  window  region,  the  major  upper 
atmospheric  emissions  are  from  O3  hot  bands  and  solar-pumped  CO2  laser  bands.  Chemiluminescent 
emissions  from  OH  rotation  lines  [Smith  et  al,  1992;  Dodd  et  al,  1994]  which  are  not  currently 
modeled  by  SAMM,  have  also  been  observed  by  the  CIRRIS-IA  experiment  in  the  mesopause 
region,  primarily  at  night.  The  O3  hot  band  intensities  tend  to  be  overpredicted  by  SAMM, 
especially  near  800  cm"';  this  is  primarily  due  to  uncertainties  in  the  spectroscopy  and  kinetics  of 
high  vibrational  states  of  O3.  Nevertheless,  the  in-band  radiance  profiles  calculated  using  SAMM  for 
filter  5  are  typically  within  a  factor  of  two  of  the  CIRRIS  radiometer  data  at  most  altitudes,  as  shown 
in  Figures  1 1  and  12.  In  the  lower  mesosphere,  the  largest  differences  occur  during  the  daytime  (see 
Figure  12) 

Below  about  30  km  tangent  height,  the  CIRRIS  spectra  show  emission  bands  near  880  cm  '  due 
to  HNO3  and  near  850  cm  '  due  to  CFC-1 1  [Bingham  et  al,  1993],  see  Figures  7-10.  The  HNO3 
band  is  prominent  in  the  calculated  spectra  up  to  about  40  km  but  only  up  to  about  30  km  in  the 
CIRRIS  data.  CFC's  are  not  modeled  currently  by  the  SAMM  code  but  have  recently  been  added  to 
MODTRAN. 
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Figure  11.  Comparisons  of  the  Nighttime  CIRRIS-IA  Radiance  Profile  in  the  11-13  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of -39° 


Figure  12.  Comparisons  of  the  Daytime  CIRRlS-lA  Radiance  Profile  in  the  11-13  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of  68° 
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8-1 1  um  (900-1200  cm~')  The  major  emissions  in  the  8- 1 1  |am  wavelength  region  are  from  the  O3  v, 
and  V3  fundamental  modes  and  solar-pumped  COj  bands.  Representative  radiance  profiles  at  day  and 
night  are  shown  in  Figures  13  and  14.  Spectral  features  are  again  illustrated  in  Figures  7-10.  As 
seen  in  the  filter  7  radiometer  data,  near  80  km  the  SAMM  code  consistently  overpredicts  the  ozone 
emission  by  about  a  factor  of  two,  which  is  likely  due  to  an  overestimation  of  the  O3  concentration 
by  SAG.  At  other  altitudes,  the  average  agreement  is  good,  and  the  diurnal  variability  is  realistically 
modeled  by  the  calculations.  Below  30  km  tangent  height,  an  emission  band  due  to  CFC-12 
[Bingham  et  al,  1993]  can  also  be  seen  in  the  CIRRIS  spectra  near  920  cm  '. 


Radiance  (W/cm2/sr) 

Figure  13.  Comparisons  of  the  Nighttime  CIRRIS-IA  Radiance  Profile  in  the  8-12  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of  -29° 
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Figure  14.  Comparisons  of  the  Daytime  CIRRIS-1 A  Radiance  Profile  in  the  8-12  pm  Bandpass 
with  SAMM  Calculations  at  a  Latitude  of  50° 


5-8  pm  (1200-2000  cm  ')  The  major  emitters  between  5  and  8  pm  are  H20(V2),  CH4(v4)  near  1311 
cm  ,  and  NO(Av-l ).  For  the  H2O  band,  the  measurement  and  model  comparison  is  consistent  with 
that  of  the  rotation  line  emission.  The  spectral  data  comparisons  also  show  that  the  methane 
emission  is  overpredicted  in  the  upper  mesosphere  but  predicted  well  at  lower  tangent  heights.  The 
source  of  the  upper  mesosphere  discrepancy  is  under  investigation  and  may  be  due  to  an  error  in  a 
SAMM  spectral  parameter  file  that  has  since  been  corrected.  In  the  mesopause  region,  the  CIRRJS 
radiometer  data  for  filters  2  and  3  typically  show  emission  layers  similar  to  those  in  the  filter  0  data, 
which  cover  the  15  pm-region  CO2  bands.  The  presence  of  this  feature  is  consistent  with  the 
explanation  of  temperature  structure.  Below  the  mesopause,  the  radiance  profiles  for  filters  2  and  3 
are  dominated  by  H20(V2)  (see  Figures  1 5-1 8). 


20 


Figure  15.  Comparisons  of  the  Nighttime  CIRRIS-1 A  Radiance  Profile  in  the  5-7  pm  Bandpass 
with  SAMM  Calculations  at  a  Latitude  of  -28° 


Figure  16.  Comparisons  of  the  Daytime  CIRRIS-1  A  Radiance  Profile  in  the  5-7  pm  Bandpass 
with  SAMM  Calculations  at  a  Latitude  of  38° 
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Figure  17.  Comparisons  of  the  Nighttime  CIRRIS-1 A  Radiance  Profile  in  the  6-9  pm  Bandpass 
with  SAMM  Calculations  at  a  Latitude  of -42° 


Figure  18.  Comparisons  of  the  Daytime  CIRRIS-1  A  Radiance  Profile  in  the  6-9  pm  Bandpass 
with  SAMM  Calculations  at  a  Latitude  of  64° 


The  5-8  nm  emission  above  the  mesopause  is  dominated  by  NO(Av=  1 ),  see  Figures  15,16,  and 
19.  Both  the  measurements  and  calculations  of  the  v=l  band  emission  vary  by  over  an  order  of 
magnitude  due  mainly  to  the  dependence  of  NO  concentration  on  solar  angle,  latitude,  and 
geomagnetic  activity.  Due  to  these  atmospheric  variations,  the  measurement  and  model  agreement  is 
erratic.  The  largest  differences  are  found  for  the  nighttime  data,  where  there  have  been  few 
measurements  for  the  development  of  an  empirical  NO  concentration  model  [Smith  et  al. ,  1 993]  for 
use  in  SAG  and  other  model  atmosphere  climatologies. 


Wavenumber  (cm-1) 


Figure  19.  Nighttime  Comparisons  of  CIRRIS-IA  Data  in  the  5-6  pm  Spectral  Region  with 
SAMM  Calculations  at  a  Tangent  Height  of  111  km  and  a  Latitude  of  -44° 


The  daytime  CIRRIS  spectra  show  a  contribution  from  high-v,J  states  of  NO  in  the  1200-1 700 
cm  '  (6-8  pm)  region  [Armstrong  et  al,  1993;  Smith  and  Ahmadjian,  1993].  The  source  of  this 
emission  is  believed  to  be  primarily  the  reaction  of  fast  N(‘'S)  atoms  with  Oj  [Duff  et  al,  1994; 
Sharma  et  al.,  1998]  and  a  model  upgrade  incorporating  this  mechanism  is  planned  for  AFRL 
radiance  models. 
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jrj  uni  (2000-2500  cm ')  The  radiation  in  the  4-5  jim  wavelength  region  is  dominated  by  the  COj 
Vj-type  bands,  with  smaller  contributions  from  CO(Av=l)  emissions  and  NO*‘Av=l) 
chemiluminescence  at  high  altitudes  during  the  day  time  [Winick  etal.,  1992],  The  in-band  radiance 
profile  (filter  1)  and  the  COj  spectral  radiance  are  modeled  well  by  the  SAMM  Code  for  both 
nighttime  (Figure  20)  and  daytime  (Figure  2 1 )  conditions,  typically  to  within  better  than  a  factor  of 
two.  Nighttime  and  Daytime  spectral  comparisons  are  shown  in  Figures  22  and  23  at  62  km  and  59 
km,  respectively.  The  emission  from  NO^  is  not  accounted  for  in  these  non-auroral  SAMM 
calculations.  Detailed  analysis  [Winick  et  al.,  1992]  of  the  CIRRIS-IA  data  has  made  the  first 
identification  of  non-auroral  NO*  emissions,  which  are  an  additional  new  source  of  radiation  in  the 
4-5  ^im  spectral  region  for  certain  atmospheric  conditions. 


Figure  20.  Comparisons  of  the  Nighttime  CIRRIS-IA  Radiance  Profile  in  the  4.1-4.5  |im 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of -13° 
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Figure  21.  Comparisons  of  the  Daytime  CIRRIS-IA  Radiance  Profile  in  the  4.1-4.5  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of  69° 
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Figure  22.  Daytime  Comparisons  of  CIRRIS-IA  Data  in  the  4-5  pm  Spectral  Region  with 
SAMM  Calculations  at  a  Tangent  Height  of  62  km  and  a  Latitude  of  70° 
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Figure  23,  Nighttime  Comparisons  of  CIRRIS-1 A  Data  in  the  8-25  pm  Spectral  Region  with 
SAMM  Calculations  at  a  Tangent  Height  of  59  km  and  a  Latitude  of -40° 


The  downward-looking  spectra,  taken  with  interferometer  filter  I,  are  dominated  by  CO^CVj) 
emissions  as  well  as  blackbody  emission  and  solar  scattering  from  the  earth's  surface  and  from 
clouds.  A  considerable  portion  of  the  COj  emission  is  from  above  30  km  due  to  optical  opacity 
effects  and  is  characterized  by  a  non-thermal  vibrational  state  distribution.  The  total  radiance  level 
varies  considerably  from  one  CIRRIS  scan  to  another.  The  major  differences  between  the  CIRRIS 
data  and  the  SAMM  calculations  are  likely  due  to  uncertainties  in  the  modeling  of  the  earth  surface 
temperature  and  the  lack  of  cloud  cover  in  the  calculation. 

Bejow  4  pm  (above  2500  cm  ')  The  radiometer  data  below  4  pm,  from  FP  2,  contain  contributions 
from  solar-pumped  COj  bands,  the  Av=  1  Meinel  bands  of  OH,  and,  at  the  lowest  altitudes,  solar 
scattering  from  aerosols  and/or  clouds.  The  OH  Meinel  bands  arise  from  the  H+0,  reaction  and  are 
a  more  prominent  at  night  when  the  O3  concentration  is  larger.  At  a  14  to  15  km  tangent  height, 
there  is  evidence  of  scattering  from  clouds,  which  are  not  modeled  in  these  calculations,  in  the 
daytime  CIRRIS  spectra  between  2600  and  2900  cm"'.  At  higher  wavenumbers,  the  spectra  have 
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poor  signal-to-noise,  and  only  the  strongest  OH  lines  are  identifiable.  The  limb  radiance  profiles 
(see  Figures  24  and  25),  which  have  good  signal-to-noise,  are  generally  modeled  reasonably  well  by 
SAMM,  although  there  is  more  variability  in  the  data  than  predicted.  Spectral  data  of  higher  S/N 
would  be  required  to  the  separate  the  radiance  contributions  from  each  individual  IR  emitter. 

23.2  Solar  Terminator  Comparisons 

IR  atmospheric  background  radiance  varies  substantially  over  the  diurnal  cycle.  In  many  spectral 
regions,  abrupt  changes  that  can  occur  at  the  solar  terminator  (earth  shadow  edge)  can  seriously 
impact  the  performance  of  many  sensing  systems.  These  terminator  radiance  variations  arise  from 
several  different  phenomena.  One  is  a  photochemically-induced  change  in  molecular  number 
density  caused  by  photochemical  reactions.  This  occurs  with  both  ozone  (O3),  which  is  photolyzed 
at  dawn  and  regenerated  at  dusk,  and  OH,  which  is  formed  from  the  reaction  of  O3  with  H  atoms.  O3 
emits  strongly  in  the  8-12  pm  region,  and  OH  chemiluminescence  is  prominent  around  the 
mesopause  at  wavelengths  below  3.7  pm.  Another  important  source  of  terminator  radiance  variation 
is  a  ehange  in  the  vibrational  excitation  of  the  molecule  with  solar  illumination.  This  non-LTE  effect 
occurs  in  a  number  of  species  including  COj,  and  dramatically  affects  the  4.3  pm  and  2.7  pm  bands 
throughout  the  mesosphere  and  thermosphere.  Data-model  comparisons  indicate  that  solar  pumping 
of  CO2  is  noticeable  in  the  4.3  pm  region  even  in  nadir  views,  which  typically  sense  only  the  lower 
atmosphere. 

Both  the  SHARC-3  and  SAMM  codes  permit  aecurate  modeling  of  limb  radiance  at  the  solar 
terminator.  They  allow  the  LOS  to  pass  through  multiple  atmospheric  regions,  corresponding  to 
different  solar  zenith  angles  (SZAs),  in  which  complete,  first-principles  ealculations  of  the  non-LTE 
kinetics  and  radiation  transport  (including  layer-layer  radiative  excitation  effects)  are  performed.  In 
a  preliminary  demonstration  of  terminator  modeling  capability  [Adler-Golden  and  Matthew,  1993], 
SHARC-3  was  used  to  model  CO2  4.3  pm  data  from  the  SPIRE  rocket  experiment  as  well  as  limited 
CIRRIS-IA  dusk  terminator  data  on  O3  in  the  8-12  pm  region  and  OH  in  the  fundamental  band 
region  (below  3.7  pm).  Good  overall  agreement  was  found. 
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Figure  24.  Comparisons  of  the  Nighttime  CIRRIS-IA  Radiance  Profile  in  the  2.6-3.3  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of -36° 


Figure  25.  Comparisons  of  the  Daytime  CIRRIS-IA  Radiance  Profile  in  the  2.6-3.3  pm 
Bandpass  with  SAMM  Calculations  at  a  Latitude  of  69° 
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The  current  investigation  presents  a  more  comprehensive  set  of  terminator  data-model  radiance 
comparisons  using  results  from  the  CIRRIS-IA  experiment.  Comparisons  of  representative 
CIRRIS-1 A  data  with  SAMM  code  predictions  are  presented  for  the  altitude  and  wavelength  ranges 
with  the  largest  terminator  variations  and  the  greatest  importance  for  sensing  systems.  The 
comparisons  include  in-band  and/or  spectral  radiances  for  the  8-12  |am,  A<3.7  jam,  and  4.3  |im 
regions  in  the  ~50-100  km  altitude  range. 

The  basic  comparisons  were  performed  using  model  atmospheres  from  the  SHARC/SAMM 
atmosphere  generator  (SAG)  as  inputs  to  the  SAMM  code.  In  addition,  for  emissions  from  the 
odd-oxygen  species  O3  and  OH  further  SAMM  calculations  were  performed  using  theoretical 
profiles  from  a  one-dimensional  photochemical  model  (1-D  PC)  [Makhlouf  et  al.,  1994].  These 
profiles  were  found  to  yield  better  agreement  with  the  CIRRIS  data. 

CIRRIS  terminator  data  also  exist  for  other  spectral  regions,  including  5-7  p.m  and  18-23  pm. 
These  data  have  been  examined  for  solar  excitation-related  variations  and  are  discussed  briefly  in 
Section  2.4.2.3. 

2.3.2. 1  8-12  um  and  A<3.7  pm  Regions 

The  8-12  pm  region  and  the  region  below  3.7  pm  contain  the  major  IR  emissions  of  the 
odd-oxygen  species,  O3  and  OH.  Since  OH(v)  is  formed  directly  from  O3  (via  reaction  with  H 
atoms),  the  O3  and  OH  emissions  track  each  other  across  the  terminator.  Dawn  and  dusk  must  be 
distinguished  since  the  concentrations  and  hence  emissions  are  not  equivalent  functions  of  solar 
angle  in  the  two  cases.  Solar-excited  bands  from  COj  also  appear  in  these  regions,  but  are  generally 
much  smaller  contributors  to  the  in-band  radiance  in  the  CIRRIS  data. 

2.3.2. LI  In-Band  Radiances  (FP3  Filter  7  and  FP2) 

2.3.2.1.1.1  Methodology 

CIRRIS-1  A  radiometer  data  from  focal  plane  3,  filter  7  (8-12  pm)  and  focal  plane  2  (A.<3.7  pm) 
were  obtained  from  the  near-horizontal  dawn  terminator  scan  in  Block  PC12B,  Mode  24-2.  The 
filter  7  emissions  are  from  the  O3  v,  and  V3  bands  and  solar-excited  CO2  laser  bands.  The  FP2 
radiance  derives  from  the  OH(Av=l )  and  solar-excited  2.7  pm  COj  bands.  The  views  are  due  west 
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across  the  terminator  and  the  tangent  heights  range  from  approximately  74  km  (detector  3-6)  to  96 
km  (detector  3-1).  Data  were  tabulated  for  specific  times  corresponding  to  approximately  2° 
increments  in  tangent  point  solar  zenith  angle  (SZA)  over  a  range  of  100°  to  88°.  During  this 

sequence  of  measurements  the  spacecraft  traveled  from  approximately  28°  N  latitude  to  41°  N 
latitude. 

The  SAMM  calculations  were  performed  using  a  set  of  seven  model  atmospheres  generated  for  a 
SZA  range  of  101-87°,  a  latitude  range  of 27-42°  N,  and  a  local  time  range  of 4.58-5.34  hrs.  These 
correspond  to  values  for  the  tangent  point  viewed  by  the  centers  of  the  CIRRIS  focal  planes  during 
the  terminator  traversal.  In  calculating  the  spectra  and  radiances  for  an  individual  detector  on  the 
focal  plane,  the  actual  tangent  height  and  tangent  point  SZA  for  that  detector  were  used;  the  latitudes 
and  local  times  of  the  model  atmospheres  along  the  LOS  differ  slightly  from  the  actual  values.  To 
exactly  match  the  actual  values,  a  separate  set  of  atmospheres  and  vibrational  state  calculations 
would  have  been  required  for  each  location  of  the  spacecraft.  A  direct  comparison  between  the 
current  approach  and  a  calculation  using  the  "exact"  atmospheric  conditions  yielded  only  1%-level 
LOS  radiance  differences. 

2.3.2.  LI. 2  Results 

The  data-model  comparisons  were  made  at  tangent  heights  of  74  to  96  km.  However, 
comparisons  are  shown  in  Figures  26  and  27  (FP3  filter  7)  for  the  highest  tangent  heights  (96  and  92 
km)  and  Figures  28  and  29  (FP2)  at  tangent  heights  of  83  and  74  km.  The  initial  SAMM 
calculations  were  based  on  model  atmospheres  from  the  SHARC/SAMM  Atmosphere  Generator 
(SAG)  and  are  labeled  SAG  .  The  radiance  changes  with  SZA  are  superimposed  on  typically 
smaller  fluctuations  associated  with  ±2  km  changes  in  the  tangent  point  during  the  measurement 
period;  the  tangent  altitude  listed  is  the  average  for  that  detector. 

In  both  bandpasses  the  largest  changes  in  radiance  with  SZA  occur  at  the  highest  altitudes;  there 
the  absolute  agreement  with  the  CIRRIS  data  is  quite  good.  The  location  of  the  terminator  edge  and 
the  radiance  trend  remain  qualitatively  reasonable  at  the  lower  altitudes;  however,  the  absolute 
radiances  are  overpredicted  by  up  to  a  factor  of  4  (in  filter  7  at  74  km). 
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Figure  26.  Dawn  Terminator  8>12  |im  Radiometric  Comparisons  for  Average  ±2  km)  at  a 
Tangent  Height  of  96  km.  CIRRIS-IA  Data  are  From  Block  PC12B  Mode  24-2. 


Figure  27.  Dawn  Terminator  8-12  pm  Radiometric  Comparisons  for  Average  ±2  km)  at  a 
Tangent  Height  of  92  km.  CIRRIS-1 A  Data  are  From  Block  PC12B  Mode  24-2. 
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Figure  28.  Dawn  Terminator  Radiometric  Comparisons  for  Focal  Plane  2  at  at  a  Tangent 
Height  of  83  km.  Data  were  Taken  Simultaneously  with  those  in  Figures  27  and  28. 
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As  indicated  earlier,  the  emissions  at  these  wavelengths  are  dominated  by  OH  and  O3,  especially 
on  the  night  side,  and  track  each  other  across  the  terminator  since  both  the  OH(v)  and  03(v) 
concentrations  are  proportional  to  the  O3  number  density.  The  radiance  overpredictions  by  SAMM 
in  the  74  to  83  km  tangent  height  range  are  associated  with  excess  O3  in  the  SAG  model 
atmospheres,  which  in  turn  is  ascribed  mainly  to  an  overestimation  of  the  atomic  oxygen 
concentration  at  these  altitudes,  particularly  at  dawn.  The  atomic  oxygen  profile  in  SAG  is 
approximate;  as  discussed  in  [Adler-Golden  1993],  the  profile  is  scaled  from  the  MSISE-90 
empirical  model  [Hedin,  1991]  using  a  simple  day/night  ratio,  an  approach  which  does  not  accurately 
treat  the  time-dependence. 

More  quantitative  modeling  of  OH  and  O3  emissions  in  the  70-85  km  altitude  range  requires  an 
accurate  simulation  of  the  time-dependent  concentration  profiles  for  both  the  odd-oxygen  and 
odd-hydrogen  species.  To  demonstrate  how  improved  species  profiles  can  be  used  with  SAMM  for 
improved  radiance  predictions,  a  second  set  of  calculations  was  performed  using  O3,  O,  and  H 
concentrations  from  a  one-dimensional  time-dependent  photochemical  model  that  has  been 
developed  at  AFRL  [Makhlouf  et  al. ,  1 994].  The  results  are  labeled  "1-D  PC"  in  Figures  26-29.  The 
CIRRIS-SAMM  agreement  is  seen  to  be  very  good  in  both  the  FP3  Filter  7  and  FP2  bandpasses  over 
the  full  altitude  range. 

It  is  concluded  that  in  order  to  provide  the  most  accurate  terminator  simulations  for  systems 
applications,  results  from  this  or  similar  time-dependent  photochemical  calculations  should  be 
incorporated  in  the  model  atmosphere.  A  simple  but  somewhat  less  accurate  alternative  would  be  to 
systematically  scale  back  the  O  atom  profile  in  SAG  near  80  km,  to  yield  reduced  O3  and  OH. 

2.3.2. 1.2  Spectra 

2.3.2. 1.2.1  Methodology 

Terminator  spectra  from  the  CIRRIS- 1 A  interferometer  in  the  wavelength  region  of  Focal  Plane 
2  are  rather  limited,  and  the  declining  response  of  the  interferometer  detectors  makes  their 
signal-to-noise  insufficient  for  separating  the  daytime  COj  2.7  pm  and  OH(Av=l)  contributions 
without  signal  averaging.  A  detailed  analysis  of  the  nighttime  OH  data  has  recently  been  published 
[Dodd  et  al. ,  1 994].  No  spectral  analysis  of  this  wavelength  region  has  therefore  been  attempted  for 
this  report. 


33 


Spectra  of  the  dawn  terminator  in  the  8-1 2  pm  region  were  taken  in  CIRRIS  during  Blocks  1 1 G, 
16,  and  42B  with  medium  (3  cm  '  FWHM)  resolution.  Similar  data  were  taken  in  the  dusk 
terminator  during  Block  1 2B.  For  the  current  comparisons,  we  selected  a  subset  of  1 4  dawn  spectra 
at  tangent  heights  near  70  km,  78  km,  and  86  km  covering  SZA's  from  87°to  101°. 

The  SAMM  simulations  were  performed  using  model  atmospheres  calculated  for  the  specific 
times  and  locations  along  the  LOS  using  the  1  -D  photochemical  model  for  odd-0  species  and  SAG 
for  the  remaining  species  and  temperature  profiles.  The  SAMM  spectra  were  degraded  to  match  the 
CIRRIS  resolution.  Unlike  the  spectral  calculations  illustrated  in  Section  2.4. 1 ,  these  calculations  do 
not  account  for  the  detector  footprint  and  therefore  may  underpredict  the  radiance  by  up  to  one-third 

at  some  wavelengths.  This  is  not  a  major  consideration  in  evaluating  the  overall  data-model 
agreement,  however. 

2.3.2.L2.2  Results 

The  data-model  comparisons  are  shown  in  detail  in  [Adler-Golden  et  al,  1 994] .  A  representative 
comparison  is  shown  in  Figure  30  at  a  tangent  height  of  71.4  km  for  the  dawn  terminator.  The  bulk 
of  the  radiance  at  these  wavelengths  is  from  the  O3  Vj  (OOl)-(OOO)  ("cold")  band.  The  SAMM  and 
CIRRIS  intensities  for  this  band  usually  agree  to  within  a  factor  of  two  or  better.  Near  86  km  (Block 
16)  the  SAMM  cold  band  radiance  is  somewhat  high,  while  near  70  km  (Block  42B)  it  is  too  low, 
especially  on  the  day  side.  This  agreement,  while  reasonable,  is  not  as  good  as  was  found  with  the 
Block  PC12B  radiometer  data  (Figures  26-29). 

The  chemiluminescent  O3  "hot"  Vj  bands  are  the  major  source  of  emissions  between  800  and 
1000  cm  '.  Although  there  are  some  problems  in  the  SAMM  O3  hot  band  model  which  affect 
wavenumbers  below  around  870  cm  '  (wavelengths  above  1 1 .5  pm),  the  overall  agreement  in  both 
spectral  shape  and  intensity  is  good  over  the  full  altitude  and  solar  angle  range  in  the  two-  to 
five-quantum  hot  band  region  (approximately  870- 1 000  cm  ')  that  is  not  overlapped  by  the  CO3  laser 
band  centered  at  961  cm  '.  At  the  resolution  shown,  the  non-LTE  OH  rotational  emission  is  not  very 
prominent  in  the  CIRRIS  data,  and  its  omission  in  SAMM  is  not  a  major  source  of  error. 

On  the  day  side  of  the  terminator,  the  CO,  (0001 1)-(  10001)  and  (0001 1)-(10002)  laser  bands, 
centered  at  respectively  961  cm  '  and  1064  cm  ',  are  underpredicted  by  SAMM  by  a  factor  of  2  to  3, 
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even  accounting  for  the  footprint  effect.  The  same  conclusion  can  be  drawn  from  the  daytime  Filter 
3  and  Filter  7  spectra  in  Section  2.4.1,  which  are  footprint-corrected.  This  observation  is  also 
consistent  with  an  underprediction  of  the  4.3  ^m  (OOOll)-(OOOOl)  band  (see  the  following 
subsection),  which  has  the  same  upper  state. 


Figure  30.  Dawn  Terminator  Comparisons  of  CIRRIS-IA  Data  in  the  8-12  |im  Region  with 
SAMM  Calculations  at  a  Tangent  Height  of  71.4  km  and  a  Solar  Zenith  Angle  of 
95.1°. 

The  greatest  data-model  differences  occur  on  the  high-wavenumber  side  of  the  spectrum  at  high 
altitudes,  where  the  v,  and  2V2-V3  bands  of  ozone  occur.  Near  70  km  the  proportional  differences  are 
not  much  greater  than  for  the  cold  band,  and  appear  to  be  due  mainly  to  an  underprediction  of  O3 
concentration  by  the  atmosphere  model.  However,  at  higher  altitudes  the  CIRRIS  data  show 
significant  excess  radiances;  this  is  more  apparent  on  the  day  side  of  the  terminator,  where  the 
predicted  v,  band  contribution  is  smaller.  The  spectral  shape  of  the  excess  radiance  somewhat 
resembles  that  of  the  2v,-V3  (“crossover”)  band  centered  at  1159  cm"',  although  it  seems  to  be 
broader  and  extended  to  higher  wavenumbers.  While  we  have  no  definitive  explanation  of  this 
excess  emission,  we  surmise  that  it  is  associated  with  ozone  chemiluminescence,  perhaps  arising 
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from  the  crossover  band  in  combination  with  other  bands  and  possibly  containing  a  rotationally  hot 
component. 

2.3.22  4.3  urn  REGION 

CIRRIS-1 A-SAMM  comparisons  of  4.3  pm-region  radiances  near  the  terminator  provide  an 
excellent  test  of  models  for  solar  excitation  of  COj  and  subsequent  energy  transfer  processes  giving 
rise  to  radiation  from  the  Vj-type  bands.  In  contrast  to  the  and  OH  emissions,  the  uncertainty  in 
the  radiance  associated  with  number  density  uncertainties  is  quite  small  due  to  the  optical  thickness 
of  the  bands,  the  uniform  COj  mixing  ratio  over  most  of  the  altitude  range  of  the  4.3  pm  data,  and 
the  relative  unimportance  of  other  trace  species  in  the  excitation  mechanisms. 

Both  the  kinetics  and  radiation  transport  processes  are  extremely  complex,  and  involve  energy 
scrambling  among  a  large  number  of  CO^  vibrational  states  (including  minor  isotopes),  transfer  to 
and  from  N2(v=  1 ),  multiple  absorption  and  re-emission  of  photons  within  and  between  atmospheric 
layers,  and  other  effects.  Considering  this  complexity,  the  ability  of  SAMM  to  reproduce  the 
two-order-of-magnitude  variation  of  this  radiation  across  the  terminator  to  within  typically  a  factor 
of  two  or  better  must  be  considered  a  major  achievement  of  the  code. 

23.2.2.1  4.3  um  In-Band  Radiance  (FP3  Filter  1) 

23.2*2. LI  Methodology 

The  SAMM-CIRRIS  database  [Adler-Golden  et  al,  1994]  contains  several  4.3  pm  radiance 
profile  comparisons  for  daytime  and  nighttime  conditions;  they  show  up  to  a  two  order-of-magnitude 
range  of  variation.  Terminator  data  are  available  in  CIRRIS  only  for  a  limited  range  of  SZAs  and 
altitudes.  For  this  investigation  we  examined  data  from  Block  1  OB,  Mode  34,  which  provide  the  best 
overall,  continuous  SZA  coverage.  The  tangent  height  range  is  65-87  km  and  the  views  are 
southeast  across  the  dawn  terminator.  The  SZA  range,  approximately  92-95°,  extends  from  the  night 
side  to  around  the  middle  of  terminator  and  thus  covers  around  half  (one  order  of  magnitude)  of  the 
day/night  variation. 

As  in  the  O3  and  OH  calculations,  seven  latitude-SZA-local  time-correlated  atmospheres  were 
used  to  simulate  all  of  the  limb  radiances  measurements  from  the  six  CIRRIS  FP3  detectors.  The 
SZA  range  in  the  calculation  is  100°to  87°. 
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2.3.2.2.1.2  Results 


The  measurement-model  comparisons  are  shown  in  Figure  31 .  There  is  excellent  agreement  in 
the  absolute  radiance  levels  (within  about  1 5%)  on  the  night  side  of  the  terminator  over  the  full 
tangent  height  range  shown,  as  well  as  over  the  full  SZA  range  up  to  79  km.  However,  at  the  higher 
altitudes  the  CIRRIS  data  fall  below  the  SAMM  predictions  by  significant  and  increasing  amounts 
(up  to  a  factor  of  1 .6)  in  proceeding  from  the  night  towards  the  center  of  the  terminator.  The  spectral 
analysis  in  the  following  section  ascribes  this  discrepancy  to  CO2  628  isotopic  bands,  which  are 
selectively  solar-pumped  in  the  terminator.  Interestingly,  in  the  non-terminator  profiles  the 
CIRRIS-SAMM  differences,  while  generally  small,  are  in  the  opposite  direction. 

2.3.2.2.2  Spectra 
2.3.2.2.2.1  Methodology 

Since  for  the  4.3  jim  region  no  single  CIRRIS  data  block  contains  spectra  for  a  wide  range  of 
SZAs  and  tangent  heights,  data  were  compiled  from  a  number  of  different  blocks  and  modes 
covering  different  mission  times,  local  times,  and  latitudes  corresponding  to  day,  night,  and 
mid-terminator.  The  terminator  data  are  for  the  dawn  and  are  from  Block  lOB,  Mode  34.  As  before, 
the  SAMM  calculations  were  performed  using  a  range  of  SAG  model  atmospheres  that 
approximately  match  the  latitudes  and  local  times  of  the  data  and  exactly  match  the  SZAs  and 
altitudes.  The  observed  and  calculated  spectra  are  presented  in  detail  in  Adler-Golden  etal.  [1994]. 
Except  for  the  ~66  km  nighttime  spectrum,  from  a  medium-resolution  scan,  the  spectra  are 
high-resolution  (I  cm  ')  and  were  degraded  to  4  cm  '  FWHM  with  the  sech  function.  The  finite  size 
of  the  detector  footprint  has  been  ignored  in  the  calculations,  but  since  at  most  altitudes  and 
wavelengths  the  radiance  does  not  fall  off  steeply  with  altitude  this  approximation  should  be 
satisfactory. 
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4.3  \im  Terminator  (CIRRIS) 


Figure  31.  Dawn  Terminator  Radiometric  Comparisons  of  CIRRIS-IA  Data  (Block  lOB 
Mode  34)  for  the  4.3  pm  Region  (Filter  I)  with  SAMM  Calculations. 


2.3.2,2.2.2  Results 


In  general,  the  results  indicate  that  the  SAMM  calculations  match  the  COj  4.3  |im  spectral 
radiance  to  within  a  factor  of  two  or  better  at  most  wavelengths  over  the  full  two-order-of-magnitude 
range  of  diurnal  variation;  this  represents  a  substantial  modeling  achievement.  The  overall 
agreement  is  best  during  the  daytime.  At  night,  the  overall  radiance  is  underpredicted  by  SAMM  by 
up  to  a  factor  of  three  (at  around  66  km),  a  small  portion  of  which  may  be  due  to  neglect  of  the 
detector  footprint.  At  47  km  and  at  72  km  (see  Figure  32)  and  above  the  nighttime  agreement  is 
excellent;  this  is  in  agreement  with  the  bandpass  radiance  comparison  discussed  in  the  previous 
subsection. 


Figure  32.  Diurnal  Comparisons  of  CIRRIS-IA  Data  in  the  4.3  pm  Region  with  SAMM 
Calculations  at  a  Tangent  Height  of  ~12  km. 

The  comparisons  show  some  data-model  discrepancies  in  the  shape  of  the  4.3  pm  CO2  spectrum. 
The  terminator  radiance  overprediction  by  SAMM  at  high  altitudes  (see  the  previous  subsection)  is 
seen  to  originate  from  the  region  around  2300  cm  '.  The  bands  in  this  region  are  mainly  from  the 
628  isotope,  which  remains  sunlit  further  towards  the  night  side  than  the  main  isotope  and  thus  is 
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preferentially  excited  at  the  terminator.  In  addition,  the  high-wavenumber  edge  of  the  spectrum, 
where  the  (0001  l)-(OOOOl)  ("cold")  band  of  the  main  isotope  is  found,  is  underpredicted  by  SAMM. 
This  is  consistent  with  the  underprediction  of  the  CO^  laser  bands  in  the  8-12  |im  region,  which  have 
the  same  upper  state.  At  the  present  time  we  have  no  explanations  for  these  data-model  differences. 
2.3.2. 3  Other  Wavelength  Regions 

In  addition  to  the  CIRRIS-1 A  terminator  data  discussed  in  Sections  2.4.2. 1  and  2.4.2.2,  in-band 
radiances  are  also  available  for  FP3  filters  2  (4.9-7.0  ^m),  5  (1 1 . 1-12.8  pm),  and  6  (1 7.9-23  pm),  and 
were  briefly  examined.  Filter  2  data  were  taken  on  the  day  side  of  the  dawn  terminator  at 
thermospheric  altitudes  near  the  peak  of  the  NO  profile;  they  show  a  strong  but  gradual  trend  of 
increasing  radiance  with  decreasing  SZA,  as  previously  reported  by  Smith  et  al.  [1992].  According 
to  the  NO  concentration  model  in  SAG,  the  observed  variation  is  only  partly  a  solar  terminator  effect 
per  se;  the  bulk  of  the  increase  is  due  to  the  movement  of  the  spacecraft  towards  higher  latitudes 
(where  the  NO  concentration  is  greater)  as  the  sun  rose. 

The  filter  5  terminator  data  cover  a  spectral  region  within  interferometer  filter  7,  the  data  from 
which  are  discussed  in  Section  2.4.2. 1.2.  Near  the  mesopause  this  bandpass  is  dominated  by 
non-LTE  emission  from  O3  hot  bands,  which  does  not  vary  across  the  terminator,  and  OH  rotational 
emissions,  which  does.  Neither  of  these  radiation  sources  are  completely  imderstood,  and  a  model 
for  the  OH  rotational  emission  is  not  currently  available  in  SAMM.  The  lower  altitudes  can  be 
modeled  reasonably  by  SAMM,  but  there  the  terminator  variation  is  relatively  small  and  is  probably 
due  mainly  to  solar  scattering  by  clouds  and  aerosols,  which  would  need  to  be  specified  correctly  in 
the  model  atmosphere. 

The  filter  6  bandpass  contains  almost  exclusively  H2O  emissions;  the  contribution  from  non-LTE 
OH  rotation  lines  is  small  and  confined  to  altitudes  above  around  80  km.  Neither  the  CIRRIS 

radiance  profiles  nor  SAMM  simulations  show  any  evidence  of  a  solar  terminator-related  variation  in 
this  bandpass. 

2.4  Conclusions  And  Future  Work 

The  CIRRIS- lA-SAMM  comparison  data  base  has  identified  many  categories  of  IR  emission 
processes  which  are  remarkably  well  modeled  by  the  DoD  standard  radiance  and  transmission  codes. 
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However,  the  results  also  indicate  significant  areas  of  uncertainty  in  IR  background  code  predictions 
requiring  new  enhanced  simulation  capabilities  to  support  the  optimization  of  space-based  optical 
sensor  systems.  In  addition,  the  database  comparisons  demonstrate  the  need  for  additional  field 
measurements  of  wider  global  and  wavelength  coverage  to  address  the  important  issues  of 
atmospheric  and  radiance  variability. 

Several  chemistry  and  atmosphere  model  upgrades  are  indicated  by  both  previous  investigations 
and  the  current  CIRRIS-1 A  data  and  SAMM  code  comparisons.  The  first  major  category  of  model 
upgrade  is  based  on  the  discovery  of  a  new  class  of  important  rotational  emission  processes.  New 
NO(v,J)  and  OH(v,N)  models  are  required  to  account  for  the  non-equilibrium  rotational  population 
distributions  producing  substantial  radiance  enhancements  in  SWIR  and  LWIR  bands.  Other 
upgrades  include  a  new  non-auroral  NO^  model  for  accurate  calculations  of  radiance  levels  in  the 
important  4.3  micron  band  used  to  support  optical  sensor  element  design  and  optimization.  The 
comparisons  of  CIRRIS-1  A  and  SAMM  data  for  the  8-11  micron  wavelength  region  indicate 
upgrades  are  required  to  the  kinetics  and  spectroscopy  models  used  to  calculate  Oj  hot  band 
emissions.  Future  SAMM  code  releases  will  also  incorporate  CFC's  and  other  trace  species  which 
are  important  IR  radiators  along  ATH  and  BTH  sensor  lines-of-sight.  Lastly,  to  support  the 
technology  transition  of  earthlimb  and  down-looking  data  from  upcoming  experiments,  including  the 
MSX  program,  to  validated  optical  background  code  simulations  for  a  wide  range  of  atmospheric 
conditions,  a  higher  spectral  resolution  code  capability  is  required. 

A  major  recurring  result  of  the  CIRRI  S-IA-SAMM  data  comparisons  is  the  uncertainty  in  global, 
seasonal,  and  other  atmospheric  background  radiance  variations.  Specifically  for  thermospheric  IR 
radiances,  the  NO  concentration  profile  used  in  the  present  work  is  in  need  of  improvement, 
particularly  for  nighttime  conditions  and  for  low  latitude  locations.  Although  a  simple  model  is 
unlikely  to  completely  account  for  the  relatively  large  variability  associated  with  NO,  an  enhanced 
capability  is  required  to  credibly  bound  the  variability  excursions.  In  addition,  upgrades  to  the  O3 
and  HNO3  atmospheric  profiles  calculated  using  the  SH  ARC/S  AMM  Atmosphere  Generator  would 
significantly  improve  the  agreement  between  measured  and  predicted  radiances  at  altitudes  below  the 
thermosphere.  An  improved  temperature  model  is  also  required  for  the  mesopause  region  to  reduce 
radiance  uncertainties  in  this  altitude  regime.  These  upgrades  and  others  indicated  by  this  study  are 
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required  to  provide  accurate,  validated  atmospheric  IR  background  radiance  codes  for  selecting 
sensor  filter  bandpasses  to  optimize  surveillance,  tracking,  and  interceptor  system  operability. 

The  CIRRIS-1 A-SAMM  solar  terminator  comparisons  have  examined  variations  in  atmospheric 
limb  radiance  in  a  number  of  important  IR  bandpasses  across  the  earth's  shadow  edge  at  tangent 
heights  between  50  and  1 00  km.  The  results,  together  with  those  of  an  earlier  study  [Adler-Golden 
and  Matthew,  1993],  provide  a  reasonably  comprehensive  assessment  of  the  ability  of  SAMM  and 
other  DoD  standard  radiation  codes  to  simulate  these  radiance  variations,  a  capability  that  is  required 
for  the  optimization  of  space-based  optical  sensor  systems.  The  comparisons  indicate  areas  of  good 
agreement,  required  code  upgrades,  and  needs  for  additional  field  measurements. 

The  SAMM  calculations  yield  reasonable  qualitative  predictions  of  the  major  variations  in 
MWIR  to  LWIR  limb  radiance  across  the  solar  terminator.  SAMM  achieves  quantitative  accuracy, 
typically  within  a  factor  of  two  or  better,  in  the  strategically  important  4.3  pm  region,  where  the 
in-band  radiance  can  change  by  a  factor  of  1 .5  with  only  a  1  “change  in  solar  angle  (i.e.,  a  1  “change 
in  latitude  at  the  equator).  While  the  overall  description  of  CO^  4.3  pm  emission  is  good,  certain 
states  and  bands  are  overpredicted  (e.g.,  CO2  628  bands  at  the  terminator)  or  underpredicted  (e.g.,  the 
solar-excited  (0001 1 )  state)  by  the  SAMM  kinetic  models  by  up  to  a  factor  of  three.  Reducing  these 
errors  will  require  a  re-examination  of  the  energy-transfer  kinetics  of  the  CO2  vibrational  states. 

For  O3  Vj  and  OH  fundamental  band  emissions,  which  occur  in  the  CIRRIS  8-12  pm  and  X<  3.7 
pm  bandpasses,  good  in-band  radiance  accuracy  (within  a  factor  of  two  or  better)  is  obtained  at  all 
solar  angles  using  model  atmospheres  based  on  a  1-D  photochemical  model.  Poorer  agreement  is 
obtained  between  70  and  85  km  using  atmospheres  generated  with  the  current  version  of  the  SAG 
code,  due  to  differences  in  the  O3  concentrations.  The  factor-of-four  maximum  errors  may  still  be 
acceptable  for  many  systems  applications,  however.  It  is  recommended  that  an  upgraded  version  of 
SAG  be  developed  with  an  improved  description  of  atmospheric  odd-oxygen  and  odd-hydrogen 
photochemistry.  In  addition,  since  the  bulk  of  the  CIRRIS-1  A  terminator  data  at  these  wavelengths 
are  for  the  dawn,  it  is  recommended  that  additional  measurements  be  performed  at  dusk  to  verify  the 
accuracy  of  1-D  photochemical  model  calculations. 
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Comparatively  poor  data-model  agreement  is  obtained  in  the  most  transparent  (lowest 
background  radiance)  portions  of  the  8- 1 2  pm  region.  This  applies  both  near  1 2  pm  where  non-LTE 
O3  and  OH  emissions  contribute  and  in  the  8-9  pm  region  above  70  km  where  significant  excess 
radiance  relative  to  SAMM  is  observed  that  is  yet  to  be  assigned.  These  conclusions  are  based  on 
both  the  terminator  comparisons  and  the  non-terminator  comparisons  presented  previously.  Further 
analysis  of  the  CIRRIS  data  and  possibly  additional  measurements  will  be  required  to  identify  the 
source  of  the  unknown  emissions  near  8  pm. 

In  the  4.9-7.0  pm,  1 1.1-12.8  pm,  and  17.9-23  pm  bandpasses,  in-band  radiance  variations 
across  the  terminator  are  either  gradual  or  absent  in  both  the  data  and  models.  Development  of  a 
model  for  rotationally  excited,  chemiluminescent  OH  will  be  required  to  simulate  the  terminator 
variation  in  the  1 1-12  pm  region.  The  5.3  pm  NO  fundamental  band  radiance  is  consistent  with 
known  NO  production  and  loss  processes;  variations  can  be  large  but  are  slow  compared  to  the 
solar  variation  across  the  terminator.  Latitude-related  variations  in  NO,  predicted  qualitatively 
by  the  SAG  concentration  model,  can  mask  these  temporal  changes  for  a  rapidly  moving  sensor. 
Aurorally-induced  fluctuations  in  NO  density  and  radiation  2ire  also  known  to  be  quite  large  but 
have  not  been  systematically  studied  with  the  CIRRIS- 1 A  database.  This  is  clearly  an  important 
area  for  future  investigation. 
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3.0  ADDITION  OF  A  CORRELATED-^  CAPABILITY  TO  MODTRAN 

3.1  Introduction 

This  section  describes  the  approach  used  to  add  a  Correlated-A:  (CK)  capability  to  the 
MODTRAN  (Berk  et  al,  1998)  atmospheric  transmission/radiance  computer  code.  The  overall 
objective  of  the  effort  is  to  provide  an  accurate  and  fast  means  for  evaluation  of  the  effects  of 
volcanic  aerosols  and  clouds  on  atmospheric  radiative  heating/cooling  and  on  retrievals  of  species 
concentration  profiles.  These  radiative  transfer  computations  require  coupling  the  effects  of  gaseous 
molecular  absorption  due  primarily  to  water  vapor  and  carbon  dioxide  and  particulate  multiple 
scattering  due  to  aerosols,  ice  crystals,  and  water  droplets.  The  molecular  absorption  band  model 
used  in  MODTRAN  is  not  suitable  for  interfacing  with  standard  multiple  scattering  algorithms.  This 
is  because  the  scattering  models  are  based  on  a  Beer’s  law  description  of  the  molecular 
absorption/transmission  (i.e.  T=exp(-ku))  whereas  molecular  band  models  exhibit  different  and  more 
complex  dependencies  on  the  absorption  coefficient  k  and  absorber  amount  u.  In  order  to  adapt  a 
band  model  approach  for  use  in  scattering  calculations  it  is  necessary  to  express  the  band  model 
transmission  function  in  terms  of  a  weighted  sum  officer's  law  exponential  terms.  Thus,  one  of  the 
fundamental  problems  addressed  in  this  work  is  the  development  of  a  method  for  determining  the 
weighting  factors  and  monochromatic  absorption  coefficients  corresponding  to  the  MODTRAN  band 
model. 

3.2  Overview  of  CK  Method 

For  simplicity,  consider  the  problem  of  determining  the  average  transmittance  for  a 
homogeneous  path  over  a  finite  spectral  interval.  The  generalization  to  inhomogeneous  paths  is 
straightforward  and  is  also  described  below.  The  path  transmittance  can  be  exactly  determined 
through  evaluation  of 

J  (1)2 

T(u)  = -  f  dto  exp(-k((o)u),  ( i  \ 

co2‘toi  ^ 

0)i 

where  to  is  frequency  and  k((o  )  is  the  monochromatic  absorption  coefficient.  The  basis  of  the  CK 
approach  is  that  evaluation  of  T(u)  by  integration  over  frequency  can  be  replaced  by  an  equivalent 
integration  over  the  distribution  of  absorption  coefficient  values  f(k)  in  the  spectral  interval 
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(2) 


00 

T(u)=  |dkf(k)exp(-ku). 

0 

The  difficulty  in  evaluation  of  T(u)  via  Equation  (2)  is  the  necessity  of  first  determining  f(k);  the 
advantage  is  that  k(co  )  is  typically  a  highly  repetitive  function  of  w  (i.e.  there  are  many  nearly 
equivalent  values  of  a  given  k  in  a  spectral  interval)  and  thus  numerical  evaluation  of  Equation  (2) 
requires  far  fewer  grid  points  than  Equation  (1).  Note  that  f(k)=0  above  the  maximum  k^^^  and 
below  the  minimum  k^j^  absorption  coefficients  in  the  spectral  interval.  Several  methods  are 
available  for  determining  f(k).  It  can  be  determined  directly  from  k(Q))  by  binning  the  k's  into 
selected  Ak  sub-intervals.  Another  method  takes  advantage  of  the  fact  that  Equation  (2)  is  in  the 
form  of  a  Laplace  transform;  thus,  f(k)  can  be  determined  from  the  inverse  Laplace  transform  of  T(u) 

f(k)  =  L-'(T(u)).  (3) 

The  utility  of  Equation  (3)  is  it  requires  that  only  T(u)  and  not  kfw)  be  known  a  priori  and  therefore 
can  be  applied  to  any  transmittance  model,  both  line-by-line  and  band  model.  The  distribution 
function  f(k)  is  not  smooth  or  monotonic;  it  generally  consists  of  a  series  of  sharp  spikes  which 
reflects  the  sharp  line  structure  of  k(to).  It  is  computationally  more  convenient  to  work  with  the 
smooth  and  monotonic  cumulative  probability  distribution  function 

k 

g(k)=  Jd'kf(k').  (4) 

0 

Physically,  g(k)  corresponds  to  a  monotonic  re-ordering  of  k(co)  within  the  finite  spectral  interval  co, 
to  (Oj.  It  assumes  the  values  g=0  at  k=k^jn  and  g=l  at  k=kn,ax-  The  transmittance  is  related  to  g  by 

1 

T(u)  =  Jdg  exp(-k(g)u),  (5) 

0 

where  k(g)  is  given  by  the  inverse  of  g(k),  k(g)  =  g  '(k).  In  practice,  the  evaluation  of  Equation  (5) 
proceeds  via  summation  according  to 
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imax  _ 

T(u)=  2]  AgjexpC-kiU),  (6) 

i=l 

where  the  Ag/s  and  icj's  are  the  sub-interval  weighting  factors  and  effective  absorption  coefficients 
we  seek  to  define  for  the  MODTRAN  band  model.  The  maximum  number  of  intervals  i^^^  and  the 
specific  selection  of  Agj  interval  boundaries  (gj.gj.,,)  are  chosen  based  on  consideration  of  a  number 
of  factors  including  computational  speed  and  accuracy,  and  the  altitude  and  spectral  ranges  of 
interest.  For  multi-layer  paths  it  assumed  that  there  is  perfect  spectral  correlation  among  the  sub¬ 
intervals  of  each  layer,  resulting  in  a  total  path  transmittance  given  by 

imax  jmax 

T(u)=  Agjexp(.  kiju),  (7) 

i=l  j=l 

where]  denotes  the  sum  over  layers. 

3.3  MODTRAN  k-Distributions 

The  MODTRAN  band  model  for  a  single  species  (multiple  species  are  discussed  later)  is  based 
on  four  parameters:  (1)  the  integrated  line  strength  S  in  a  spectral  interval  Ato  (A(o=l  cm  '  in 
MODTRAN),  (2)  the  effective  number  of  equivalent  lines  n  (non-integer  values  of  n  are  acceptable) 
in  the  interval,  (3)  the  average  pressure  broadening  Lorentz  line  width  y^,  and  (3)  the  Doppler  line 
width  Yd.  These  parameters  are  determined  directly  from  the  HITRAN  {Rothman  et.  al. ,  1998)  line 
parameter  compilation.  The  molecular  transmittance  model  is  given  by 

Ati>  ^  ^ 

where  Ss,=S/n  is  the  strength  for  a  single  effective  line,  and  W  is  the  Voigt  equivalent  width  of  the 
line  taking  into  consideration  only  the  portion  of  the  line  shape  which  falls  within  the  spectral 
interval.  This  functional  form  for  the  band  model  transmittance  derives  from  the  assumption  that  all 
the  lines  in  the  spectral  interval  are  positionally  uncorrelated  and  thus  accounts  for  line  overlap  in  a 
statistically  average  sense.  The  transmittance  contribution  due  to  the  tails  of  lines  which  originate 
from  outside  the  spectral  interval  are  included  as  a  single  multiplicative  Beer's  law  term. 
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Detailed  formulas  for  the  equivalent  width  W  are  available  elsewhere  {Berk  et.  al,  1998).  As 
discussed  earlier  the  k-distribution  f(k)  for  the  band  model  could  be  determined  directly  from  the 
inverse  Laplace  transform  of  T(u)  given  by  Equation  (8).  Alternatively,  and  the  method  applied 
here,  it  can  be  determined  from  line-by-line  (LBL)  simulations  which  are  physically  constrained  by 
the  same  assumptions  used  to  derive  Equation  (8).  The  use  of  real  time  LBL  simulation  or  inverse 
Laplace  transformation  to  derive  f(k)  during  a  MODTRAN  run  is  much  too  slow  to  be  of  practical 
value;  however,  the  approach  adopted  here  is  to  use  LBL  simulation  off-line  to  pre-compute  a  table 
of  k-distributions  which  can  be  rapidly  accessed  and  interpolated  during  a  MODTRAN  run. 


3.3.1  Monte  Carlo  Line-by-Line  Simulation  of  f(k) 

For  a  given  choice  of  n,  Yl,  Yd>  S=1  cm'^  atm  '  (f(k)  is  directly  proportional  to  S  so  any 
convenient  reference  value  can  be  used)  a  series  of  LBL  absorption  coefficient  spectra  are  generated. 
The  positions  of  the  n  lines  for  each  spectrum  are  selected  at  random  with  uniform  probability  over 
the  spectral  interval.  The  f(k)  for  each  spectrum  is  determined  and  all  the  f(k)'s  are  averaged  with 
equal  weighting  to  determine  the  final  band  model  equivalent  k-distribution.  The  final  f(k)  is  then 
transformed  into  the  more  conveniently  used  k(g)  distribution  as  discussed  earlier.  The  procedure  is 
illustrated  in  Figures  33-34  for  an  example  using  n=5  lines.  Figure  33  shows  several  of  the  many 
LBL  spectra  used  to  determine  f(k).  In  determining  the  f(k)  distribution,  each  trial  spectrum  is 
interrogated  over  a  fine  grid  of  equally  spaced  frequency  points  5co  and  the  k  value  at  each  point  is 
binned  into  a  logarithmic  sequence  of  k-bins 

kj  =  kminlO“j,  (9) 

where  a  controls  the  spacing  of  adjacent  k-bins.  In  this  study  a=0.005  was  used  which  results  in  a 
spacing  of  10‘’°”^=1.012.  The  Monte  Carlo  averaged  k(g)  distribution  is  then  calculated  from 


nj  5co 
HMC 


(10) 


where  Oj  is  the  number  of  binned  k  values  between  kj  and  kj+,,  and  n^c  is  the  number  of  Monte  Carlo 
trial  spectra.  The  function  gj(kj)  is  constructed  from  the  Agj's  by 


gj  =  E^8i’ 


i=l 


(11) 
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which  also  defines  the  final  desired  distribution  function  k/gj). 


WAVENUMBER  (cm‘‘) 

Figure  33.  Sample  Monte  Carlo  Trial  Spectral  Calculations.  The  bottom  spectra  are  offset 
by  factors  of  10"^  and  lO"*  from  the  top  spectrum  for  display  purposes 

The  convergence  of  the  k(g)  distribution  as  a  function  of  the  number  of  Monte  Carlo  trial  spectra 
is  shown  in  Figure  34.  Reasonable  convergence  is  obtained  after  20  trials,  although  some 
differences  are  noted  between  20  and  100  trials.  The  number  of  trials  needed  to  obtain  a  good 
distribution  depends  on  the  number  of  lines  n;  as  a  general  rule  of  thumb  at  least  n^c^l  00/n  trials  are 
required  to  obtain  a  useful  result.  It  is  difficult  to  see  the  convergence  of  k(g)  in  the  region  of  g=l  in 
Figure  34;  this  region  is  shown  in  an  expanded  scale  in  Figure  35.  It  is  seen  that  the  peak  value  of  k 
approximately  doubles  when  20  or  more  trials  are  considered.  This  reflects  the  probability  that  any 
two  of  the  n=5  lines  will  closely  overlap  for  a  given  trial  spectrum.  If  enough  trials  were  considered, 
eventually  overlaps  of  3,  4,  and  5  lines  would  be  observed;  however,  the  cumulative  probability 
interval  Ag  associated  with  these  cases  would  be  insignificant  and  would  have  an  inconsequential 
impact  on  radiative  transfer  calculations. 
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CUMULATIVE  PROBABILITY  (g) 

Figure  34.  Convergence  of  the  k(g)  Distribution  as  a  Function  of  the  Number  of  Monte 
Carlo  Trial  Spectra 
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CUMULATIVE  PROBABILITY  (g) 

Figure  35.  Same  as  Figure  33  But  Showing  the  Convergence  in  the  Region  of  g=l  in  More 
Detail 
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Several  examples  of  k-distributions  and  their  variability  with  n,  Yl,  and  7^  are  shown  in  Figures 
36-38.  Figure  36  shows  the  variation  with  the  Lorentz  line  width  where  it  is  seen  that  the  line  wing 
absorption  coefficients  (i.e.  g<0.8)  are  proportional  to  Yl.  Figure  37  shows  that  in  the  limit  of  a  very 
large  number  of  lines  the  absorption  coefficients  asymptote  to  a  constant  value  for  the  entire  spectral 
interval.  Figure  38  illustrates  the  changes  in  the  k-distribution  when  Doppler  broadening  is  added  to 
a  pure  Lorentz  broadened  line.  Only  the  region  near  the  peak  absorption  coefficients  is  shown  (i.e. 
near  g-1.0)  since  the  line  tail  absorption  coefficient  distribution  is  unaffected  by  the  inclusion  of 
Doppler  broadening. 

A  non-integer  number  of  lines  n.x  (x  denotes  a  fractional  line)  is  treated  by  consideration  of  a 
distribution  consisting  of  n  lines  of  strength  S,  and  a  single  line  of  strength  S^.  The  strengths  are 
constrained  to  both  preserve  the  total  interval  strength  (i.e.  S^S^+S^)  and  the  effective  number  of 

lines  n.x.  These  lines  are  used  in  the  Monte  Carlo  simulation  described  above  to  determine  f(k)  for  a 
non-integer  number  of  lines. 


CUMULATIVE  PROBABILITY  (g) 


Figure  36.  Dependence  of  k(g)  on  the  Lorentz  Line  Width  Yl  for  Yd=0.0  cm  '  and  n=l 
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Figure  37.  Dependence  of  k(g)  on  the  Number  of  Lines  n  for  yL=0.01  cm’*  and  yj,  cm’’ 
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Figure  38.  Effect  on  k(g)  for  Inclusion  of  Doppler  Broadening  for  n=l  and  y^  =0.001  cm’* 
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3.3.2  Determination  of  the  MODTRAN  k(g)  Table 

A  table  of  k(g)  distributions  was  constructed  for  real  time  use  during  a  MODTRAN  run.  Each 
tabulated  k(g)  distribution  corresponds  to  a  specific  set  of  values  for  n,  and  and  a  reference 
interval  strength  of  S=1  cm‘^  atm  '.  The  ranges  of  values  considered  for  these  parameters  were 
chosen  to  be  compatible  with  radiative  transfer  calculations  covering  the  altitude  range  0-70  km,  the 
spectral  range  10-10,000  cm  ',  and  including  all  the  MODTRAN  atmospheric  species.  This 
translated  into  the  ranges  ofO.2-1  O'®  cm'  for  y^,  0.02-10'®  cm'  foryD,and  1-90  form  Each  of  these 
ranges  was  spanned  by  a  logarithmic  scale  of  values  with  a  spacing  factor  of  2'^^.  This  resulted  in 
approximately  1 0,000  tabulated  distributions  where  each  distribution  was  represented  at  33  g  values. 
The  overall  size  of  the  k(g)  data  base  is  1.4  mbytes.  The  g  values  were  selected  so  that  the  k(g) 
distribution  with  the  largest  dynamic  range  in  k  values  would  have  each  decade  of  k  values  covered 
by  approximately  3  (ki,gi)  points.  It  would  be  preferable  from  a  computational  speed  point  of  view  to 
span  each  decade  with  fewer  points;  however,  as  shown  later  3  points  per  decade  is  a  reasonable 
trade  off  between  computational  speed  and  accuracy.  The  largest  dynamic  range  occurs  for  n=l, 
yL=10'®  cm  ',  and  yD=10'®  cm  '.  Most  of  the  g,  points  were  selected  based  on  the  logarithmic  scale 

gi  =  l-10-^(1.46780y,  (12) 

where  i  ranged  from  1  to  30.  Three  additional  points  were  added,  g=l  for  the  peak,  and  g=0.05  and 
g-0. 1  to  resolve  structure  near  g=0.0.  The  selection  of  g  values  for  this  stressing  case  is  illustrated  in 
Figure  39  which  shows  that  the  distribution  spans  12  orders  of  magnitude,  and  in  Figures  40-42 
which  provide  a  more  detailed  look  at  various  portions  of  the  distribution.  Each  (n,yL,yD) 
distribution  is  tabulated  at  the  same  gj  grid  points  as  indicated  in  these  figures.  During  a 

MODTRAN  calculation  the  (ki,g;)  points  for  each  atmospheric  layer  are  found  by  a  three 
dimensional  linear  interpolation  over  n,  y^,  and  y^. 
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CUMULATIVE  PROBABILITY  (g) 


Figure  39.  Selection  of  gi  Table  Values  for  the  k(g)  Distribution  with  the  Largest  Dynamic 
Range  of  k  Values 


CUMULATIVE  PROBABILITY  (g) 


Figure  40.  Same  as  Figure  39  But  Showing  an  Expanded  k  range  near  k„i„ 
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CUMULATIVE  PROBABILITY  (g) 


Figure  41.  Same  as  Figure  39  But  Showing  an  Expanded  k  and  g  Scales 


CUMULATIVE  PROBABILITY  (g) 


Figure  42.  Same  as  Figure  39  But  Showing  an  Expanded  k  and  g  Scales 
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3.3.3  Determination  of  Effective  Sub-interval  Absorption  Coefficients 

We  adopt  the  approach  of  Lacis  and  Oinas  {Lads  and  Oinas,  1 99 1 )  for  determining  an  effective 
absorption  coefficient  kj  for  the  sub-interval  Agj=gi+|-gj  which  preserves  the  average  transmittance 
for  each  atmospheric  layer.  It  is  given  by 


exp(-kiu) 


1 


Agi 


Si+I 

jdgexp(-k(g)u). 

gi 


(13) 


Since  we  have  tabulated  the  k(g)  distribution  over  a  relatively  coarse  grid,  an  assumed  functional 
form  for  the  k  variation  over  Ag  is  required  in  order  to  evaluate  Equation  (13).  An  exponential 
variation  is  assumed 


k  =  kiexp(b(g-gi)). 


(14) 


where  b  is  given  by 

^P(ki+i/ki) 

Agj 


(15) 


This  variation  is  a  reasonable  representation  of  k(g)  as  can  be  seen  in  Figures  40-42  where  the 
tabulated  points  are  connected  with  straight  lines  (i.e.  this  is  equivalent  to  Equation  (14)).  There  are 
differences  of  several  percent  with  the  actual  k(g)  in  some  places;  however,  this  case  is  the  most 
stressing  one.  When  eqs.(14)  and  (15)  are  substituted  into  Equation  (16)  the  final  expression  for  kj 
becomes 


u  bAgj 

where  E,  is  a  first  order  exponential  integral.  From  Equation  (16)  it  is  clear  that  the  effective 
absorption  coefficient  is  a  function  of  the  absorption  column  density  u.  While  this  approach 
preserves  the  average  transmittance  for  each  layer  in  a  multi-layer  calculation,  it  does  not  preserve 
the  total  path  transmittance. 
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In  general,  the  CK  transmittance  T^k  and  the  bandmodel  transmittance  Tb^  for  the  same  values  of 
S,  Yl,  Yd.  and  n  will  not  be  identical.  This  is  illustrated  in  Figure  43  for  a  500  km  horizontal  path  at 
15  km  altitude;  differences  in  the  model  transmittances  as  large  as  approximately  V  10%  are  evident. 
Although  both  transmittance  formulations  are  based  on  the  same  underlying  physical  assumptions 
(i.e.  equal  strength  positionally  uncorrelated  lines),  the  implementation  of  each  involves  different 
numerical  approaches  and  approximations.  In  order  to  maintain  internal  consistency  in  MODTRAN 
between  the  two  approaches,  the  CK  absorption  coefficients  are  adjusted  to  exactly  preserve  the 
bandmodel  transmittance  for  each  layer.  This  is  achieved  by  determining  the  fractional  change  in  the 
iq’s  from  Equation  ( 1 6)  required  to  change  Tck  into  Tbm-  It  can  be  shown  that  the  fractional  change 
5  can  be  related  to  the  transmittance  difference  AT=Tck-Tbm  through  a  Taylor  expansion  of  the  CK 
transmittance  function  Equation  (6)  and  is  given  by 

_  _  2  imax  _  _ 

AT  =  -5  Y,  Agjkiuexp(-kiu)  +  —  Y  ^gikf u^exp(-kiu),  (17) 

i=l  ^  i=l 

where  the  value  of  6  of  interest  corresponds  to  the  negative  root  of  this  quadratic  equation.  This 
approach  does  not  preserve  the  total  path  transmittance  for  multi-layer  paths  as  is  evident  in  Figure 
44  for  a  vertical  path  from  5  km  to  space;  the  largest  differences  are  approximately  ±7%. 

In  summary,  a  new  CK  approach  has  been  developed  and  integrated  into  MODTRAN.  We  are 
currently  using  this  new  capability  to  analyze  infrared  observations  of  clouds  in  the  thermal  and  solar 
multiple  scattering  regimes. 
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Figure  43.  Comparison  of  Computed  Transmittances  for  the  Correlated-k  (CK)  and 

bandmodel  (BM)  Approaches  for  a  500  km  Horizontal  Path  at  15  km  Altitude 
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Figure  44.  Same  as  Figure  43  Except  for  a  vertical  Path  from  5  km  to  Space 
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4.0  VERY  NARROW  BAND  MODEL  CALCULATIONS  OF  ATMOSPHERIC  FLUXES 

AND  COOLING  RATES 

Because  of  the  large  computational  speed  advantage  afforded  by  band  models  over  line-by- 
line(LBL)  models  the  development  and  application  of  band  models  to  the  calculation  of  atmospheric 
radiative  fluxes  and  cooling  rates  continues  to  be  an  active  area  of  research.  A  vvide  variety  of  nanow 
band  model(NBM)  and  wide  band  model(WBM)  approaches  have  been  developed.  Typically  the 
calculational  spectral  resolution,  Aw,  of  the  NBMs  falls  within  the  5<Aw^0  cm'  range,  while  for  the 
WBMs  Aw>100  cm Most  of  the  band  model  approaches  are  based  on  one  of  just  a  few  simple 
analytical  functions  that  approximately  represent  the  average  transmission  of  a  finite  spectral  interval 
for  a  homogeneous  path.  Application  of  these  functions  to  the  atmospheric  problem  requires  further 
approximations  in  order  to  handle  the  inhomogeneous  atmospheric  path  properties  as  well  as  the  effects 
of  spectrally  overlapping  bands  from  different  species.  Direct  application  of  these  band  models  to  the 
calculation  of  atmospheric  fluxes  and  cooling  rates  usually  results  in  poor  to  only  fair  agreement  with 
LBL  calculations.  In  order  to  circumvent  this  problem  the  band  models  are  treated  as  parameterizations 
rather  than  as  first  principles  models.  One  or  more  model  parameters  are  then  adjusted  to  yield  a  best 
fit  to  LBL  reference  transmittance  or  atmospheric  flux  calculations.  These  tuned  band  models  can  yield 
atmospheric  flux  and  cooling  rate  predictions  of  LBL-like  quality. 

A  new  very  narrow  band  model(VNBM)approach  [Bernstein  et  al,  1996]  has  been  developed 
which  yields  atmospheric  flux  and  cooling  rate  predictions  of  near  LBL  quality  without  the  necessity  of 
tuning  the  VNBM  to  a  specific  set  of  LBL  reference  calculations.  The  VNBM  is  incorporated  in  the 
Phillips  Laboratory  MODTRAN  [Berket.  al,  1998]  atmospheric  transmittance-radiance  code  which 
enables  general  specification  of  the  atmospheric  species,  concentration  and  thermodynamic  profiles, 
calculational  spectral  interval,  and  path  geometry,  and  incorporated  into  the  MODTRAN  atmospheric 
transmittance-radiance  code.  The  VNBM  approach  includes  a  computational  spectral  resolution  of  1 
cm  ,  a  single-line  Voigt  equivalent  width  formalism  based  on  the  Rodgers-Williams  approximation 
which  accounts  for  the  finite  spectral  width  of  the  interval,  explicit  consideration  of  line  tails,  a 
statistical  line  overlap  correction,  a  new  sub-layer  integration  approach  which  treats  the  effect  of  the 
sub-layer  temperature  gradient  on  the  path  radiance,  and  the  Curtis-Godson(CG)  approximation  for 
inhomogeneous  paths.  An  additional  and  new  correction  for  exact  and  near-exact  line  overlap  is 
introduced  as  a  modification  to  the  determination  of  the  band  model  line  density  parameter,  1/d.  The 
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standard  two  parameter  CG  approximation  is  used  for  HjO  and  CO2,  while  the  Goody  three  parameter 
CG  approximation  is  used  for  O3.  Atmospheric  flux  and  cooling  rate  predictions  using  MODTRAN  are 
presented  for  HjO  (with  and  without  the  continuum),  COj,  and  O3  for  several  model  atmospheres.  The 
effect  of  doubling  the  CO2  concentration  is  also  considered.  These  calculations  are  compared  to  line- 
by-line(LBL)  model  calculations  using  the  AER,  GLA,  GFDL,  and  GISS  codes. 

The  availability  of  these  LBL  calculations  is  an  outgrowth  of  the  Intercomparison  of  Radiation 
Codes  used  in  Climate  Models  (ICRCCM)  study.  The  results  from  the  ICRCCM  study  indicated  that 
band  model  calculations  of  atmospheric  fluxes  and  cooling  rates  which  included  all  the  atmospheric 
species  demonstrated  the  best  overall  eigreement  with  LBL  calculations.  However,  the  agreement  was 
much  poorer  when  only  single  species  were  considered  and  also  when  calculations  of  the  sensitivities 
of  the  fluxes  and  cooling  rates  to  changes  in  species  concentrations  were  performed.  The  focus  of  this 
study  is  on  the  individual  species  H2O  (with  and  without  continuum),  CO2,  and  O3.  Several  model 
atmospheres  are  considered  as  well  as  the  sensitivity  to  a  doubling  of  the  CO2  mixing  ratio.  Each  of 
these  molecules  challenge  different  aspects  of  an  atmospheric  radiative  transfer  model.  H2O  dominates 
the  clear  sky  tropospheric  atmospheric  cooling  and  its  concentration  can  vary  extremely  rapidly  with 
altitude.  CO2  dominates  the  cooling  above  the  tropopause  and  thus  requires  that  the  effects  due  to  the 
Voigt  line  shape  function  be  reflected  in  the  model.  O3  is  of  interest  because  its  mixing  ratio  increases 
strongly  with  increasing  altitude.  This  is  a  troublesome  behavior  for  the  Curtis-Godson  approximation 
often  employed  in  band  model  codes  to  account  for  the  non-uniform  distribution  of  absorbers  along  the 
line  of  sight  path. 

The  MODTRAN  calculations  are  based  on  a  0-3,000  cm  '  spectral  interval,  a  60  layer  atmosphere, 
and  a  two  angle  quadrature  integration.  It  was  demonstrated  for  H2O,  CO2,  and  O3  that  the  MODTRAN 
predictions  fall  within  the  spread  of  the  LBL  results.  We  estimate,  based  on  previous  comparisons  of 
MODTRAN  to  F  ASCODE  transmittance  calculations,  that  the  computational  speed  of  MODTRAN  is 
in  the  range  of  10  to  100  times  faster  than  that  for  the  LBL  models.  The  unique  combination  of 
computational  speed,  1  cm’'  spectral  resolution,  wide  spectral  coverage,  wide  selection  of  and  ease  of 
adding  new  species,  and  near  LBL  accuracy  afforded  by  MODTRAN  makes  it  an  attractive  alternative 
to  LBL  models  for  applications  in  which  absolute  accuracy  is  not  the  sole  criterion  for  code  selection. 
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5.0  CREATION  OF  A  MODTRAN  INPUT  FILE  FROM  RADIOSONDE 

MEASUREMENTS 

Local  radiosonde  measurements  of  temperature,  pressure  and  atmospheric  absorbers  such  as 
water  and  ozone  can  be  used  in  MODTRAN  to  model  the  local  radiative  environment.  For  this 
purpose  MODTRAN  must  be  run  with  the  MODEL  =  7  option,  enabling  it  to  read  in  radiosonde  data 
from  the  input  file,  called  TAPES.  This  section  describes  the  workings,  requirements  and  I/O,  of  a 
program,  called  SNDTP5,  which  converts  radiosonde  measurements  into  a  TAPES  format. 

The  input  to  SNDTPS  consists  of  radiosonde  measurements  of  pressure,  temperature  and  water 
vapor  content  at  various  altitude  levels.  Typically,  radiosonde  measurements  are  performed  at 
several  thousand  altitudes  ranging  from  the  surface  of  the  earth  to  tens  of  kilometers.  The  number  of 
atmospheric  altitude  levels  acceptable  to  MODTRAN  is  currently  61.  This  number  is  set  by  a 
FORTRAN  parameter  and  can  be  increased  at  will.  However,  if  MODTRAN  is  to  be  executed 
repeatedly,  there  is  a  need  for  a  program  like  SNDTPS  which  "compresses"  data  at  thousands  of 
altitude  levels  into  a  reasonable  number.  This  is  also  determined  by  a  parameter,  and  currently  set  to 
MODTRAN’s  61.  It  also  can  be  increased  at  will;  of  course,  all  FORTRAN  codes  must  be 
recompiled  once  a  parameter  is  changed. 

5.1  Summary  of  Important  Features 

•  SNDTPS  can  be  executed  interactively  or  in  the  batch  mode. 

•  The  radiosonde  data  can  be  in  any  arbitrary  units  permitted  in  MODTRAN.  TAPES  outputs 
are  in  millibars  for  water  and  pressure,  and  Kelvin,  for  temperature. 

•  MODTRAN  default  profile  values  are  written  out  for  altitudes  beyond  the  highest  sonde 
altitude.  The  MODTRAN  default  profiles  are  chosen  based  on  a  least-squares  match  of 
water  and  temperature  values  to  data.  Consequently,  there  are  two  TAPES's  written  out, 
one  corresponding  to  best-match  temperature  and  the  other,  to  H\20.  The  user  also  has  the 
choice  of  smoothly  blending  in  MODTRAN  default  profiles  with  data  at  higher  altitudes. 

•  The  user  may  have  up  to  1 0  layers  of  data  closest  to  the  ground  retained  "as  is"  in  TAPES 

•  The  user  may  discard  noisy/inaccurate  layers  of  data  both  near  the  ground  and  at  the  top. 
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•  In  order  to  maximize  the  number  of  layers  in  the  measurement  region,  the  user  has  control 
over  the  highest  MODTRAN  altitude  in  TAPES,  in  lieu  of  simply  accepting  the  standard 
100  km.  Most  users  probably  will  set  this  altitude  to  70  km.  Of  course,  more  altitudes  in 
the  sonde-region  means  fewer  in  the  upper  atmosphere,  and  vice  versa. 

5.2  General  Description  of  the  SNDTP5  Program 

The  SNDTP5  program  incorporates  and  satisfies  several  common  sense  ideas  and  requirements. 

PROGRAM  EXECUTION 

SNDTP5  is  a  rather  general  program.  Therefore,  it  queries  the  user  for  several  inputs  to  be 
described  later.  There  are  two  ways  to  run  this  program,  interactively  or  in  the  batch  mode.  Helpful 
information  is  displayed  during  interactive  sessions.  Two  aimotated  sample  batch  input  files  are 
included  for  ease  of  editing.  SNDTP5  can  be  executed  in  the  batch  mode  on  PC's  and  UNIX  systems 
by  typing:  "SNDTP5  <  batchfile". 

ALTITUDE  REGIONS 

The  radiosonde  is  often  inaccurate  at  higher  altitudes,  typically  above  20  km  or  so.  The  surface 
measurements  also  appeeir  to  be  unstable  for  the  first  few  altitude  levels.  Therefore,  the  program 
queries  the  user  for  two  altitudes  (SNDBOT  and  SNDTOP),  below  and  above  which  the  data  are 
ignored.  The  region  bounded  by  SNDBOT  and  SNDTOP  is  referred  to  as  the  sonde-region. 

Below  SNDBOT,  there  are  several  choices.  Two  of  these  are  either  to  do  nothing,  in  which  case 
SNDBOT  becomes  the  ground  level,  or  to  input  surface  meteorological  information.  The  other 
choices  are  to  use  MODTRAN  model  profiles  or  to  extrapolate  (via  interpolation)  radiosonde  data 
down  to  0  km. 

We  have  modified  MODTRAN  to  accept  negative  ground  altitudes  for  the  atmosphere  and  lines- 
of-sight  (LOS's)  as  the  the  radiosonde  data  can  actually  start  from  below  the  sea  level. 

Above  SNDTOP,  MODTRAN  builtin  profiles  are  used.  The  choice  of  the  builtin  MODTRAN 
profiles,  of  which  there  are  six,  is  based  on  least-squares-matches  of  radiosonde  data  to  MODTRAN 
profiles.  Two  separate  TAPE5's  are  output  based  on  best  match  to  temperatures  and  to  water  values. 


61 


Depending  on  user  inputs,  the  least-squares  procedure  does  not  necessarily  match  the  entire 
sonde-region  values  to  the  MODTRAN  models.  The  reason  is  that  at  higher  altitudes,  the  radiosonde 
measurements  (particularly,  relative  humidity)  tend  to  be  inaccurate.  Also,  close  to  the  ground,  the 
measurements  probably  reflect  the  local  terrain,  geography  and  weather  characteristics.  Therefore,  a 
subregion  of  the  region  is  used  for  matching  against  MODTRAN  models  which  are  essentially 
averages  in  time  and  space  over  large  areas  of  the  globe.  The  user  is  asked  to  supply  the  bottom 

(SUBBOT)  and  top  (SUB  TOP)  of  a  subregion  contained  inside  the  sonde-region  for  the  least-squares 
match. 

As  mentioned  in  the  above  paragraphs,  MODTRAN  model  profiles  are  used  at  altitudes  above 
SNDTOP.  However,  in  order  to  avoid  abrupt  discontinuities  in  the  atmospheric  quantities,  the  user 
is  queried  to  supply  an  intermediate  value:  SNDMID,  SNDBOT  <  SNDMID  <  SNDTOP.  From 
SNDMID  to  SNDTOP,  a  linear  combination  of  radiosonde  measurements  and  MODTRAN  profiles 
are  used.  At  SNDMID  the  linear  combination  coefficient  is  unity  for  the  measured  data  and  0  for  the 
MODTRAN  values;  at  SNDTOP,  the  coefficients  are  reversed. 

UNITS  AND  CONTENTS  OF  THE  RADIOSONDE  DATA  FILE 

Typically,  radiosonde  water  measurements  are  in  relative  humidity.  However,  SNDTP5  can 
accept  water  content  in  any  of  these  eight  MODTRAN  units:  ppmv,  cm g/kg,  g/m-\  dew  point  in 
K  and  C,  and  relative  humidity  in  percent.  Pressure  values  can  be  in  mb,  torr  or  atm.  The 
temperature  units  can  be  in  C  or  K.  However,  the  TAPES  outputs  of  SNDTP5  will  be  in  units  of  K, 
mb  and  mb  for  temperature,  pressure  and  water  vapor. 

It  is  assumed  that  radiosonde  data  are  stored  in  a  file  each  line  of  which  is  in  this  format: 
altitude,  pressure,  temperature,  water  content,  ...  The  altitudes  can  be  in  either  increasing  or 
decreasing  order.  There  must  be  one  set  of  data  per  altitude  per  line.  Any  additional  information 
present,  for  example,  time,  is  ignored  by  SNDTP5.  The  numbers  are  read  using  list-directed,  that  is, 
unformatted  read.  Sometimes  the  balloon  performing  measurements  bobs  up-and-down  near  the 
highest  altitude.  Therefore,  SNDTP5  stops  reading  data  at  the  first  sign  of  a  decreasing  altitude. 

The  SNDTP5  program  uses  data  without  any  error  correction.  Bad  data  points  in  the  sonde- 
region  are  not  dropped.  If  the  user  suspect  that  a  data  point  is  bad,  s/he  can  either  edit  that  line  out  of 
the  file  or  replace  the  bad  data  by  a  negative  number  less  than  or  equal  to  -999.0.  The  SNDTP5 
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program  will  treat  that  entry  as  a  "missing"  value  and  will  replace  it  by  an  interpolated  value.  The 
user  has  the  choice  of  doing  this  automatically  for  the  first  1 0  layers. 

Throughout,  SNDTP5  uses  linear  interpolation  for  temperature.  Exponential  interpolation  of 
P/T  and  partial  pressure  divided  by  T  is  used  for  atmospheric  pressure  and  water.  Division  by  T 
converts  these  quantities  to  number  densities. 

FILENAMES 

The  user  is  queried  for  a  file  root  name,  which  for  illustration  purposes  in  this  writeup  is  assumed 
to  be  ROOT.  SNDTP5  assumes  that  the  radiosonde  data  resides  in  a  file  called  ROOT.SND.  If 
these  is  no  such  file,  SNDTP5  looks  for  a  file  called,  simply,  ROOT.  The  TAPES  output  files  will 

be  called  ROOT.TP5T  and  ROOT.TP5W.  ROOT.TP5T  is  the  TAPES  which  has  the  MODTRAN 
built-in  profile  values  based  on  the  best  least-squares-match  to  temperatures;  ROOT.TPSW  has 
values  based  on  best  match  to  water  values.  Miscellaneous  information  is  output  to  a  file  called 

ROOT.OUT.  Information  for  plotting  purposes  is  contained  in  the  file  ROOT.DAT.  It  has  both  the 
radiosonde  data  (in  TAPES  units)  and  the  TAPES  information,  separated  by  a  blank  line. 

CONSERVATION  OF  COLUMN  DENSITIES  AND  TAPES  ALTITUDES 

The  most  important  quantities  that  determine  transmittance  and  radiance  are  the  column  densities 
of  species  along  the  LOS.  Therefore,  it  is  important  how  both  the  TAPES  altitudes  and  the  species 
concentrations  at  these  altitudes  are  chosen  to  represent  the  radiosonde  data. 

How  the  TAPES  altitudes  £ire  chosen  is  described  later.  Once  these  altitudes  are  chosen,  the 
species  concentrations  cannot  simply  be  chosen  by  interpolating  radiosonde  measurements.  If  the 
number  of  TAPES  layers  were  not  restricted  to  61  or  so,  interpolation  may  actually  preserve  column 
densities.  But  61  layers  is  not  sufficient  for  this  purpose.  So  the  TAPES  concentrations  and  altitudes 
are  a  compromise  between  conservation  of  vertical  column  density  and  reproduction  of  the  detailed 
humidity  and  temperature  profiles. 

The  water  concentration  at  a  TAPES  altitude  is  determined  from  the  cumulative  water  amount 
which  is  computed  by  adding  up  the  layer  amounts  in  data.  The  adding  process  starts  from  the  top  as 
it  is  much  less  error  prone  that  starting  from  the  bottom.  At  certain  altitudes,  this  method  may  result 
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in  a  TAPES  concentration  which  differs  significantly  from  the  radiosonde  data.  This  is  particularly 
true  if  the  data  changes  sharply  in  a  very  thin  interval.  If  this  occurs,  we  set  the  TAPES  value  to  the 
interpolated  value. 

Pressure  can  be  calculated  similarly.  However,  experience  suggests  that  simple  exponential 

interpolation  of  P/T  is  sufficient  to  conserve  column  density  of  air,  and  the  more  elaborate  scheme  is 
not  necessary. 

Temperature  is  determined  by  assuming  linear  variation  at  all  altitudes. 

DETERMINATION  OF  ALTITUDES  FOR  TAPES 

SNDTPS  selects  TAPES  altitudes  based  on  one  of  three  schemes.  The  MODTRAN  code  has  a 
FORTRAN  parameter  LAYDIM  which  can  be  set  to  accommodate  arbitrary  number  of  altitudes;  the 
code  must  be  recompiled  if  LAYDIM  is  changed.  LAYDIM  is  also  used  in  SNDTPS  and  set  to  the 
current  MODTRAN  value.  If  any  of  the  schemes  comes  up  with  a  number  of  altitudes  which 
exceeds  LAYDIM,  the  program  will  stop  with  an  appropriate  message. 

SCHEME  1  The  TAPES  altitudes  (in  kilometers)  are  individually  read  from  a  file.  The  altitudes 
must  be  in  the  file  one  per  line  and  in  increasing  or  decreasing  order.  The  SNDBOT  and  SNDTOP 

will  be  inserted  in  the  altitude  grid  if  they  differ  from  grid  points.  Altitudes  below  SNDBOT  are 
ignored. 

SCHEME  2  The  TAPES  altitudes  (in  km)  will  be  computed  from  a  file  containing:  EndPtl, 
EndPt2,  Deltal,  EndPtS,  Delta2, ...,  EndPtN,  DeltaN-1 

EndPtl  is  the  earth  surface  and  EndPtN  is  the  top  of  the  atmosphere  (TOA).  Deltal  is  the 
altitude  increment  for  the  region  from  EndPtl  to  EndPt2,  Delta2  is  the  increment  in  the  region  from 
EndPt2  to  EndPt3, ....  The  values  must  be  one  per  line  or  in  the  following  format: 

EndPtl,  EndPt2,  Deltal 
EndPt3,  Delta2 

EndPtN- 1,  DeltaN-2 
EndPtN,  DeltaN-1 

The  SNDBOT  and  SNDTOP  will  be  inserted  in  the  altitude  grid  if  they  differ  from  grid  points. 
Altitudes  below  SNDBOT  are  ignored. 
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SCHEME  3  This  is  the  most  useful  and  automatic  altitude  scheme.  Most  users  are  expected  to 


use  this  method. 

Radiosonde  pressure  and  water  concentration  (in  practice  in  percent  relative  humidity)  arrays  are 
converted  into  mb  (millibar)  units.  Temperature,  pressure  and  H\20  arrays  are  then  examined  for 
local  minima  and  maxima.  In  practice,  pressure  and  temperature  do  not  possess  significant  extrema 
points,  only  the  water  profile  does.  Finding  extrema  array  in  water  is  not  an  exact  science.  We  have 
tried  to  detect  L-  and  elbow-shaped  extrema  in  addition  to  the  usual  ones.  This  involved  some 
experimentation  on  using  thresholds.  The  HjO  data  usually  has  lot  of  bumps  in  it.  Therefore  prior  to 
finding  extrema,  the  data  is  averaged  over  an  width  of,  say,  100  meters. 

All  the  extrema  are  then  combined  to  form  an  array  of  altitudes.  The  extrema  points  are 
candidates  for  TAPES  altitudes.  But  there  may  not  be  enough  of  them.  So,  the  user  is  queried  for  an 
increment  (DELSND)  which  is  then  used  to  form  an  altitude  grid  in  the  sonde-region.  By  combining 
this  grid  with  the  extrema  array,  a  new  altitude  grid  is  formed.  Still  there  may  not  be  enough 
altitudes  in  the  sonde  region.  There  SNDTP5  also  requires  an  input  for  the  minimum  number  of 
altitudes  (NMNSND)  acceptable  in  the  sonde-region.  The  actual  number  of  layers  will  exceed 
NMNSND  if  possible.  The  input  NMNSND  must  equal  or  exceed  {(SNDTOP-SNDBOT)/DELSND 
+  1). 

(As  mentioned,  the  minima  and  maxima  are  found  by  using  thresholds  on  smoothed  raw  data. 
The  thresholds  and  smoothing  windows  are  both  gradually  increased,  if  needed,  until  the  total 
number  altitudes  in  TAPES  equals  the  MODTRAN  limit.) 

Next  the  altitudes  fi'om  the  top  of  the  sonde-region  to  the  TOA  are  chosen.  For  this  purpose,  the 
user  can  either  choose  the  built-in  MODTRAN  altitudes  or  input  individual  altitudes  explicitly. 

The  total  number  of  MODTRAN  input  altitudes  are  now  computed.  The  program  will  then 
output  a  MODTRAN  input  file  if  the  total  number  of  altitude  does  not  exceed  LAYDIM.  If  this 
condition  is  not  met,  the  program  will  terminate  with  an  appropriate  message. 

5.3  SUMMARY  OF  USER  INPUTS 

1 .  Radiosonde  data  units  for  altitudes,  pressure,  temperature  and  H\20. 

2.  Root  name  of  the  radiosonde  data  file. 
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3.  SNDBOT  and  SNDTOP  (optional). 

4.  SUBBOT,  SUBTOP  and  SNDMID  (for  matching  radiosonde  data  with  MODTRAN  profiles,  see 
comment  below). 

5.  TAPES  altitude  scheme. 

6.  File  names  (optional  for  altitude  schemes  1  and  2). 

7.  DELSND,  NMNSND  (for  altitude  scheme  3) 

8.  Options  and  data  (optional)  for  action  below  SNDBOT. 

All  FORTRAN  parameters  reside  in  the  file  PARMS.H.  These  need  not  be  changed  often;  the 
program  must  be  recompiled  if  they  are  changed. 

In  addition  to  inputting  SUBBOT,  SUBTOP  and  SNDMID  explicitly,  there  are  two  other 
choices.  The  user  may  input  "1"  which  sets  SUBBOT  =  SNDBOT,  SUBTOP  =  SNDTOP  and 
SNDMID  =  SNDTOP.  An  input  of  "2"  is  identical  to "  1 "  in  all  respects  except  that  SNDMID  is  one- 
tenth  of  the  way  (one-tenth  of  the  width  of  the  sonde-region)  below  SNDTOP.  An  input  of  "3" 
requires  all  three  inputs  explicitly. 

5.4  TWO  ANNOTATED  SAMPLE  INPUT  FILES 

Both  input  files  are  for  the  altitude  scheme  3  where  the  TAPES  altitudes  are  selected  by 
examining  extrema  of  the  radiosonde  data. 

The  input  file  below  drops  data  below  0.33  and  above  22  km.  Note  that  the  input  line  containing 
TAPEl  cannot  be  annotated  as  this  line  is  read  as  a  character  string.  See  the  previous  section  for 
understanding  the  input  for  defining  the  subregion  for  matching  and  blending  with  the  MODTRAN 
internal  profiles. 

TAPEl 

1  (1  KM,  2  M) 

1  (1  MB,  2  ATM,  3  TORR) 

2  (1  K,  2  C) 

8  (1  PPMV,  2  CM- 3,  3  G/KG,  4  G/M3 ,  5  MB,  6  DP  (K) ,  7  DP  (C) ,  8  RH) 

(DROP  DATA  BELOW  AN  ALT?  1/0.  NEXT  LINE:  ALT  IF  1,  NO.”oF  "AS  IS"  LYRS  IF  0) 

1^  (IGNORE  DATA  ABOVE  AN  ALTITUDE?  1=YES,  0=NO.  IF  1,  NEED  ALTITUDE  BELOW) 

3  (1,2,3  FOR  ALTITUDE  SCHEME,  ISCHEM) 

2.0  (DELSND,  FOR  ISCHEM=3) 

32  (NMNSND,  FOR  ISCHEM=3) 

1  (1,  2,  3  FOR  ALTITUDES  ABOVE  THE  SONDE-REGION,  ITOP) 

70.  (MAXIMUM  TAPES  ALT) 

1  (1,  2,  3  FOR  DEFINING  SUBREGION  FOR  MATCHING  AND  BLENDING  W/  MODTRAN  MODEL) 

I  (OPTION  BELOW  SONDE  l=DO  NOTHING.  2=SURF  DAT,  3=INTERP  TO  0  KM.  4=MODTRAN  MODEL) 
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This  input  file  does  not  drop  any  data.  It  prints  the  three  lowest  layers  "as  is".  Note  that  one  cannot 
have  any  "as  is"  layers  while  at  the  same  time  ignore  data  below  a  certain  altitude. 

TAPEl 

1  (1  KM,  2  M) 

1  (1  MB,  2  ATM,  3  TORR) 

2  (1  K,  2  C) 

8  (1  PPMV,  2  CM-3,  3  G/KG,  4  G/M3,  5  MB,  6  DP  (K)  ,  7  DP  (C)  ,  8  RH) 

0  (DROP  DATA  BELOW  AN  ALT?  1/0.  NEXT  LINE:  ALT  IF  1,  NO.  OF  "AS  IS"  LYRS  IF  0) 

3 

0  (IGNORE  DATA  ABOVE  AN  ALTITUDE?  1=YES,  0=NO.  IF  1,  NEED  ALTITUDE  BELOW) 

3  (1,  2,  3  FOR  ALTITUDE  SCHEME,  ISCHEM) 

2.0  (DELSND) 

40  (NMNSND) 

1  (1,  2,  3  FOR  ALTITUDES  ABOVE  THE  SONDE-REGION,  ITOP) 

70.  (MAXIMUM  TAPES  ALT) 

1  (1,2,3  FOR  DEFINING  SUBREGION  FOR  MATCHING  AND  BLENDING  W/  MODTRAN  MODEL) 

1  (OPTION  BELOW  SONDE  l=DO  NOTHING,  2=SURF  DAT,  3=INTERP  TO  0  KM,  4=MODTRAN  MODEL) 

5.5  Example  Calculations 


Here  we  show  several  plots  of  radiosonde  data  and  TAPES  values  for  a  typical  set  of  radiosonde 
measurements.  The  abscissa  values  of  the  circles  denote  the  TAPES  altitudes.  As  can  be  seen  from 
these  figures,  the  SNDTPS  program  captures  the  shape  of  the  radiosonde  data  reasonably  well. 
These  were  generated  using  the  input  files  shown  above.  Note  that  in  the  radiosonde  data  P  and  T 
become  zero  (0  mb  and  0  K)  suddenly  after  22.8SS  km.  The  relative  humidity  is  also  essentially  zero 
beyond  22.8SS  km.  The  first  three  plots  (Figures  4S-47)  process  the  entire  sonde  region  in  the  data 
file,  from  0.3 1 S  km  to  68  km.  The  first  three  altitudes,  at  0.3 1 5, 0.320  and  0.330  km,  are  retained  "as 
is"  in  TAPES.  Notice  the  sharp  drop  in  the  H\20  value  at  0.320  km. 

In  the  second  set  of  three  plots  (Figures  48-50),  data  below  0.33  km  and  above  22  km  are 
dropped. 
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Figure  45.  Comparison  of  Sonde  Data  and  TAPES  File  for  Temperature.  No  Data  are 
Dropped,  and  the  Number  of  "as  is"  Layers  is  3.  The  Appropriate  MODTRAN 
Model  is  Chosen  by  Matching  the  Entire  Data.  However,  the  Data  Abruptly 
Become  Zero  (0  mb  and  0  K)  Above  22.855  km.  This  Figure  is  Generated  by 
Using  the  Second  Input  File  of  the  Preceding  Section. 


Figure  46.  Comparison  of  TAPES  vs.  Data  for  Pressure,  Generated  as  in  Figure  45. 
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Figure  47.  Comparison  of  TAPES  vs.  Data  for  Water  Vapor,  Generated  as  in  Figure  45 


Figure  48.  Comparison  of  Sonde  Data  and  TAPES  File  for  Temperature.  Data  Below  0.33 
and  22  km  are  Dropped  Because  the  Temperature  Measurements  Fluctuate  for  the 
First  Few  Measurements  and  All  Measurements  Become  Abruptly  Zero  Above  22 
km.  This  Figure  is  Generated  by  Using  the  First  Input  File  of  the  Previous  Section. 
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Figure  49.  Comparison  of  TAPES  vs.  Data  for  Pressure,  Generated  as  in  Figure  48 


Figure  50.  Comparison  of  TAPES  vs.  Data  for  Water  Vapor,  Generated  as  in  Figure  48 
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6.0  IMPROVEMENTS  TO  MODTRAN  ACCURACY  FOR  REMOTE  SENSING 


Real-time  or  near-real-time  optical  remote  sensing  of  the  atmosphere  requires  a  fast  yet  accurate 
computer  code  for  calculating  transmittance  and/or  radiance  at  the  spectral  resolution  of  the  sensor. 
MODTRAN  [Berk  et  al.,  1998],  a  moderate-resolution,  band-model-based  atmospheric  radiation 
code,  has  proven  valuable  for  a  number  of  such  applications,  including  remote  sensing  of 
temperature  and  humidity  [Anderson  et  al.,  1993]  and  the  calculation  of  radiative  heating  and 
cooling  rates  [Bernstein  et  al.,  1996]. 

In  the  current  investigation,  we  compared  MODTRAN3  with  "exact"  line-by-line  calculations  for 
several  representative  remote  sensing  applications.  We  also  developed  improvements  to 
MODTRAN  which  provide  increased  accuracy.  Most  of  the  improvements  have  been  incorporated 
in  a  new  beta  test  code.  Upgrades  were  made  to  the  O2  vibrational  continuum,  the  collisional 
linewidths,  the  equivalent  width  formula,  and  the  band  model  formalism.  The  band  model  has  been 
revised  to  address  three  problem  areas: 

( 1 )  Line  strength  distribution  effects.  A  more  sophisticated,  2-line-group  model  was  developed 
which  accounts  for  effects  of  a  wide  line  strength  dynamic  range  within  a  spectral  bin.  As 
discussed  in  this  report,  the  improvements  in  accuracy  are  significant  for  molecules  such  as 
ozone  that  have  both  a  high  density  of  lines  and  a  wide  dynamic  range  of  line  strengths. 

(2)  Effects  caused  by  non-random  line  locations  within  the  spectral  bins.  Such  positional 
correlation  effects  are  noticeable  in  a  limited  number  of  spectral  regions,  such  as  in  the  NOj 
band  near  1600  cm  '.  A  practical  and  effective  method  for  adjusting  the  band  model 
parameters  to  compensate  for  these  effects  was  developed  and  tested  in  a  stand-alone  code 
for  uniform  paths. 

(3)  The  quasi-continuum  representing  the  line  tails  that  "leak"  into  or  out  of  the  l-cm  '-wide 
spectral  bins.  An  improved  method  for  calculating  this  continuum  has  been  developed. 

The  results  of  this  work  are  illustrated  with  transmittance  and  spectral  radiance  comparisons 
between  MODTRAN3,  the  new  beta  code,  and  "exact"  FASCODE  [Clough  et  al. ,  1 988]  calculations 
for  both  horizontal  and  vertical  atmospheric  paths.  The  overall  agreement  between  the  spectra  from 
all  three  codes  is  found  to  be  quite  good  for  resolutions  down  to  several  cm  '. 
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6.1  Band  Model  Upgrade 

6.1.1  Standard  MODTRAN  Approach 

For  the  discrete  line  portion  of  the  spectrum,  MODTRAN's  calculation  of  transmittance  for  a  given 
species  has  been  based  on  the  simple  power-law  expression, 

T  =  (1  -  W/Av)^'"'*  .  (18) 

This  equation  assumes  that  the  spectrum  can  be  represented  by  a  set  of  n=Av/d  lines  of  identical 
strength,  S,  randomly  located  in  the  bin  of  width  Av  (=1  cm');  d  represents  the  average  line  spacing 
and  W  the  equivalent  width  of  an  individual  line.  W  is  based  on  standard  formulas  applicable  to  an 
infinitely  wide  bin;  an  approximate  correction  for  the  finite  bin  width  is  made.  Following  Goody 
and  Yung  [1989],  the  S  and  n  values  may  be  chosen  by  solving  the  simultaneous  equations 

nS"'  =  (19) 

nS  =  ZSj  (20) 

where  the  Sj  are  the  actual  line  strengths.  This  solution  insures  that  the  sum  of  the  equivalent  widths 

for  the  model  lines  matches  the  sum  for  the  real  lines  in  both  the  optically  thin  (linear  curve  of 
growth)  and  thick-Lorentzian  (square-root  curve  of  growth)  limits.  As  a  result,  the  correct 
transmittance  is  obtained  in  those  same  limits  as  long  as  the  individual  line  equivalent  widths  are 
small  compared  to  the  width  of  the  bin,  Av. 

A  basic  problem  with  this  approach  is  that  there  is  usually  a  wide  range  of  line  strengths  in  a 
given  bin,  so  it  can  be  a  poor  approximation  to  represent  the  actual  lines  by  equal-strength  lines 
*  when  the  all  of  the  lines  are  neither  very  thin  nor  very  thick.  The  problem  is  most  pronounced  at 
high  altitudes,  where  the  lines  are  closer  to  Doppler  than  Lorentzian  in  shape  and  the  limiting  square- 
root  regime  for  which  the  model  was  designed  is  reached  only  at  extremely  high  optical  depths  that 
are  not  encountered  in  the  atmosphere. 

At  lower  altitudes,  where  the  lines  are  nearly  pure  Lorentzian,  the  errors  are  smaller  but  still 
significant  (up  to  several  percent  in  transmittance)  when  the  equivalent  widths  are  intermediate 
between  the  thin  and  thick  limits.  To  find  the  conditions  under  which  this  will  occur,  we  examine 
the  equivalent  width  expression  for  Lorentzian  lines. 
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W  =  Su/(1  +  0.25Su/a)''' 


(21) 


where  u  is  the  column  density  and  a  is  the  pressure-broadened  half-width.  Intermediate  equivalent 
widths  occur  when  the  two  terms  in  the  denominator  of  Equation  (2 1 )  are  comparable  in  size.  The 
effect  of  W  errors  on  fractional  transmittance  is  greatest  when  the  transmittance  is  of  order  1/e  -that 
is,  when  W  is  approximately  equal  to  the  average  spacing  d.  Setting  Su/4a  =  1,  W  =  d,  and  a  = 
0.07P  where  P  is  total  pressure  in  atm  leads  to  the  condition  1/d  =  5/P  under  which  sizable 
transmittance  errors  are  expected.  Since  P  is  always  1  atm  or  less,  this  occurs  only  for  molecular 
bands  having  more  than  around  5  lines  per  wavenumber.  Generally  speaking,  these  are  bands  of  Oj, 
NjO,  NOj,  and  CH4. 

To  reduce  the  overall  errors,  the  1/d  values  may  be  empirically  reduced.  In  MODTRAN3  this  is 
accomplished  by  combining  lines  within  sub-bins.  This  does  a  reasonable  job  of  accounting  for  the 
line  density-dependence  of  the  error  and  gives  good  results  for  the  test  case  of  a  500  km  horizontal 
path  at  15  km.  However,  no  single  1/d  value  is  correct  at  all  optical  depths,  so  substantial  errors 
occur  for  certain  atmospheric  paths.  For  example,  as  shown  in  the  following  subsection 
MODTRAN3  can  give  up  to  a  10%  transmittance  error  in  the  ozone  fundamental  band  at  40  km 
altitude. 

6.1.2  Modified  Approach 

To  address  these  problems  of  the  standard  MODTRAN  band  model,  we  have  developed  a  new 
model  that  uses  two  groups  of  lines  of  different  strengths  to  better  approximate  the  actual  line 
strength  distribution.  The  use  of  four  parameters  (S,,  Sj,  n„  and  nj)  rather  than  two  (e.g..  S/d  and 
1/d)  allows  (and  requires)  the  incorporation  of  two  extra  mathematical  constraints.  The  drawback  is 
that  the  data  storage  and  computation  requirements  are  increased  by  almost  a  factor  of  2. 

In  the  current  investigation,  the  constraints  on  the  parameters  were  taken  as 

n,S,P''  +  n2S2'"  =  ESi'"  ,  p  -  1  -4  (22-25) 

where  the  subscripts  "1"  and  "2"  label  the  line  group  and  the  subscript  "i"  labels  the  actual  line 
strengths.  The  p  =  1  and  p  =  2  equations  correspond  to  Equations  ( 1 9)  and  (20).  The  p  =  4  equation 
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insures  that  the  equivalent  width  sums  for  the  actual  and  model  line  distributions  are  matched 
through  second  order  terms  in  line  strength.  With  the  addition  of  the  p  =  3  equation  there  is  a  simple, 
analytical  solution  which  is  found  to  give  good  overall  agreement  with  calculations  that  use  the  exact 
line  strength  distribution. 

Figure  5 1  shows  comparisons  of  MODTRAN3  and  the  new  algorithm  with  FASCOD3  for  a  500 
km  horizontal  path  at  15  km  altitude.  (For  clarity,  only  the  difference  vs.  FASCOD3  is  shown  for 
the  new  algorithm.)  The  agreement  with  FASCOD3  is  quite  good  overall  for  both  band  model 
codes,  the  RMS  transmittance  errors  are  only  1 .8%.  It  should  be  noted  that  these  comparisons  are 
superior  to  ones  performed  previously  [Anderson  et  al,  1993]due  to  the  Oj  continuum  correction, 
discussed  in  Section  6.4.  The  new  algorithm  provides  noticeable  improvements  in  a  few 
troublesome  regions,  such  as  on  the  red  sides  of  the  CO^  and  O3  bands  near  600  and  1000  cm  ', 
respectively.  The  region  near  600  cm '  is  significant  for  remote  temperature  sensing,  and  is 
discussed  further  in  Section  6.5. 

For  a  1 000  km  path  at  40  km  altitude,  MODTRAN3  gives  good  results  at  most  wavelengths  but 
produces  up  to  10%  errors  in  the  O3  band  near  1000  cm-1,  as  shown  in  Figure  52. 

Initial  calculations  with  the  new  algorithm  gave  a  smaller  but  still  sizable  (around  5%)  error  in 
the  O3  band  in  the  opposite  direction  from  MODTRAN3 .  The  sizable  error  indicates  that  at  veiy  low 
pressures  even  two  line  groups  have  difficulty  in  matching  the  curve  of  growth  generated  by  the 
extremely  wide  range  of  ozone  line  strengths.  The  error  is  reduced  when  the  range  of  S  is  artificially 
restricted  by  combining  lines  that  are  always  optically  thin  in  the  atmosphere.  To  do  this,  ozone 
lines  that  are  weaker  than  a  certain  cutoff  strength  are  combined  into  a  smaller  number  of  lines  that 
have  approximately  the  cutoff  strength;  the  sum  of  the  strengths  is  preserved.  A  reasonable  value  for 
the  cutoff  strength  is  the  Doppler  half-width  divided  by  a  maximum  column  density  of  around  3  atm- 
cm.  The  band  model  parameters  were  recomputed  using  the  combined  lines  and  the  result  is  shown 
in  Figure  52.  The  agreement  with  FASCODE  is  now  reasonable  in  the  ozone  band.  The  RMS  error 
for  the  entire  spectrum  is  around  20%  smaller  than  with  MODTRAN3.  The  introduction  of  the 
strength  cutoff  has  no  effect  on  the  spectrum  at  15  km  altitude  or  below. 
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HORIZONTAL  PATH.  15  KM  ALT. 


Figure  51.  Comparisons  of  MODTRAN3,  the  New  Band  Model,  and  ''Exact"  FASCOD3 
Transmittance  Calculations  for  a  500  km  Long  Horizontal  Path  Through  the  US 
Standard  Atmosphere  at  15  km  Altitude.  1  cm*'  Spectra  were  Convolved  with  a  5 
cm*'  Wide  Slit.  The  Offset  Traces  at  the  Top  Show  the  Difference  vs.  FASCOD3 
for  the  New  Model  (Upper  Trace)  and  MODTRAN3  (Lower  Trace). 

6.1.3  Path  Averaging 

For  treating  non-uniform  paths,  both  MODTRAN  and  the  new  algorithm  use  Curtis-Godson 
averaging  of  the  band  model  parameters  over  the  varying  temperatures  and  pressures  along  the  path. 
Separate  averages  of  all  parameters,  including  the  Doppler  and  pressure-broadened  half-widths,  are 
computed  for  each  line  group  of  each  species.  The  averaging  method  implicitly  assumes  that  the 
spectrum  for  a  given  species  is  comprised  of  the  same  actual  lines  over  the  entire  path,  an 
assumption  which  may  not  hold  over  a  wide  temperature  range.  To  come  closest  to  this  ideal  with 
the  two-group  algorithm,  the  most  closely  corresponding  line  groups  should  be  averaged  together. 
Since  the  two  line  group  solutions  generally  do  not  correspond  to  actual  lines,  there  is  nothing  to 
identify  the  groups  except  their  n  and  S  values.  For  the  present,  we  identify  corresponding  groups  by 
ranking  them  in  order  of  their  S  values.  This  should  be  adequate  for  a  small  temperature  range  and 
when  the  lines  in  one  group  are  consistently  stronger  than  in  the  other. 
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HORIZONTAL  PATH,  40  KM  ALT. 


Figure  52.  Comparisons  of  MODTI^N3,  the  New  Band  Model,  and  "Exact"  FASCOD3 
Transmittance  Calculations  for  a  500  km  Long  Horizontal  Path  Through  the  US 
Standard  Atmosphere  at  15  km  Altitude.  1  cm '  Spectra  were  Convolved  with  a  5 
cm  '  Wide  Slit.  The  Offset  Traces  at  the  Top  Show  the  Differences  vs.  FASCOD3 
for  the  New  Model  (Upper  Trace)  and  MODTRAN3  (Lower  Trace) 


If  one  set  of  actual  lines  dominates  at  low  temperatures  and  another  dominates  at  high 
temperatures,  the  s-ranking  approach  might  not  properly  identify  corresponding  groups,  causing  the 
averaging  method  to  break  down.  The  single-group  MODTRAN  method  as  well  as  the  correlated-k 
method  would  also  break  down  in  this  situation  because  of  the  changing  line  identities.  For  a  wide 
range  of  temperature  along  the  path,  we  expect  that  the  actual  transmittance  will  be  smaller  than 
predicted  due  to  reduced  line  position  correlation.  The  best  way  to  eliminate  this  problem  would  be 
to  adopt  a  more  time-consuming  line-by-line  equivalent  width  formulation,  such  as  is  used  in  the 
SHARC  [Sundberg  et  al,  1996]  and  SAMM  [Sharma  et  al,  1998], 

Spectral  calculations  for  a  non-uniform,  vertical  path  are  presented  and  discussed  in  Section  6.5. 
In  this  example,  the  differences  between  MODTRAN  and  the  2-line-group  algorithm  turn  out  to  be 
qualitatively  the  same  as  for  a  uniform  path. 
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6.2  Compensation  for  Line  Positional  Correlation 

Like  many  other  band  models,  MODTRAN  implicitly  assumes  that  individual  lines  are 
positionally  uncorrelated,  so  that  the  net  transmittance  can  be  approximated  as  a  product  of 
individual  line  transmittances.  Positional  correlation  among  the  lines  (e.g.,  due  to  spin-orbit  near¬ 
degeneracies,  band  heads,  Q  branches,  etc.)  can  increase  or  (occasionally)  decrease  line  overlap  and 
thereby  alter  the  effective  number  of  lines  in  the  interval,  n=Av/d.  Since  the  degree  of  overlap 
increases  with  the  equivalent  widths,  correlation  effects  are  more  pronounced  at  low  transmittances 
than  at  high  transmittances.  To  some  extent,  at  different  transmittances  different  line  spacings  enter 
into  play.  However,  in  the  important  middle  range  of  transmittance  from  around  0.3  to  0.7,  where 
the  transmittance  is  most  sensitive  to  equivalent  width  errors,  the  line  spacings  of  interest  are  on  the 
order  of  the  average  line  spacing,  d. 

Through  extensive  experimentation  with  simulated  spectra,  we  have  developed  an  efficient  and 
reasonably  accurate  approach  for  compensating  for  line  positional  correlation.  The  basic  approach 
was  developed  for  the  standard  Goody  (equal  line  strength)  band  model;  it  involves  adjusting  the  1/d 
value  (which  represents  the  effective  number  of  lines  in  the  spectral  bin)  by  interpolating  between 
the  standard  value,  which  is  used  for  transmittances  near  unity,  and  a  pre-calculated  correlation- 
corrected  value,  which  is  used  for  transmittances  of  around  1/e  and  below.  To  perform  this 
interpolation,  an  approximate  transmittance  must  first  be  calculated  using  the  uncorrected  1/d.  A 
simpler  alternative  that  avoids  the  transmittance  iteration  and  1/d  interpolation  is  to  use  the 
correlation-corrected  1/d  value  at  all  times  and  accept  the  resulting  errors  at  high  transmittances. 
Since  we  have  not  yet  decided  on  the  optimal  form  for  the  interpolation  function,  for  the  present  we 
shall  discuss  only  the  calculation  of  the  correlation-corrected  1/d  value,  1/d^.  With  the  new  2-line- 
group  algorithm,  the  same  1/d  correction  factor  is  applied  to  both  the  n,  and  nj  values. 

Ideally,  1  /d^  should  be  determined  for  each  spectral  bin  by  ( 1 )  generating  an  exact  spectrum  with 
a  transmittance  of  around  1/e  using  an  optically  thick  Lorentzian  line  shape  and  the  actual  line 
strengths  and  locations,  (2)  comparing  the  resulting  transmittance  with  that  obtained  using  randomly 
positioned  (but  otherwise  identical)  lines,  and  finally  (3)  adjusting  1/d  (at  a  constant  line  strength 
sum  s/d)  for  the  random  case  until  the  exact  transmittance  is  matched.  However,  we  have  found  that 
surprisingly  good  results  can  be  obtained  with  the  following,  more  computationally  efficient 
procedure: 
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( 1 )  Assign  each  line  in  the  bin  an  equivalent  width  Wj  that  is  proportional  to  Sj  and  is  scaled  to 
give  a  net  interval  transmittance  T  of  1/e.  The  scaling  factor  may  be  found  by  using  Equation 
(1)  to  solve  for  W,  which  is  taken  as  the  average  equivalent  width. 

(2)  Assuming  that  the  transmittance  line  shapes  Ti(v)  are  square  with  equivalent  width  W,  and 
are  centered  on  the  actual  line  centers,  evaluate  the  net  transmittance  in  the  bin.  Since  the 
line  shapes  are  square,  this  can  be  accomplished  using  an  exact,  analytical  method  which  is 
much  faster  than  a  line-by-line  grid  calculation. 

(3)  The  1/d  correction  factor  (l/d^  to  1/d  ratio)  equals  the  ratio  of  the  logarithms  of  the 
transmittances  calculated  with  the  randomly  located  lines  (i.e.,  T  =  1/e)  and  with  the  actual 
line  centers  (Step  2). 

At  this  time,  we  do  not  have  a  complete  understanding  of  why  this  approach  works  as  well  as  it 
does,  but  we  have  verified  its  success  using  a  variety  of  actual  and  simulated  molecular  line 
distributions.  In  typical  tests,  the  above  definition  of  1/d,  together  with  an  empirical  interpolation 
function  were  found  to  eliminate  the  bulk  of  the  line  correlation  effects  and  reduce  overall 
transmittance  errors  around  a  factor  of  2  for  the  case  of  equal -strength  lines.  In  addition,  it  is  easy  to 
show  mathematically  that  accurate  results  (in  some  cases  exact)  are  obtained  for  a  set  of  equal- 
strength  lines  in  which  some  or  all  of  the  lines  are  multiply  degenerate. 

Since  there  are  many  other  sources  of  error  in  MODTRAN,  positional  correlation  effects  are 
usually  noticeable  in  only  a  few  spectral  regions  where  there  is  very  strong  line  clustering.  The 
largest  effect  that  we  have  observed  is  in  the  NO^  band  at  around  1600  cm  ',  and  leads  to 
transmittance  errors  of  up  to  around  10%.  This  would  impact  the  1 5  km  horizontal  spectrum  ifNO^ 
were  included  in  the  species  list  for  that  test  case. 

Comparisons  of  the  2-line-group  band  model  and  "exact"  line-by-line  calculations  for  0. 1  atm-cm 
of  NO2  are  shown  in  Figure  53.  Without  the  correlation  correction,  the  band  model  is  very 
inaccurate.  Results  from  MODTRAN3  are  slightly  better  but  still  very  poor.  With  the  correlation 
correction,  the  errors  are  reduced  by  around  a  factor  of  3,  yielding  reasonable  agreement  with  the 
exact  calculation. 

At  the  present  time,  the  correlation  correction  has  not  been  incorporated  into  the  beta  test  code; 
we  plan  to  do  at  a  future  date. 
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N02,  Total  P  =  0.13  atm,  Full  Resolution 


Wavenumber  (cm-1 ) 

Figure  53.  Comparisons  of  "Exact"  and  New  Band  Model  (2-Group  and  2-Group 
Correlation-Corrected)  Transmittances  for  NO2  at  1  cm  '  Spectral  Resolution. 
Column  Density  is  2.5  x  10'^  molec/cm^,  Pressure  is  0.13  atm.  Temperature  is  296  K 

6.3  Quasi-Continuum  Treatment 

MODTRAN  uses  a  simple  quasi-continuum  treatment  to  account  for  the  tail  portions  of  the  line 
that  extend  beyond  the  l-cm  '-wide  spectral  bin  into  which  the  line  center  falls.  It  is  assumed  that 
these  tail  portions  are  sufficiently  flat  that  they  may  be  represented  by  an  average,  continuum 
absorption  coefficient.  This  is  a  good  approximation  for  distant  bins,  and  is  adequate  for  the 
neighboring  bins  except  when  the  line  center  is  close  to  the  bin  edge  and  the  optical  depth  is  very 
large,  leading  to  excess  absorption.  In  this  situation,  it  is  sensible  to  shift  the  line  towards  the  center 
of  the  bin  in  order  to  give  a  more  accurate  result  for  the  equivalent  width  sum  over  all  the  bins. 

To  calculate  the  continuum  and  the  center  bin  equivalent  width,  MODTRAN  has  traditionally 
assumed  that  each  line  center  is  at  a  standard  location  in  the  bin  that  is  not  too  close  to  the  edge,  such 
as  at  the  midpoint  or,  in  a  recent  version,  1/4  of  the  way  from  one  edge.  The  first-order  continuum 
may  be  calculated  by  assuming  a  Lorentz  line  shape  at  1  atm  pressure.  A  correction  to  account  for 
non-Lorentzian  line  tails  is  added  for  several  species.  In  the  radiation  transport  algorithm,  this  pre¬ 
calculated  continuum  is  scaled  by  the  column  density  and  the  total  pressure. 
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In  the  current  investigation,  we  used  an  alternative  approach  for  the  continuum  calculation  in 
which  actual  line  center  locations  (within  the  limits  discussed  below)  are  used  and  the  half-width 
term  in  the  denominator  of  the  lineshape  expression  is  ignored  (this  gives  the  correct  line  shape  in 
the  zero  pressure  limit).  To  maintain  consistency  with  the  center  bin  equivalent  width  calculation, 
the  center  bin  continuum  is  adjusted  to  account  for  the  difference  between  the  absorption  coefficient 
calculated  with  the  actual  line  location  and  that  calculated  with  the  standard  line  location,  which  is 
taken  as  1/4  of  the  way  from  an  edge.  (Because  of  this  adjustment,  the  continuum  for  a  bin 
containing  one  or  more  strong  lines  can  sometimes  be  negative;  these  values  may  be  set  to  zero.)  To 
prevent  excessive  absorption  from  lines  that  are  close  to  the  edge  of  the  bin  and  from  neglect  of  the 
half-width  term,  lines  are  relocated  as  needed  to  insure  that  the  distance  from  the  line  center  to  the 
bin  edges  is  a  certain  minimum  value,  d^|„. 

Comparisons  with  exact  calculations  using  synthetic  line  distributions  showed  that  with  a  d„i„ 
value  of  0.15  to  0.20  cm  the  new  continuum  definition  improves  the  transmittance  accuracy  by  a 
small  amount,  typically  2%  or  less  in  individual  1  cm''  bins.  The  improvement  is  noticeable  when 
the  resolution  and  total  pressure  are  both  high  (see  Section  5.5).  At  low  resolution,  results  are 
virtually  identical  to  those  obtained  using  the  continuum  from  MODTRAN3,  since  both  continua 
have  nearly  the  same  average  values. 

6.4  Other  Upgrades 

In  addition  to  the  band  model  improvements,  the  beta  code  incorporates  three  other  upgrades. 
The  first  is  a  rescaling  of  the  Oj  continuum  (i.e.,  the  pressure-induced  component  of  the  fundamental 
absorption  band)  to  match  that  in  FASCOD3.  The  coefficients  are  multiplied  by  0.78  based  on  the 
work  of  Rinsland  et  al  [1989].  The  omission  of  this  factor  in  earlier  versions  of  MODTRAN  had 
been  a  major  source  of  MODTRAN  vs.  FASCODE  differences  for  the  1 5  km  test  case.  The  scale 
factor  is  included  in  the  December,  1 995  version  of  MODTRAN3  that  was  used  in  the  current  work. 

The  second  upgrade  is  a  more  accurate  temperature  scaling  of  the  collision-broadened  half 
widths.  The  actual  temperature  exponents  from  the  HITRAN  line  list  rather  than  a  preset  value  are 
now  used  to  compute  the  273  K  half-widths.  In  addition,  based  on  a  survey  of  the  exponent  data  we 
decided  to  replace  the  "standard"  exponent  value  of  0.5  with  the  value  0.75,  which  was  previously 
used  only  for  COj. 
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The  third  upgrade  is  the  replacement  of  the  Lorentz  equivalent  width  approximation,  Equation 
(21),  with  a  more  accurate  one,  such  as  the  expression  from  the  NASA  Handbook  [Ludwig,  1973] 

W  =  Su/(1  +0.1768(Su/a)'^y''.  (26) 

The  differences  between  the  two  expressions  can  amount  to  several  percent. 

6.5  Water  Vapor  and  Temperature  Sounding  Simulations 

To  see  what  effects  the  various  code  upgrades  might  have  on  ground-based  remote  sensing  of 
water  vapor  cmd  temperature,  we  ran  a  spectral  radiance  case  using  MODTRAN3,  the  new  algorithm 
(which  includes  all  upgrades  except  the  correlation  correction),  and  FASCOD3  for  a  vertical  line-of- 
sight  through  the  US  Standard  atmosphere.  The  spectrum  was  convolved  with  a  3  cm  '  wide  slit. 

Results  for  the  1900  -  2300  cm  '  (4.3  -  5.3  pm)  region,  which  has  been  studied  by  Cipperly 
[private  communication],  are  shown  in  Figure  54.  While  the  average  error  is  40%  smaller  with  the 
new  algorithm  than  with  MODTRAN3,  the  agreement  with  FASCOD3  is  very  good  for  both  band 
model  codes.  Most  of  the  differences  between  MODTRAN3  and  the  new  code  are  due  to  differences 
in  the  quasi-continuum  from  the  water  line  tails.  At  1  cm''  resolution  these  line  tail  effects  are  more 
pronoxmced.  At  lower  resolution  they  nearly  disappear,  and  the  main  difference  between  the  codes  is 
in  the  NjO  band  near  2200  cm'' .  The  new  algorithm  is  more  accurate  than  MODTRAN3  on  the  low 
wavenumber  side  of  the  N2O  band  but  slightly  less  so  on  the  high  wavenumber  side. 

Figure  55  shows  the  corresponding  calculations  for  the  500  -  900  cm''  (11-20  pm)  region,  which 
has  been  studied  for  remote  temperature  sensing.  Here  the  emission  is  dominated  by  CO2  bands.  As 
elsewhere,  the  agreement  with  FASCOD3  is  good  for  both  band  model  codes;  the  average  error  is 
30%  smaller  with  the  new  algorithm.  The  improvement  in  the  700  -  750  cm''  region  is  due  mainly  to 
use  of  the  2-group  band  model,  while  the  improvement  in  the  500  -  600  cm  '  region  is  due  to  the 
quasi-continuum  upgrade. 
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Figure  54.  Comparisons  of  FASCOD3,  New  Band  Model,  and  MODTRAN3  Spectral 
Radiances  for  a  Vertical  Ground-to-Space  Path  Through  the  US  Standard 
Atmosphere.  Spectra  were  Convolved  with  a  3  cm  '  Wide  Slit.  The  Offset  Traces 
at  the  Bottom  Show  the  Differences  vs.  FASCOD3. 


Figure  55.  As  in  Figure  54,  but  for  the  500  -  900  cm  '  Region 
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6.6  Conclusions 


The  main  conclusions  of  this  investigation  are: 

(1)  Improvements  in  the  accuracy  of  MODTRAN3  have  been  achieved  in  a  beta  test  version  of 
the  code,  primarily  by  upgrading  the  line  tail  quasi-continuum  and  Oj  continuum  models  and 
replacing  the  standard  band  model  with  a  more  sophisticated  two-line-group  treatment . 
Over  wide  spectral  ranges  the  RMS  errors  are  reduced  by  up  to  40%  (not  including  the  effect 
of  the  Oj  continuum  modification),  with  larger  improvements  in  accuracy  for  certain 
wavelengths  and  optical  paths. 

(2)  Comparisons  with  FASCOD3  for  typical  remote  sensing  sceneries  show  that  both  the  new 
code  and  MODTRAN3  have  good  overall  accuracy.  At  5  cm’’  resolution  the  average  RMS 
transmittance  errors  are  2%  or  less,  and  average  spectral  radiance  errors  are  in  a  similar 
proportion  to  radiances  at  nearby  optically  thick  wavelengths.  These  results  lend  support  to 
the  use  of  MODTRAN  for  near-real-time  remote  measurement  of  atmospheric  properties. 

(3)  An  effective  method  has  been  developed  and  tested  for  correcting  band  model  parameters  for 
non-random  (i.e.,  correlated)  positions  of  the  spectral  lines.  Errors  due  to  line  correlation  are 
small  at  most  wavelengths  but  were  found  to  be  large  for  NOj. 

(4)  The  MODTRAN  upgrades  developed  during  this  investigation  lead  to  a  better  match  to  exact 
spectral  calculations  in  the  4  -  5  Dm  region,  but  their  effect  on  water  vapor  determinations 
should  be  negligible  since  the  errors  are  small  to  begin  with  and  fairly  random.  There  may 
be  some  improvement  in  temperature  profiles  determined  from  CO2  radiation  at  long 
wavelengths,  where  the  errors  in  the  MODTRAN3  spectra  are  more  systematic.  This  could 
have  a  small  indirect  effect  on  water  vapor  determinations  and  data-model  comparisons  since 
the  relationship  between  absolute  humidity  (which  is  measured  spectrally)  and  relative 
humidity  (measured  by  most  in-situ  sensors)  is  temperature-sensitive. 

Since  the  accuracy  of  an  approximate  model  such  as  MODTRAN  can  depend  strongly  on  the 
spectral  region  and  the  optical  path,  further  comparisons  with  FASCODE  should  be  performed  for  a 
wider  variety  of  remote  sensing  test  cases.  For  best  accuracy,  a  correction  for  line  position 
correlation  should  be  implemented  for  NOj  and  perhaps  other  species.  Finally,  the  increase  in 
computation  time  due  to  the  2-line-group  band  model  needs  to  be  quantified. 
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7.0  NEW  SOLAR  IRRADIANCE  DATA  SETS  IN  MODTRAN 

MODTRAN  requires  accurate  values  of  solar  irradiance  at  the  top  ofthe  atmosphere  for  the  computation 
of  directly  transmitted  solar  irradiance  and  the  solar  contribution  to  total  path  radiance.  Currently, 
MODTRAN  uses  a  table  of  solar  irradiances  computed  by  Kurucz,  consisting  of  values  from  51  to  50000 
cm-J  at  I  cm-*  intervals.  There  is  a  continuing  effort  towards  solar  spectrum  irradiance  measurements.  For 
example,  recently  Woods  et  al  have  examined  the  Upper  Atmosphere  Research  Satellite  (UARS)  solar  UV 
irradiances.  MODTRAN  has  now  been  upgraded  to  include  several  solar  irradiance  data  sets  in  addition  to 
Kurucz  s  calculations.  Additionally,  the  user  now  has  the  ability  to  use  his/her  own  solar  irradiance  file. 
This  feature  will  be  useful  as  new  calibrated  upper  atmosphere  solar  data  become  available. 

The  five  files  now  available  to  MODTRAN  have  irradiance  values  from  5 1  to  -50000  cm'*  in  step  sizes 
of  1  cm-1.  These  files  are  described  very  briefly  below  and  in  more  detail  in  section  7.1. 

•  newkur.dat;  This  is  the  same  as  oldkur.dat  (the  original  Kurucz  data  currently  used  in  MODTRAN) 
but  with  a  deficiency  between  310-340  nm  corrected 

•  oldkur.dat:  This  is  the  original  Kurucz  data  currently  used  in  MODTRAN. 

•  clikur.dat:  This  data  (from  43478  to  1 25000  cm*  * ),  measured  on  Kitt  Peak  and  elsewhere,  is  from 
K.  Chance  [private  communication,  1997];  however,  the  data  are  replaced  by  the  Kurucz  data  where 
O2  lines  occur  and  supplemented  by  the  Kurucz  data. 

•  cebckkundat:  This  is  the  ATLAS  SSBUV  data  from  24603-49986  cm-*,  supplied  by  R.  Cebula 
[private  communication,  1 997].  Below  24603  cm-  * ,  chkur.dat  values  are  used. 

•  thkur.dat:  This  is  the  SOLSPEC  data  from  G.  Thuillier  [private  communication,  1 997]  and  ranges 
from  1 1404-50221  cm"* .  Below  1 1404,  values  from  newkur.dat  are  used. 

For  the  above  files,  in  the  spectral  region  where  data  are  available,  values  at  consecutive  integral 
wavenumbers  are  obtained  either  by  interpolation  or  by  integrating  data  within  ±0.5  cm**  of  the 

wavenumber.  In  what  follows,  MODTRAN's  solar  data  are  discussed  in  some  detail  and  compared  against 
each  other. 

7.1  New  Solar  Irradiance  Files 

7.1.1  Kurucz  Calculations 

This  is  the  original  MODTRAN  solar  irradiance  "data",  but  these  are  actually  the  result  of  extensive 
line-by-line  calculations  performed  by  Kurucz.  (Gail,  we  should  have  a  couple  more  sentences  about 
Kurucz's  method.)  These  values  were  recently  corrected  in  the  3 1 0  to  340  nm  region  using  measurements 
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carried  out  by  an  UARS  instrument  Solar  Ultraviolet  Spectral  Irradiance  Monitor  (SUSIM).  As  the  figures 
below  will  show,  Kurucz's  calculation  generally  holds  up  very  well  against  measurements.  Furthermore,  it  is 
the  most  comprehensive  irradiance  set  available  with  irradiances  at  wavenumbers  all  the  way  down  to  5 1 
cm'l .  The  tabulated  irradiances  have  values  from  51-50000  cm'^  (200  nm  and  greater)  at  1  cm'^  intervals. 
The  corrected  Kurucz  data  is  the  default  MODTRAN  choice. 

7.1.2  Corrected  Chance  Data  -  Kitt  Peak  Measurements 

The  data,  between  800-230  nm  (125000-43478  cm‘l),  consist  of  two  solar  spectra  which  are  joined 
together.  For  wavelengths  greater  than  305  nm  the  data  is  from  the  McMath  Instrument  (Kitt  Peak)  and  was 
collected  by  Kurucz.  The  data  has  been  mostly  corrected  for  atmospheric  absorption  but  some  atmospheric 
effects  still  remained  (see  below).  For  wavelengths  less  than  300  nm  the  data  is  taken  from  Anderson  and 
Hall.  Between  300-305  nm  (33333-32787  cm'l),  the  spectrum  is  a  combination  of  the  two. 

Chance  recalibrated  (in  wavelength)  the  Anderson  and  Hall  data  using  a  set  of  solar  reference  lines.  The 
Kurucz  data  was  wavelength  calibrated  using  a  single  O2  line.  Chance  believes  that  the  Kurucz  calibration  is 
only  good  to  about  0.002  nm,  and  he  plans  to  recalibrate  for  accuracy  to  better  than  0.001  nm.  A  triangular 
scanning  function  with  a  FWHM  of  0.01  nm  was  used  to  interpolate  the  spectrum  to  a  fixed  grid.  Currently 
the  spectrum  goes  from  about  230-800  nm,  but  the  data  exists  out  to  1300  nm  (but  would  require  some 
calibration  work  to  make  available). 

As  the  Chance  data  still  were  not  corrected  for  O2  absorption,  the  relevant  O2  lines  were  replaced  with 
Kurucz  values.  This  is  an  eminently  reasonable  compromise  given  the  close  general  agreement  of  the 
Chance  data  with  the  Kurucz  calculations. 

7.1.3  Cebula  Data  -  ATLAS  SSBUV  Measurements 

This  data,  supplied  by  Richard  Cebula  of  NASA,  one  of  the  coauthors  of  the  paper  by  Woods  et  al,  is  the 
result  Atmospheric  Laboratory  for  Applications  and  Science  (ATLAS)  Shuttle  Solar  Backscatter  Ultraviolet 
(SSBUV)  experiment.  The  SSBUV  data  covers  200.055-406.453  nm  (49986-24603  cm'^)  with  step  size 
varying  from  0.152  to  0.136  nm.  The  uncertainty  in  wavelength  measurements  is  ±0.026  nm.  The  20 
uncertainly  in  irradiance  varies  from  6%  at  200  nm  to  2.8%  at  250  nm  to  2.6%  at  300-400  nm. 

7.1.4  Thuillier  Data  -  SOLSPEC  Measurements 

This  data,  obtained  using  the  SOLSPEC  instrument,  is  from  199.12-876.86  nm  (50221 .0-1 1404.3  cm'^). 
It  is  comparable  to  the  ATLAS  SSBUV  data  but  is  taken  at  a  much  coarser  step  size;  the  step  size  is  0.36  nm 
at  the  lower  wavelength  end  and  changes  to  0.9  nm  at  the  higher  wavelength  end. 
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7.2  Bandpass  Intercomparisons 


The  overlap  of  the  wavelength  regions  for  different  data  sets  suggests  five  wavelength  regions  for 
intercomparison.  Between  876.76-800.00  nm  (interval  1)  the  Kurucz  and  SOLSPEC  data;  between 
800.00-406.45  nm  (interval  2)  the  Kurucz,  Chance  and  SOLSPEC  data;  between  406.45-230.00  nm  (interval 
3)  all  data  sets;  between  230.00-200.06  nm  (interval  4)  the  Kurucz,  SOLSPEC  and  SSBUV  data;  and  finally, 
between  800.00-230.00  nm  the  Kurucz,  Chance  and  SOLSPEC  sets  can  be  compared. 

Table  6  shows  the  bandpass  integrals  in  these  intervals.  These  values  lead  one  to  conclude  that  all  four 
data  sets  are  consistent  with  each  other.  However,  a  more  rigorous  comparison  must  be  made  by  examining 
the  spectral  variations  between  the  data  sets. 


Table  6.  Wavelength  Regions  for  Different  Data  Sets 


877-800  nm 
Interval  1 

800-406  nm 
Interval  2 

406-230  nm 
Interval  3 

230-200  nm 
Interval  4 

800-230  nm 
Interval  55 

KR80 

647 

122 

0.95 

768 

CH 

645 

118 

764 

CB 

121 

1.02 

TH  82 

654 

_ I  _ _ j  y 

124 

1.02 

778 

(KR  =  corrected  Kurucz,  CH  =  corrected  Chance,  CB  =  SSBUV,  TH  =  SOLSPEC) 


73  Spectral  Comparisons 

7.3.1  SSBUV  and  SOLSPEC  Data  Comparison  Below  400  nm 

The  SSBUV  and  SOLSPEC  measurements,  both  of  which  were  carried  out  in  the  upper  atmosphere,  are 
compared  against  each  other  in  the  region  where  they  overlap  (200-407  nm)  in  Figure  56.  The  SSBUV  data 
was  degraded  into  0.5  nm  bins  in  the  lower  wavelength  half  of  the  region  and  into  0.75  nm  bins  in  the  upper 
half.  The  agreement  between  SSBUV  and  SOLSPEC  is  much  better  in  the  shorter  wavelength  region. 

In  the  rest  of  the  comparisons,  we  use  the  ATLAS  SSBUV  (from  Cebula)  data  as  benchmark 
measurements  below  400  nm  and  the  SOLSPEC  data  (from  Thuillier)  as  benchmark  measurements  above 
400  nm  where  the  SSBUV  data  is  not  available.  This  is  consistent  with  Cebula’s  claim  that  the  accuracy  of 
the  SSBUV  data  is  better  at  400  nm  than  at  shorter  wavelengths  and  with  the  fact  that  the  SSBUV  and  the 
SOLSPEC  data  agree  reasonably  well  between  200-400  nm. 
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Figure  56.  Comparison  of  the  SSBUV  and  SOLSPEC  Measurements  in  the  200-407  nm 
Spectral  Region 


7.3.2  Correction  to  Kurucz  and  Comparison  to  SSBUV  Below  400  nm 

As  mentioned  before  the  Kurucz  data  was  deficient  in  the  310-340  nm  region.  This  was  corrected  by 
using  SUSIM  data.  Figure  57  shows  the  old  and  new  Kurucz  data  plotted  against  the  SSBUV  benchmark 
data.  Clearly,  the  newer  Kurucz  set  is  much  closer  to  the  SSBUV  data.  The  same  figure  shows  that  in  the 
3 1 0-340  region,  the  new  Kurucz  set  is  also  closer  to  the  Chance  measurements  than  the  older  Kurucz  set. 
But  when  compared  against  the  SSBUV  data,  the  new  Kurucz  data  reproduces  the  upper  atmospheric 
SSBUV  data  better  than  the  Chance  set;  the  Chance  data  may  still  may  contain  some  atmospheric  effects. 

The  comparisons  in  Figure  57  establish  the  need  to  validate  the  correction  to  older  Kurucz  data  set.  They 
also  show  that  in  the  wavelength  region  shown  in  Figure  57  the  Kurucz  calculation  is  more  reliable  than  the 
Chance  measurements.  This  is  confirmed,  as  we  shall  see,  for  the  entire  wavelength  region  of  MODTRAN. 
See  Figure  58  for  comparison  in  the  200-400  nm  region;  ratio  plots.  Figure  59,  spanning  this  region  nm  also 
show  that  the  corrected  Kurucz  data  agree  well  with  the  SSBUV  data  whereas  the  Chance  data  are  up  to  ± 
20%  off  the  measurements  around  325  nm.  (In  all  the  radiance  plots,  data  was  degraded  such  that  all  the 
curves  have  similar  undulations.  This  entailed  degrading  different  set  of  data  by  triangular  slits  of  differing 
halfwidth.  The  Chance  and  Kurucz  numbers  were  degraded  by  1 .2  nm  in  Figure  56.  The  ratio  plot  of  Figure 
57  uses  ratios  evaluated  at  every  2  nm.) 
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Figure  57.  Comparison  of  the  New  and  Old  Kurucz  Data  with  the  SSBUV  Measurements 
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Figure  58.  Comparison  of  the  New  Kurucz  and  Chance  Data  with  the  SSBUV  Measurement 
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Figure  59.  Comparison  of  the  Ratio  of  the  New  Kurucz  and  Chance  Data  Sets  to  the  SSBUV 
Measurements 


7.3.3  Kurucz,  Chance  and  SOLSPEC  Data  Above  400  nm 

Figure  60  displays  the  Kurucz,  Chance  and  the  SOLSPEC  data  sets.  Figure  61  shows  the  ratios  of  the 
data  sets.  With  the  exception  of  data  around  720  nm,  as  shown  in  Figure  6 1 ,  Kurucz  and  Chance  agree  well 
with  the  SOLSPEC  data.  Around  this  value  the  Chance  data  are  off  by  about  1 5%  from  both  Kurucz  and  the 
SOLSPEC  data.  This  can  be  attributed  to  atmospheric  effects  due  to  H2O. 

7.4  Conclusions 

The  analysis  and  discussion  in  this  section  strongly  support  that  the  best  available  data  below  400  nm  in 
terms  of  resolution  and  accuracy  the  SSBLfV  data.  The  Kurucz  calculations,  after  the  3 1 0-340  nm  correction, 
agree  well  with  this  data  set.  Beyond,  400  nm,  the  only  available  data  set  is  the  SOLSPEC  data.  This  data 
set  is  of  low  resolution,  however.  The  Kurucz  set  also  agrees  with  the  SOLSPEC  data  remarkably  well.  The 
Kitt  Peak  measurements  still  have  some  unaccounted  for  atmospheric  effects.  Therefore,  the  corrected 
Kurucz  data  set  is  chosen  as  the  MODTRAN  default. 
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Figure  60.  Comparison  of  the  New  Kurucz,  Chance,  and  SOLSPEC  Data  Sets 
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Figure  61.  Comparison  of  the  Ratio  of  the  New  Kurucz  and  Chance  Data  Sets  to  the 
SOLSPEC  Measurements 
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8.0  INCORPORATION  OF  NOVAM  INTO  MODTRAN  3.7/4.0 

The  most  recent  compilation  of  the  NOVAM  (Navy  Oceanic  Vertical  Aerosol  Model)  profiles 
offers  a  new  set  of  aerosol  descriptions,  providing  both  optical  and  size  distributions  appropriate 
from  the  shipboard  surface  to  6  km,  covering  the  spectral  range  from  0.2  micron  to  40  microns  at 
relatively  sparse  spectral  resolution.  Since  the  ozone  retrievals  currently  implemented  in  the  UV 
encompass  an  accounting  of  the  aerosol  background  as  stated  above,  the  addition  of  NOVAM 
profiles  to  MODTRAN  was  deemed  critically  important. 

8.1  NOVAM  Code 

SSI  obtained  the  NOVAM  code  from  NRaD  through  S.  Gathman  [Gathman  and  Davidson, 
1993].  This  code  was  released  by  R.  A.  Paulis  under  the  authority  of  J.  H.  Richter,  Oceanic  and 
Atmospheric  Sciences  Division,  Naval  Command,  Control  and  Ocean  Surveillamce  Center,  San 
Diego.  The  NOVAM  code  is  an  upgrade  to  NAM  (Navy  Aerosol  Model)  which  is  already  in 
MODTRAN.  NOVAM  is  based  on  extensive  direct  shipboard  measurements  carried  out  by  several 
different  agencies  specializing  in  the  marine  environment.  The  inputs  to  the  NOVAM  code  are 
radiosonde  data  consisting  of  altitude,  temperature,  pressure  and  relative  humidity  (RH),  and  other 
surface  observation  parameters  such  as  optical  visibility,  wind  speeds  and  surface  IR  extinction 
(1/km)  at  10.6  microns;  not  all  the  inputs  are  required  for  implementation. 

NOVAM  recognizes  three  types  of  meteorological  profiles  characterized  by  existence  or 
non-existence  of  temperature  inversions.  The  cases  are  denoted  numerically:  1 .  No  inversion,  2. 
Two  inversions,  and  3.  One  inversion.  The  wavelength  spectrum  ranges  fi'om  0.2  to  40  microns;  the 
actual  spectral  grid  (in  microns)  is:  0.2, 0.3, 0.3371, 0.55, 0.6943, 1.06, 1.536, 2.0,  2.25, 2.5, 2.7, 
3.0, 3.3923, 3.75, 4.5, 5.0, 5.5, 6.0, 6.2, 6.5, 7.2, 7.9, 8.2, 8.7, 9.0, 9.2, 10.0, 10.591 , 1 1 .0, 1 1.5, 12.5, 
14.8,  15.0,  16.4,  17.2,  18.5,  21.3,  25.0,  30.0,  40.0.  The  model  contains  four  classes  of  marine 
aerosols  with  three  mode  radii  of 0.03, 0.24  and  2.0  microns,  where  the  mode  radius  is  the  "size"  of 
the  most  populous  part  (i.e.,  the  peak)  of  the  distribution  at  the  RH  of  80%.  The  0.03  micron  aerosol 
consists  of  two  classes:  soluble  and  insoluble.  The  other  two  sizes  consist  of  soluble  aerosols  only. 

The  version  of  NOVAM  from  NRaD  outputs  surface  layer  altitudes,  and  the  net  extinction, 
absorption  and  asymmetry  coefficients  by  combining  the  effect  of  all  four  aerosols.  The  output  of 
NOVAM  consists  of  aerosol  size  distribution  pareimeters,  and  total  extinction,  absorption  and 
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asymmetry  values  as  a  function  of  wavelength.  In  this  study,  NOV  AM  was  modified  to  output  this 
information  as  a  function  of  wavelength  for  a  series  of  altitude  values  beginning  at  the  lowest 
“significant”  radiosonde  altitude  (usually  a  few  meters),  extending  into  the  lower  troposphere.  The 
NOV  AM  model  is  claimed  to  be  valid  up  to  6  km.  However,  in  consultation  with  Gathman  [private 
communication,  1 998],  we  have  restricted  the  NOVAM  aerosol  profiles  to  reach  no  higher  than  2km. 

The  set  of  NOVAM  routines  consists  of  about  6000  lines  of  FORTRAN  code  written  in 
non-standard  FORTRAN  77.  NOVAM,  however,  needs  only  minimal  modification  so  as  to  be 
acceptable  to  most  FORTRAN  compilers.  Extensive  modification  of  the  code  was  ruled  out  in  order 
to  maintain  an  easily  discernible  correspondence  between  the  modified  and  original  versions. 

The  user  should  familiarize  herself/himself  with  the  NOVAM  input  files  of  which  there  are  three: 
(i)  the  Surface  Observation  Data  File,  (ii)  the  Radiosonde  Profile  File,  and  (iii)  a  file  called 
novamdn.  For  purposes  of  familiarizing  with  NOVAM,  it  is  highly  recommended  that  the  user 
consult  the  above  referenced  NOVAM  manual.  In  this  report  only  a  very  brief  description  of  the 
inputs  and  output  are  given.  Questions  regarding  NOV  AM's  use  in  MODTRAN  should  be  directed 
to  the  authors  of  this  report. 

Note  that  the  NOVAM  code  supplied  with  this  delivery  has  13  inputs  in  the  Surface 
Observation  File  as  opposed  to  9  as  stated  on  page  9,  Table  2,  of  the  NOVAM  manual.  These 
inputs  are  the  same  as  stated  for  positions  1  to  7.  The  revised  Table  2  (Table  7)  is  described  below. 
Values  outside  the  stated  range  make  the  code  use  built-in  default  values.  It  is  suggested  that  the 
user  employ  the  default  values  when  any  of  the  specific  data  items  are  not  available. 


Table  7.  Inputs  for  the  NOVAM  Surface  Observation  File 


1 .  Sea  Surface  Temperature  (^C) 

2.  Air  Temperature  (^C) 

3.  Relative  Humidity  (%) 

4.  Optica]  V isibil ity  (km) 

5.  Current  Real  Wind  Speed  (m/s) 

6.  Averaged  Wind  Speed  (24  hours,  m/s) 

7.  Air  Mass  Parameter  (1  to  30) 

8.  Cloud  Cover  Fraction  (0  to  1 ) 

9.  Cloud  Type  (0  to  9) 

10.  Surface  IR  Extinction  at  10.6  micron  (I/km,  0.001  to  100.0) 

1 1 .  Weather  (0  to  11) 

12.  Height  of  Lowest  Cloud  (meters,  negative  value  uses  default) 

13.  Zonal/Seasonal  Category  (1  to  6) 
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The  Radiosonde  Profile  Data  File  is  in  either  of  the  formats  described  on  page  15,  Table  4  and 
Table  5,  of  the  NOV  AM  manual.  Table  4  contains  data,  each  line  of  which  consists  of  an  altitude 
(m),  potential  temperature  (°C)-  and  aerosol  mixing  ratio  (g/kg).  The  relationship  between  the 
potential  temperature  (0)  and  the  usual  air  temperature  (T)  is  given  by  the  formula: 

0  =  T.  (Po/Pr,  K  =  (Cp-C,)/Cp  =  0.288 

where  the  C's  are  heat  capacities  at  constant  pressure  and  constant  volume,  Pq^IOI  3.25  mb  and  both 
temperatures  are  in  Kelvin.  Potential  temperature  is  the  temperature  attained  by  air  at  pressure  P  and 
temperature  T  where  it  is  brought  adiabatically  (i.e.,  at  constant  entropy)  to  a  standard  pressure  P^ 
[Houghton,  1986].  Table  5  contains  data,  each  line  of  which  consists  of  a  line  number  (an  integer), 
log  (base  1 0)  of  pressure  in  millibars  multiplied  by  1  O'*,  the  air  temperature  in  °C,  RH  in  percent  and 
pressure  in  millibars  multiplied  by  10.  As  stated  above,  one  needs  the  profile  data  either  in  the 
format  of  Table  4  and  Table  5.  Table  4  is  said  to  be  in  ‘n’  format  whereas  Table  5  is  said  to  be  in  ‘r’ 
format,  presuming  that  ‘n’  denotes  ‘number’  defined  by  mixing  ratio,  while  ‘r’  denotes  ‘relative 
humidity.’ 

In  addition  to  these  files,  NOV  AM  needs  another  file  called  novam.in.  An  example  of  novam.in 
is  reproduced  below: 

1905 sops 

1905proftxt 

n 

Here,  1905sops  is  the  Surface  Observation  File  and  1905prof.txt  is  the  Profile  File  in  the  ‘n’ 
format  as  indicated  by  the  last  line.  This  file  then  specifies  for  the  program  where  the  necessary  data 
files  can  be  found. 

The  output  of  NOV  AM,  novam.out,  now  in  a  form  suitable  for  MODTRAN,  typically  looks  as 
follows.  The  italicized  text  will  not  appear  in  the  output.  The  first  number  is  40  which  is  the  number 
of  wavelengths  (in  microns)  which  are  then  individually  listed.  The  number  10  is  the  number  of 
altitudes  (in  meters)  which  are  then  individually  listed.  Then  the  temperatures  (in  K)  for  each 
altitude  are  listed,  followed  by  the  pressures  (in  mb)  and  relative  humidity  (RH  in  %).  Then  for  the 
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first  wavelength  (0.2  micron),  the  extinction  coefficients  (in  1/km)  for  each  altitude  is  listed.  The 
absorption  coefficients  (in  1/km)  for  each  altitude  is  followed  by  the  asymmetry  parameters  for  each 
altitude.  Then  the  same  set  of  information  of  the  second  wavelength  (.3  micron)  are  listed.  This 
pattern  continues. 


40  (number  of  wavelengths  and  wavelengths  in  microns) 


.2000 

.3000 

.3371 

.5500 

2.2500 

2.5000 

2.7000 

3.0000 

5.5000 

6.0000 

6.2000 

6.5000 

9,0000 

9.2000 

10.0000 

10.5910 

15.0000 

16.4000 

17.2000 

18.5000 

10  (number  of  altitudes  and  altitudes  in  m  ) 

20.9  123.6  226.3  329.1  393.8  458. 

(temperature  in  K) 

287.65  286.49  285.57  284.85  285.37 

(pressures  in  mb) 

1010.70  999.40  988.10  976.80  969.66 

(RH) 

88.80  91.41  95.39  95.60  81.88 


(spectral  data  for  0.2  microns) 


. 156E+00 

. 146E+00 

.145E+00 

.145E+00 

.377E-01 

.377E-01 

. 377E-01 

(Mtlnction) 

.224E-03 

. 140E-03 

.133E-03 

.132E-03 

.  635E-06 

. 635E-06 

. 635E-06 

(absorption) 

. 801E+00 

.798E+00 

.797E+00 

.797E+00 

.758E+00 

. 758E+00 

.758E+00 

(asynaaetry) 

(spectral  data  for  0.3  microns) 

.150E+00 

.140E+00 

. 139E+00 

. 139E+00 

.283E-01 

.283E-01 

.283E-01 

.377E-05 

.255E-05 

.245E-05 

.243E-05 

.488E-06 

.488E-06 

. 488E-06 

.804E+00 

.800E+00 

.799E+00 

. 799E+00 

.777E+00 

. 777E+00 

.777E+00 

.6943 

1.0600 

1.5360 

2.0000 

3.3923 

3.7500 

4.5000 

5.0000 

7.2000 

7.9000 

8.2000 

8.7000 

11.0000 

11.5000 

12.5000 

14.8000 

21.3000 

25.0000 

30.0000 

40.0000 

6  523.4 

572.0  620 

.7  669.3 

285.95 

285.65  287 

.65  288.91 

288.45 

962.55  955.50  949. 

60  943.73 

937.90 

66.69 

65.60  50. 

08  37.44 

35.80 

. 144E+00 

f 

. 142E+00 

. 140E+00 

.130E-03 

1 

. 128E-03 

.125E-03 

.797E+00 

.797E+00 

.797E+00 

.137E+00 

.135E+00 

.133E+00 

.240E-05 

.236E-05 

.233E-05 

.799E+00 

.799E+00 

.799E+00 

94 


8.2  Incorporation  Into  MODTRAN 

First  all  structure  variables  were  eliminated  and  all  non-standard  system  routines  (such  as  gettim) 
were  also  eliminated  from  NOV  AM.  Several  non-standard  (i.e.,  non-FORTRAN  77)  features  were 
left  intact.  These  include  the  DO ...  ENDDO  structure,  longer  than  six  character  variable  names  and 
the  use  of  the  INCLUDE  statement  as  these  are  acceptable  by  almost  all  modem  compilers.  The 
goal  was  to  minimize  changes  to  NOV  AM  and  to  use  it  almost  "as  is".  The  changes  to  the  NOV  AM 
code  are  briefly  stated  later. 

Extensive  changes  to  MODTRAN  code  had  to  be  made  to  accommodate  the  way  NOV  AM  treats 
its  four  aerosols.  The  reason  changes  were  extensive  is  that,  unlike  MODTRAN's  current 
requirement,  NOV  AM  does  not  output  an  aerosol  profile  (varying  with  altitude)  and  spectral 
extinction  and  absorption  coefficients  (varying  with  wavelength  but  not  with  altitude).  Instead 
NOV  AM  outputs  both  altitude  and  spectrally  varying  quantities  which  are  products  of  profile  and 
spectral  parameters.  Changes  to  NOV  AM  code  itself,  however,  were  kept  to  a  minimum.  This 
meant  that  in  order  to  use  NOV  AM  in  MODTRAN  the  user  must  supply  the  required  radiosonde 
input  data  to  NOV  AM,  separate  from  the  MODTRAN  inputs.  NOV  AM  is  executed  off-line  and 
creates  a  file  called  novam.out  (lower  case  in  UNIX)  which  is  used  as  input  to  MODTRAN 
(uppercase  filename  in  UNIX).  Note  that  NOV  AM  input  files  are  currently  separate  and  in  addition 
to  MODTRAN's  usual  input  file  (which  is  named  tapeS).  If  the  altitudes  in  tapeS  overlap  with  those 
in  the  NOV  AM  output  file,  the  meteorological  parameters,  such  as  humidity,  pressure  and 
temperature,  used  by  MODTRAN  will  be  those  provided  by  NOV  AM. 

In  a  future  upgrade,  the  requirement  for  NOV  AM  to  have  a  separate  input  file  can  be  eliminated; 
both  MODTRAN  and  NOV  AM  will  then  use  the  information  contained  in  the  MODTRAN  input 
file,  tapes.  This  process  will  be  facilitated  by  the  prior  development  of  a  sonde  compression 
scheme;  SSI  and  PL/GPO  have  collaborated  to  write  a  program,  called  SNDTP5,  which  can 
compress  radiosonde  measurements,  consisting  of  hundreds  of  altitude  layers  (such  as  those  used  by 
NOV  AM),  into  a  form  more  suitable  for  the  finite  layering  appropriate  (and  generally  just  as 
accurate  for  transmittance  and  radiance  calculations)  for  a  MODTRAN  tapeS. 

As  mentioned  NOV  AM  actually  can  model  altitudes  as  high  as  6000  meters.  However,  in 
consultation  with  E.  Shettle  [private  communication,  1998]  and  S.  Gathman  [private  communication. 
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1998],  the  maximum  NOVAM  altitude  relevant  for  MODTRAN  was  determined  to  be  2  km.  In 
reality,  for  most  applications  it  will  be  less  than  2  km.  NOVAM  distinguishes  between  three 
different  temperature  inversion  cases.  The  code  was  modified  to  output  these  inversion  layers 
explicitly  which  are  then  used  in  MODTRAN.  This  enables  MODTRAN  to  use  only  a  few  layers 
and  still  accurately  model  the  temperature  effects.  If  the  aerosol  does  not  contain  inversion  layers, 
currently  MODTRAN  will  introduce  layers  which  are  at  most  100  m  apart.  Although  adequate,  this 
scheme  may  be  improved  so  that  layers  are  more  closely  spaced  nearer  to  the  surface  (where  the 
scale  height  is  smaller/steeper)  and  are  farther  apart  towards  the  top  of  the  boundary  layer  (where  the 
scale  height  is  generally  larger).  This  may  allow  using  fewer  layers  without  loss  of  accuracy. 

In  summary,  NOVAM  is  simply  used  to  generate  a  database  of  marine  aerosol  profiles  and 
spectral  information  for  MODTRAN.  NOVAM  does  not  at  present  generate  angular  phase 
functions.  Instead,  it  has  a  database  of  asymmetry  parameters  from  which  Henyey-Greenstein  phase 

functions  can  be  computed.  In  principle,  a  Mie  code  can  be  used  to  generate  the  phase  functions  for 
NOVAM. 

8.3  Some  Results 

Three  typical  (as  provided  in  the  NOVAM  package)  profiles  of  aerosol  extinction  and  coincident 
temperature  are  shown  if  Figures  62a  and  62b.  Figure  63a,  b,  and  c  shows  the  simulated 
backscattered  UV  signatures  associated  with  these  profiles,  as  might  be  measured  by  a  potential 
ozone  monitor  staring  down  from  a  space  platform.  These  calculations  use  all  three  types  of 
temperature  inversions  modeled  in  NOVAM.  The  calculation  with  no  aerosol  includes  only  the 
Rayleigh  scattering  component  and  is  used  as  the  measure  of  change  imparted  to  the  backscattered 
signature  by  low-lying  aerosols.  No  attempt  was  made  to  smoothly  incorporate  these  profiles  into  a 
total  profile.  Rather,  the  "default"  US  Standard  temperature,  pressure,  and  constituent  (primarily 
ozone)  profiles  and  background  rural  (23  km  visibility)  aerosols  were  employed  above  0-2  km,  the 
acceptable  vertical  range  for  the  NOVAM  input.  The  spectral  range  presented  is  only  that  reaching 
the  surface  and  near-surface,  as  wavelengths  shortward  of  300  nm  will  be  absorbed  (in  general)  at 
higher  altitudes.  MODTRAN  will  accommodate  simulations  from  200  nm  to  the  far-IR.  including 
the  aerosol  impact,  so  the  short  spectral  range  depicted  in  these  calculations  is  not  a  restriction. 
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Figure  62a  and  62b.  The  3  aerosol  and  coincident  temperature  profiles  (in  extinction  at  0.3  um 
and  K,  respectively)  as  a  function  of  altitude.  These  profiles  were  chosen  to  capture  the 
number  of  temperature  inversions  used  as  a  parameter  in  NOV  AM,  1  or  no  inversion,  2  or 
two  inversions,  and  3  for  1  inversion.  There  was  no  attempt  to  find  the  most  perturbing  case, 
so  these  can  be  considered  typical.  Note  that  MODTRAN  merges  these  profiles  into  those 
describing  the  rest  of  the  atmospheric  profile  from  whatever  source  has  been  specified, 
‘default’  or  ‘user-defined’.  This  can  lead  to  very  coarse  discontinuities  whose  impact  might 
need  to  be  further  explored. 
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Figure  63a,  b,  c.  As  denoted,  these  represent  typical  sensitivity  to  the  new  NOV  AM  aerosol 
profiles  shown  in  Figure  62.  They  plots  are  shown  linearly  to  emphasize  the  impact  at  the 
longer  wavelengths  that  see’  to  the  surface,  and,  therefore,  would  be  impacted  by  boundary- 
layer  variability.  At  shorter  wavelengths,  <0.3  um),  the  stratospheric  aerosol  component  might 
be  important  under  extremes  of  volcanic  loading.  That  sensitivity  requires  a  logarithmic  plot 
and  has  not  been  explored  in  this  study. 


8.4  Novam  Input  And  Modtran  Input  Files 


The  NOVAM  files  were  described  earlier.  So  they  are  not  reproduced  here.  In  the  delivered 
code,  there  are  several  novam.in,  Surface  Observation  and  Radiosonde  Profile  files. 
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The  tapes  used  to  run  MODTRAN  with  NOVAM  aerosols  for  the  calculations  in  this  report  is 
shown  below.  The  'N'  in  the  third  line  invokes  the  NOVAM  aerosol  option  in  MODTRAN. 


6 

2 

2 

1 

0  0 

0  0  0 

0 

1  0 

OF 

0 

IN 

0 

0 

0 

0  0 

0.000 

.000 

.000 

0050.1 

0000 

.10 

180.00000 

2 

2 

0 

0 

45. 

60. 

.3 

.4 

.0001 

.0010  $ 

Ml 

1  .0500 


.000 


0 


First  NOVAM  is  executed  to  produce  the  novam.out  file.  This  file  then  should  be  copied  to  the 
directory  containing  the  MODTRAN  executable  as  NOVAM.OUT;  MODTRAN  requires  this  file 
with  the  uppercase  name. 


8.5  Future  Upgrades  To  Novam  Implementation 


There  are  at  least  six  general  areas  in  which  the  aerosol  product  in  MODTRAN  can  be  improved. 

1 .  The  first  is  to  enable  NOVAM  to  run  from  MODTRAN's  input  file,  tapeS.  This  task 
will  enable  MODTRAN  to  use  radiosonde  data  consisting  of  several  hundred  altitude 
layers  several  of  which  can  even  be  redundant.  This  will  alleviate  the  need  for 
NOVAM  to  have  its  own  input  file  as  is  required  in  the  current  input  scheme.  Note 
that  there  still  may  be  a  need  for  the  NOVAM  input  file,  for  example,  to  input  surface 
observations. 

2.  MODTRAN  does  not  now  have  phase  functions  for  several  aerosols  (e.g.,  the  desert 
aerosols)  and  for  none  of  the  cloud  models.  In  the  future  this  can  be  rectified  by 
generation  of  the  phase  functions  using  the  Mie  code  and  incorporating  them  in 
MODTRAN. 

3.  The  phase  functions  for  NOVAM  are  also  not  available.  In  consultation  with  S. 
Gathman,  they  can  be  generated  for  the  NOVAM  aerosols  and  incorporated  in 
MODTRAN. 

4.  The  output  of  the  Mie  code  can  be  put  in  a  format  so  that  user  can  include  them  in  the 
MODTRAN  input  file  without  extensive  editing. 

5.  Based  on  the  El  Chichon  and  Mt.  Pinatubo  erruptions,  the  content,  size,  type,  and 
H2SO4  component  of  fresh  and  aging  volcanic  aerosols  need  to  be  altered  fi’om  the 
default  profiles  now  available  within  MODTRAN  (Shettle,  private  communication). 

6.  MODTRAN  currently  merges  NOV AM-generated  profiles  (e.g.,  xtinction  and 
temperature)  into  those  describing  the  rest  of  the  atmospheric  profile  fi'om  whatever 
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source  has  been  specified,  'default’  or  'user-defined*.  This  could  lead  to  very  coarse 
discontinuities  whose  impact  might  need  to  be  explored.  General  validation  against 
real  sonde  data  will  provide  additional  confidence  in  the  procedure. 

8.6  Modifications  To  Novam  To  Code 

As  mentioned  NOVAM  modifications  were  kept  to  a  bare  minimum.  Here  is  a  list  of  types  of 
coding  changes  to  NOVAM. 

1 .  All  structure  variables  were  replaced  using  this  scheme: 
stmcture.member  was  replaced  by  structure_member 

This  of  course  meant  that  numerous  corresponding  changes  to  subroutine  arguments  had  to  be  made. 

2.  The  driverS.f  routine  was  substantially  changed  to  output  the  novam.out  file 
described  earlier. 

3.  The  assyml  routine  in  the  file  optics2.f  was  substantially  rewritten  to  fix  an 
interpolation  problem  with  the  asymmetry  parameters. 

4.  The  calls  to  gettim  were  eliminated  as  it  is  not  available  on  all  machines. 

5.  potential_temperature  was  replaced  by  potential_temp  as  this  variable  and  routine 
name  is  too  long. 

6.  The  file  drivesub2.f  was  renamed  drivesb2.f  so  that  the  new  prefix  has  no  more  than 
eight  characters  which  is  the  maximum  for  the  PC  environment. 

7.  As  before  the  sigfile  is  created  by  calling  it  with  repeatflag  equal  to  .false..  In  the 
same  call,  a  new  file  called  invfile  is  created  with  inversion  and  other  extra  layers  to  be 
used  as  MODTRAN  layers.  This  file  also  contains  pressure,  air  temperature  (not 
potential  temperature)  and  RH.  It  is  created  by  modifying  the  routine 
make_rdataary.  Later  the  driver  (with  repeatflag=.true.)  reads  the  invfile  and 
creates  the  novam.out  file  at  these  altitudes. 

8.  The  driver  checks  to  see  that  all  altitudes  in  the  invile  which  are  greater  than  2  km  are 
discarded.  Also  discarded  are  the  set  of  all  top  altitudes  if  the  first  altitude  in  the  set 
has  a  relative  humidity  which  is  below  50%.  That  is  because  the  NOVAM  aerosols 
appear  to  be  restricted  to  be  in  an  environment  of  50%  humidity  or  higher. 
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9.0  MODTRAN  3.7/4.0  AEROSOL  AND  CLOUD  UPGRADES 

MODTRAN  3.7  and  4.0  contain  several  upgrades  to  the  current  aerosol  and  cloud  models: 

The  four  builtin  aerosol  models  are  no  longer  confined  to  fixed  and  essentially  non-overlapping 
regions.  Each  can  be  independently  moved  to  any  region,  and  additionally,  can  be  stretched, 
compressed  and  scaled.  Because  they  can  be  moved  independently,  they  can  occupy  overlapping 
regions. 

The  user-supplied  spectral  parameter  input  schemes  for  aerosols  have  been  improved.  Now  the 
spectral  grid  can  be  arbitrary  (i.e.,  not  limited  to  the  47  fixed  spectral  points)  and  different  for  each 
aerosol.  The  scattering  phase  function  can  have  wavelength  dependence  (not  just  angular 
dependence).  Additionally,  there  can  be  now  up  to  four  (not  just  one)  user-defined  aerosol  profiles. 
These  upgrades,  used  in  conjunction  with  the  stand-alone  Mie  scattering  code  supplied  with 
MODTRAN  3. 7/4.0,  allow  the  user  to  model  the  effects  of  aerosol  more  realistically  than  was 
previously  possible. 

A  stand-alone  Mie  code  was  used  to  recompute  the  spectral  parameters  (extinction,  absorption 
and  the  asymmetry  parameters)  of  the  builtin  clouds  on  a  spectral  grid  consisting  of  788  points 
between  0.2  to  300  microns.  This  is  a  much  needed  improvement  on  the  older  grid  of  47  points  in 
the  same  spectral  interval. 

The  upgrades  were  coded  while  still  maintaining  backward  compatibility  with  earlier 
MODTRAN  input  files  (tapeS’s).  In  the  following  each  upgrade  and  ways  to  access  it  are  described. 

9.1  Moving,  Stretching/Compressing  &  Scaling  Builtin  Aerosols 

The  user  now  has  the  ability  to  move  MODTRAN’ s  builtin  aerosols  from  their  original  positions 
to  arbitrary  altitudes  and  also  to  compress  and  stretch  them  using  two  additional  input  lines.  One 
important  gain  from  this  new  scheme  is  the  ability  to  move  the  tropopause  to  a  more  appropriate 
height  rather  than  having  it  fixed  at  1 0  km.  (Note  that  this  flexible  aerosol  scheme  cannot  be  used  in 
conjunction  withNOVAM,  the  new  Navy  Oceanic  Vertical  Aerosol  Model;  however,NOVAM  can 
be  used  with  other  upgrades.  NOV  AM  refers  to  the  Navy  Oceanic  Vertical  Aerosol  Model,  which  is 
part  of  the  suite  of  aerosol  models  now  in  MODTRAN  3.7  and  higher.) 
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The  new  flexibility  is  invoked  by  having  the  string  (‘A+’)  in  CARD  2  in  columns  1-2  as  shown 
in  the  example  below  (‘A+’  stands  for  AEROSOL  PLUS); 


A+  1  OUSS  0 

0 

0  0 

0.000 

.000 

.000 

.000 

(CARD  2) 

.70 

3.7 

1.0 

15.0 

1.0 

35.5 

65.5 

1.0 

70.0 

90.0 

(CARD  2A+) 

The  string  ‘  AH-’  triggers  reading  a  new  card,  CARD  2A+,  which  actually  consists  of  two  lines  of 
input  and  which  are  inserted  immediately  after  the  current  CARD  2  and  before  CARD  2A  (if  clouds 
are  being  input).  This  new  two-line  card  is  necessary  if  and  only  if  ‘  A+’  is  present. 


CARD  2:  APLUS,  IHAZE,  CNOVAM,  ISEASN,  ARUSS,  IVULCN,  ICSTL,  ICLD,  TVSA, 

VIS,  WSS,  WHH,  RAINRT,  GNDALT 
FORMAT  (A2, 13,  Al,  14,  A3, 12, 3(15),  5F10.5) 

CARD  2A+:  ZAERll,  ZAER12,  SCALEl,  ZAER21,  ZAER22,  SCALE2,  ZAER31,  ZAER32, 

SCALE3,  ZAER41,  ZAER42,  SCALE4  (If  APLUS='A+') 
FORMAT((3(lX,  F9.0),  20X,  3(1X,  F9.0))) 

Only  APLUS,  CNOVAM  and  USS  are  new  to  CARD  2;  all  others  have  the  same  position  and 
meaning  as  before.  CMOVAM  is  a  single  character  immediately  after  IHAZE  and  is  set  to  ‘N’  for 
using  NOV  AM  output  by  MODTRAN.  Although  CNOVAM  and  APLUS  are  both  in  CARD  2, 

NOV AM  aerosols  cannot  be  used  in  conjunction  with  APLUS  aerosol  enhancements. 

USS,  the  other  new  variable  in  CARD  2,  is  a  three-character  string,  for  which  the  only 
menaingful  value  is  ‘USS’  (for  USER-SUPPLIED-SPECTRA).  Previous  to  this  upgrade,  a  user  had 
to  rely  on  builtin  spectral  databases  for  the  builtin  aerosols.  The  builtin  spectra  are  provided  on  a 
very  sparse  spectral  grid  which  is  inadequate  for  many  purposes.  Now  the  user  can  supply  a  denser 
set  of  spectral  data  for  the  builtin  aersols.  The  scheme  for  the  spectral  input  is  essentially  the  same 
as  is  currently  done  using  CARD  2D2.  This  is  expanded  in  more  detail  below.  ' 

There  are  1 2  variables  in  the  two  lines  of  CARD  2A+  as  enumerated  above.  The  first  set  of  three 
is  for  aerosol  number  1 ;  the  second  set  of  three,  for  aerosol  2;  the  third  set,  for  aerosol  3  and  the 
fourth  set,  for  aerosol  4 .  The  meanings  of  the  numerical  values  for  ZAERi  1 ,  Z  AERi2  and  SC  ALEi, 
i=l,  2, 3  and  4,  are  as  follows: 
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ZAERil 


The  base/bottom  of  aerosol  I 


ZAERi2  >  ZAERil 
<  ZAERil 
=  ZAERil 


The  top  of  aerosol  I 

Translate  original  profile  to  new  base,  ZAERil 

Set  values  to  default,  ignore  SCALEi 

(Also  set  to  default  when  both  ZAERil  and  ZAERil  are  blank) 


SCALEi  >0.0  Multiply  vertical  profile  by  SCALEi 

=  0  or  blankMultiply  vertical  profile  by  1 .0  (i.e.,  preserve  column  density) 

The  aerosols  are  linearly  mapped  into  the  new  region  and  the  column  densities  are  preserved  if 
SCALEi  is  unity.  Note  that  since  the  cards  are  read  using  fixed  formats,  blanks  are  interpreted  as 
zeros.  By  default,  SCALEi  is  set  to  unity  if  blanks  or  0.0  are  input.  Note  that  if  the  APLUS  option 
is  used,  the  two  lines  of  CARD  2  A+  must  be  present  even  if  any  of  these  lines  are  intended  to  consist 
of  all  blanks. 

The  oft-used  current  definition  of  an  aerosol  region  leads  to  some  inevitable  confusion.  The 
definition  really  should  refer  to  the  eontiguous  region  where  the  aerosol  concentration  is 
non-negative  and  outside  which  it  is  zero.  By  this  definition,  the  region  of  aerosol  1 ,  for  example,  is 
from  0  to  3  km;  the  profile  linearly  decreases  from  a  non-zero  value  at  2  km  to  zero  at  3  km.  Instead, 
in  MODTRAN  documentation  this  region  is  said  to  be  from  0-2  km.  In  the  MODTRAN  3.7  upgrade, 
the  ZAERil  and  ZAERi2  values  refer  to  the  bounding  altitudes  which  sandwich  the  entire  region 
where  the  aerosol  concentration  is  non-negative.  Table  8  lists  the  default  values  of  these  boimding 
altitudes  along  with  the  commonly  (but  confusingly)  referred-to  region  boundaries  for  each  aerosol. 


Table  8.  Default  Aerosol  Region  Boundaries. 


Aerosol  Region 

Common  Region  Definition 

Actual  ZAERil 

Actual  ZAERi2 

1 

0-2  km 

0  km 

3  km 

2 

2-10 

2 

11 

3 

10-30 

10 

35 

4 

30-100 

30 

100 
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9.2  User-Supplied  Aerosol  Parameter  Upgrades 

This  section  contains  instractions  for  the  MODTRAN  3. 7/4.0  options  that  provide  flexible 
wavelength-dependent  specification  of  extinction,  absorption,  and  asymmetry  parameters  and  phase 
functions.  These  upgrades,  used  in  conjimction  with  the  stand-alone  Mie  code,  allow  aerosols  to  be 
modeled  more  realistically.  The  spectral  grid  can  be  arbitrary  (i.e.,  not  limited  to  the  default  47  fixed 
spectral  points  of  the  previous  MODTRAN)  and  different  for  each  aerosol.  The  scattering  phase 
function  can  have  wavelength  dependence  in  addition  to  angular  dependence.  There  can  be  up  to 
four  user-defined  aerosol  profiles.  In  addition,  utility  programs  are  provided  which  allow 
MODTRAN  to  be  run  with  the  Navy  Oceanic  Vertical  Aerosol  Model  (NOV AM). 

9.2.1,  User-Supplied  Aerosol  Spectral  Parameters  (ARUSS  Option) 

Previous  to  this  upgrade,  the  user  could  provide  extinction,  absorption  and  asymmetry 
parameters  only  for  user-supplied  aerosol  profiles  (IHAZE=7  or  ICLD=1 1)  which  are  in  fact  the 
extinction  values  at  0.55  pm.  Furthermore,  the  spectral  parameters  were  limited  to  the  47 
wavelengths.  This  was  done  using  CARDs  2D,  2D1  and  2D2  with  IHAZE=7  or  ICLD=1 1. 

There  have  been  two  generalizations  to  user-supplied  aerosol  spectral  data: 

•  Now  the  user  can  supply  spectral  data  on  an  arbitrary  grid  for  IHAZE=7  or  ICLD=1 1 .  For 
this  ARUSS  (in  CARD  2)  needs  to  be  set  to  the  three-character  string  'USS'.  Additionally, 
the  meaning  of  the  IREG(N),  N=1 , 2, 3  and  4,  variables  in  CARD  2D  has  been  generalized; 
when  >  1,  they  now  specify  the  number  of  wavelengths  at  which  data  is  supplied. 

•  The  user  can  also  supply  spectral  data  for  the  default  aerosol  profiles,  as  selected  by  IHAZE, 
ISEASN  and  IVULCN  (IHAZE;67  and  ICLD  1),  instead  of  relying  on  the  sparse  built-in 
databases  of  MODTRAN.  This  is  also  done  by  setting  ARUSS  to  the  character  string  'USS'. 
The  USS  option  can  also  be  used  in  conjunction  with  the  APLUS  option. 

The  relevant  cards  for  these  upgrades  are  CARD  2D,  2D1  and  2D2  as  described  below. 

Note  the  units  of  extinction  and  absorption  parameters  in  MODTRAN.  THE  EXTINCTION 
AND  ABSORPTION  COEFFICIENTS  ARE  ACTUALLY  UNITLESS,  BECAUSE  THEY  ARE 
OBTAINED  BY  DIVIDING  THE  ACTUAL  VALUES  BY  EXTINCTION  AT  0.55  pm: 
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Ke„(X)  =  EXT(A)/EXT(0.55  urn);  K^bsW  =  ABS(X)/EXT(0.55  urn). 

CARD  2D:  IREG(l),  IREG(2),  IREG(3),  IREG(4)  (if  IHAZE=7  or  ICLD=11  or 
ARUSS='USS') 

FORMAT(4I5) 

CARD  2D1:  AWCCON,  TITLE  (CARDs  2D1  and  2D2  needed  if  IREG(l)  >  0) 

FORMAT(E10.3, 18A4) 

CARD  2D2:  (VX(1, 1),  EXTC(1, 1),  ABSC(1, 1),  ASYM(1, 1),  1=1,  IREG(l)  or  47) 

(if  ARUSS  is  not  set) 

FORMAT(3(F6.2, 2F7.5,  F6.4)) 

CARD2D1:  AWCCON,  TITLE  (if  IREG(2)  >  0) 

FORMAT(E10.3, 18A4) 

CARD  2D2;  (VX(2, 1),  EXTC(2, 1),  ABSC(2, 1),  ASYM(2, 1),  1=1,  IREG(2)  or  47)  (if 

ARUSS  is  not  set) 

FORMAT(3(F6.2, 2F7.5,  F6.4)) 


CARD2D1:  AWCCON,  TITLE  (ifIREG(3)>  0) 

FORMAT(E10.3, 18A4) 

CARD  2D2:  (VX(3, 1),  EXTC(3, 1),  ABSC(3, 1),  ASYM(3, 1),  1=1,  IREG(3)  or  47)  (if 

ARUSS  is  not  set) 

FORMAT(3(F6.2, 2F7.5,  F6.4)) 

CARD  2D1 :  AWCCON,  TITLE  (if  IREG(4)  >  0) 

FORMAT(E10.3, 18A4) 

CARD  2D2:  (VX(4, 1),  EXTC(4, 1),  ABSC(4, 1),  ASYM(4, 1),  1=1,  IREG(4)  or  47)  (if 

ARUSS  is  not  set) 

FORMAT(3(F6.2, 2F7.5,  F6.4)) 

CARDs  2D1  and  2D2  are  repeated  up  to  four  times,  one  pair  for  each  aerosol.  However,  the  two 
cards  for  aerosol  i  are  needed  if  and  only  if  IREG(N)  >  0.  The  only  differences  between  the  present 

and  prior  forms  are  in  CARD  2D  and  CARD  2D2.  Now  CARD  2D  has  four  integer  values  denoting 
the  number  of  spectral  grid  points  for  each  of  the  four  aerosols;  IREG(N)  =  number  of  spectral  grid 
points  for  aerosol  N.  CARD  2D2  is  the  list  of  the  spectral  parameters:  VX  is  the  wavelength  in 
microns,  EXTC  is  the  extinction  coefficient,  ABSC  is  the  absorption  coefficient  and  ASYM  is  the 


asymmetry  parameter.  Previously  IREGs  were  all  set  to  1  or  0.  A  value  of  1  meant  that  spectral 
parameters  had  to  be  read  using  CARD  2D2  and  the  number  of  spectral  points  were  fixed  at  47. 
(Actually  the  VX  array  was  not  used  at  all  because  the  47  wavelengths  were  already  fixed  in  the  code 
at  an  earlier  point.)  Now,  VX  is  a  2D-array;  the  first  dimension  identifies  the  aerosol  and  the  second 
is  the  wavelength  index.  The  user  must  input  VX  values  in  microns  and  in  increasing  order;  that  is, 
the  first  VX  must  be  the  lowest  wavelength.  The  VX  array  may  differ  for  each  aerosol. 

The  meaning  of  IREG  is  summarized  and  further  clarified  below: 

VALUE  OF  ARUSS  VALUE  AND  MEANING  OF  IREG(i) 

ARUSS-  '  IREG(N)  =  0  No  user-supplied  spectral  data 

IREG(N)  =  1  User-supplied  data  for  the  earlier  set  of  47  fixed  grid 
points.  Although  VX  array  is  read,  they  are  not  used;  instead  the  47 
wavelengths  are  used. 

ARUSS- USS'  IREG(N)  =  0  No  user-supplied  data 

IREG(N)  =  M  User-supplied  data  for  M  arbitrary  wavelengths 

9.2.2.  User-Supplied  Aerosol  Phase  Functions  (CARDS  3bl,  3b2, 3cl-3c6) 

The  user-supplied  phase  function  input  scheme  has  also  been  upgraded.  As  before,  the 
user-supplied  phase  functions  are  read  in  if  IPH  (CARD  3A2)  is  set  to  1.  Now,  the  user-supplied 
phase  functions  can  vary  with  wavelength  in  addition  to  angle.  This  upgrade  is  actually  independent 
of  the  A+  and  USS  upgrades,  and  it  necessitates  a  "generalized"  form  of  CARD  3B 1 : 

CARD  3B1 :  NANGLS,  NWLF  (if  IPH=1 ) 

FORMAT(2(I5)) 

NWLF  is  the  new  variable  which  can  be  either  0  or  a  positive  integer;  0  means  that  the  phase 
function  has  no  wavelength  dependence  whereas  a  positive  integer  means  that  the  phase  function 
will  be  specified  on  a  wavelength  grid  with  that  many  points.  The  phase  function  array,  F,  now  has 
three  indices:  aerosol  index,  angle  index  and  the  wavelength  index. 

If  NWLF  =  0  or  blank,  CARD  3B2  is  used  as  before: 

CARD  3B2:  (ANGF(I),  F(l,  1, 1),  F(2, 1, 1),  F(3, 1, 1),  F(4, 1, 1),  1=1,  NANGLS) 
FORMAT(5E10.3)  (ifIPH=l,NWLF=0) 

If  NWLF  >  0,  CARD  3B2  is  replaced  by  CARDs  3C1-3C6: 
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CARD  3C1:  (ANGF(I),  1=1,  NANGLS) 
FORMAT(8(1X,F9.0)) 


CARD3C2:  (WLF(J),  J=l,  NWLF) 
FORMAT(8(1X,F9.0)) 


CARD3C3:  (F(l,  I,  J),  J=l,  NWLF) 
FORMAT(8(lX,  E9.3)) 

CARD3C4:  (F(2, 1,  J),  J=l,  NWLF) 
FORMAT(8(lX,  E9.3)) 

CARD3C5:  (F(3, 1,  J),  J=l,  NWLF) 
FORMAT(8(lX,E9.3)) 

CARD3C6:  (F(4, 1,  J),  J=l,  NWLF) 
FORMAT(8(lX,  E9.3)) 


(read  angles  if  IPH=1,  NWLF  >  0), 


(read  wavelengths,  in  microns, 
ifIPH=l,NWLF>0) 

(read  phase  function  for  aerosol  1 

if  IPH=1,  NWLF  >  0;  repeat  NANGLS  times) 

(read  for  aerosol  2, 
repeat  NANGLS  times) 

(read  for  aerosol  3, 
repeat  NANGLS  times) 

(read  for  aerosol  4, 
repeat  NANGLS  times) 


In  this  upgrade,  the  wavelength  grid  and  the  angle  grid  is  the  same  for  each  of  the  four  aerosols. 
Furthermore,  the  phase  function  must  be  supplied  either  for  all  aerosols  or  no  aerosol.  For  each,  all 
CARDS  3C3  are  supplied  first,  then  all  CARDs  3C4,  all  3C5,  and  finally  all  3C6;  the  CARDs  for  the 
subsequent  aerosol  then  follow. 


9.2.3.  User-Supplied  Aerosol  Profiles  (CARD  2C3) 

Prior  to  these  upgrades,  the  user  could  only  input  one  aerosol  profile  by  using  the  user-selected 
profile  option,  MODEL=7,  IRD2=1 .  Now  the  user  can  have  up  to  four  user-defined  aerosol  profiles 
with  MODEL=7,  IRD2=2.  (MODEL=0  is  not  allowed.) 

This  upgrade  cannot  be  used  with  the  A+  upgrade  option;  the  APLUS  option  is  ignored  if 
MODEL=7  and  IRD2=2  or  1.  The  A+  option  allows  the  built-in  aerosols  to  be  shifted  around, 
whereas  this  upgrade  allows  the  user  to  input  aerosol  profiles  (up  to  all  four)  with  greater  control. 
The  four  profiles  can  only  be  input  as  altitude-dependent  aerosol  extinction  coefficients  at  0.55  pm. 
Previously  the  single  user-defined  aerosol  profile  could  be  either  the  altitude-dependent  extinction 
coefficient  or  the  altitude-dependent  liquid  water  content  (in  gm/m^).  For  backward  compatibility 
the  previous  option  for  the  single  aerosol  profile  is  maintained. 
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This  upgrade  is  achieved  by  a  generalization  of  CARD  2C3.  For  this  purpose  AHAZE  was 
changed  from  a  scalar  variable  to  an  array,  AHAZE(4).  The  two  versions  of  CARD  2C3  are  shown 
below: 


CARD  2C3:  AHAZE(l),  EQLWCZ,  RRATZ,  IHAl,  ICLDl,  IVULl,  ISEAl,  ICHRl 
(ifIRD2=l) 

FORMAT(10X,  3F10.0, 515) 

CARD2C3:  AHAZE(l),  RRATZ,  AHAZE(2),  AHAZE(3),  AHAZE(4)  (ifIRD2=2) 
FORMAT(10X,  FIO.O,  lOX,  4F10.0) 

The  variables  missing  in  the  newer  version  of  CARD  2C3  (IRD2=2)  are  not  needed  for 
specifying  aerosols.  However,  ICLD 1  (IRD2=  1 )  allows  the  user  to  specify  cloud  profiles  in  addition 
to  aerosol  profiles  with  the  restriction  that  a  cloud  extinction  and  an  aerosol  extinction  cannot  be 
specified  at  the  same  altitude  using  CARD  2C3.  The  price  of  the  current  upgrade  is  the  elimination 
of  the  cloud  extinction  at  an  altitude  for  having  the  luxury  of  inputing  up  to  four  aerosol  extinctions. 
However,  user-specified  cloud  profiles  may  be  entered  using  CARD  2E1. 

As  mentioned,  in  lieu  of  extinction  an  aerosol  profile  could  also  be  given  as  liquid  water  content 
in  gm/m^.  The  conversion  factor  for  converting  liquid  water  content  (gm/m^)  to  extinction 
coefficient  is  given  by  AWCCON.  AWCCON  can  also  be  user-specified  by  using  CARD  2D 
(IHAZE=7  or  ICLD=1 1).  However,  since  in  the  present  upgrade,  the  aerosol  profiles  cannot  be 
stated  in  terms  of  liquid  water  content,  AWCCON  values  in  CARD  2D1  are  not  used. 


9.3  Example  tape5  File 

Here  is  an  example  of  a  tapeS  which  has  both  the  A+  and  ARUSS  aerosol  options.  Notice  the 
CARD  2A+  following  CARD  2  (which  contains  'A+'  as  its  first  two  characters).  Also  note  that 
user-supplied  spectral  data  are  used  for  a  built-in  aerosol  profile. 


M430100000010  -1  .0500 

F  OF  0 

A+  1  lUSS  0000  0.000  0.000  0.000  0.000 

.5  1.0  (CARD  2A+) 

(This  blank  line  must  be  here  or  this  line  should  have  zeros)  (CARD  2A+) 

40  0  0  0  (CARD  2D) 

O.OOOe+OOregion  #1  desert  summer  aerosol 
0.201.016580.434950.8797  0.301.016680.434950.8797  0.341.019360.447350.8857 

0.551.000000.219350.7980  0.691.037040.167430.7666  1.061.114910.037210.7143 
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1.541.208350.043480.7689 
2.500.770220.050250.9116 
3.390.812440.122180.8623 
5.000.608420.104880.8551 
6.200.337160.150810.9065 
7.900.285240.129920.8881 
9.000.829650.339030.6736 
10.590.692100.275960.6814 
12.500.526480.201000.7023 
16.400.466220.177020.6539 
21.300.485200.268970.5959 
40.000.368010.239470.4688 
004.90000  180.00 

2500  02600 


2.001.047120.042120.8557 

2.700.667040.086210.9281 

3.750.788880.100130.8493 

5.500.511680.115510.8706 

6.500.311720.152880.9079 

8.200.301080.188320.8855 

9.200.831530.346750.6684 

11.000.669310.250000.6748 

14.800.493950.190370.6708 

17.200.461220.172020.6408 

25.000.457050.223520.5460 


25  5 


0 


2.250. 905020 . 035770 .8936 
3.000.628860.114680.9247 
4.500.677650.104040.8524 
6.000.362390.150330.9038 
7.200.670350.206630.7748 
8.700.600290.258340.7717 
10.000.808380.344870.6558 
11.500.625310.232900.6805 
15.000.487910.188070.6700 
18.500.442030 . 181610.6366 
30.000.401790.258470.5494 


9.4  Cloud  Spectral  Upgrades 


A  Mie  code  was  used  to  recompute  the  sepectral  parameters  of  the  five  builtin  clouds,  cumulus, 
altostratus,  stratus,  stratus/strato-cumulus  and  nimbostratus,  on  a  dense  grid  of 788  spectral  points 
between  0.2  and  300  microns,  with  the  density  of  the  spectral  points  heavily  weighted  towards  the 
shorter  wavelengths.  This  new  grid  is  a  much  needed  improvement  on  the  older  grid  of  47  points  in 
the  same  wavelength  region.  The  aerosol  optical  properties  were  also  retabulated  at  the  788  grid 
points,  but  the  data  is  simply  an  interpolation  of  the  lower  resolution  data;  in  other  words,  no  new 
Mie  calculation  was  performed  for  the  aerosols. 
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